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The new modulating fuel valve offers flow rates up to 
200 GPM with low power input. This redesigned version 
of Consolidated Controls very successful FLEXBELT® 
valve may be used in control of a wide range of fluids 
and fuels .. . the latest being hydrogen peroxide. 


The flexbelt valve can be provided with linear, square, 
logarithmic or other unusual flow contours to match 
control systems requirements. The complete unit which 
contains a two-phase servo motor, gear train, adjustable 
stops, feedback, potentiometer, flow control belt and 
drum, and seat, weighs just over four pounds. The 
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valve is constructed of stainless steel and the belt 
material for most applications is Teflon impregnated 
glass fibre. Because the belt is self-seating as it rolls 
off the drum, operation is insensitive to contamination. 
This offers users a high degree of reliability in conjun- 
tion with versatility. 


The unique FLEXBELT® principle enables the valve to 
handle a wide variety of flow characteristics with a 
minimum of input power. Write for complete specifica- 
tions to Manager—Sales: Consolidated Controls Cor- 
poration, 750 So. Isis Avenue, Inglewood, California. 


=a CONSOLIDATED CONTROLS CORPORATION 


BETHEL, CONN., INGLEWOOD, CALIF. 
A SUBSIDIARY OF CONSOLIDATED DIESEL ELECTRIC CORPORATION 
STAMFORD, CONNECTICUT 
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This Journal is devoted to the advance- 
ment of astronautics through the dissemina- 
tion of original papers disclosing new scientific 
knowledge and basic applications of such 
knowledge. The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, space flight 
mechanics, high-speed aerodynamics, flight 
guidance, space communications, atmospheric 
and outer space physics, materials and striuc- 
tures, human engineering, overall system 
analysis, and possibly certain other scient fic 
areas. The selection of papers to be printed 
will be governed by the pertinence of the tonic 
to the field of astronautics, by the current or 
probable future significance of the resear:h, 
and by the uparnns of distributing the in- 
formation to the members of the Society «nd 
to the profession at large. 
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The quiet, solid feeling inside 
the 1959 Oldsmobile results from 
patient testing by a highly spe- 
cialized ‘‘jury’’— experts in the 
field of sound analysis. 


Unlike most other phases of auto- 
motive engineering, there are no 
instruments which measure the 
“annoyance” factor of sounds inside 
an automobile. Subjective, personal 
evaluation known as the “jury sys- 
tem” is the most important means of 
determining how quiet it is inside 


an Oldsmobile. 


With the “jury system”, several ex- 
perienced engineers assemble in an 
acoustically treated room and listen 
to binaural tape recordings made 
inside a test car. Careful editing of 
the tape recording allows rapid side- 
by-side comparison between two 
different test conditions and elimi- 
nates reliance on memory since con- 
siderable time may elapse between 
actual tests. By pushing a button, 
each engineer casts a secret “vote” 
and then the tabulation of the results 
determines which test condition is 
superior. 


For the greatest possible accuracy, 
considerable care is taken making 
the recordings. Special microphones 
are placed at points corresponding 
to the passengers’ ears. The recording 
equipment is in a specially equipped 


THE SOUND... 


AND THE JURY 


vehicle which follows the test car 
over a “ride road” that is carefully 
constructed to excite the vibrations 
and sounds that occur under normal 
driving conditions. 


Quality, like quietness, is not acci- 
dental. It is the deliberate effort by 
Oldsmobile engineers to design quality 


into every automobile they produce. 
You'll like the 1959 Oldsmobile 
with its solid Body by Fisher— 
it’s the quietest car built! Experience 
this ae for yourself by visiting 
your local Authorized Oldsmobile 
Quality Dealer as soon as possible. 


OLDSMOBILE DIVISION 
GENERAL MOTORS CORPORATION 


OLDSNMO BILE >-— Where Quality is Standard Equipment! 
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COMMAND DESTRUCT 


The flight testing of second generation missiles—more versatile 
and powerful than their predecessors—requires a device for sure 
termination of any missile flight that might endanger the test 
range or surrounding area. 


Ramo-Wooldridge engineers, under a United States Army Signal 
Corps contract, have successfully developed and delivered the 
first sub-miniature, completely transistorized radio “command 
destruct” receivers. 


Specifically designed for missile flight safety operations, the 
receiver (AN/DRW-11) can actuate safety mechanisms or destruct 
devices. It has three command channels, each of which actuates 
a control relay. 


The “command destruct” receiver accepts frequency modulated 
signals in the UHF radio command control band. It is. designed 
to operate with closer radio frequency and command frequency 
channel spacing than has been used to date, thus making possib]« 
more efficient use of the available radio spectrum. 


Compact and rugged, the radio receiver's modular construction 
permits rapid and complete accessibility to all components. One 
module houses the basic receiver. The second module contains 
the three command channels and relays. This integrated package 
occupies 115 cubic inches, and weighs 4 pounds. The receiver 
requires no pressurization and operates reliably under the adverse 
environmental conditions encountered in missile flight testing. 


Engineers and scientists interested in being associated with some 
of the nation’s most advanced research and development programs 
are invited to acquaint themselves with current opportunities at 
Ramo-Wooldridge. The areas of activity listed below are those 
in which R-W is now engaged and in which openings exist. 


Missile electronics systems 

Advanced radio and wireline communications 
information processing systems 

Electronic language translation 

Anti-submarine warfare 

Air navigation and traffic control 

Analog and digital computers 

infrared systems 

Electronic reconnaissance and countermeasures 
Basic and applied physical research 


For a copy of our brochure, An Introduction to Ramo-Wooldridge, 
or other additional information write to Mr. Donald L. Pyke. 


RAMO-WOOLDRIDGE 


P. O. BOX 90534, AIRPORT STATION * LOS ANGELES 45, CALIFORNIA 


a division of Thompson Ramo Wooldridge Inc. 
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Recent Developments in Convective 
Heat Transfer 


ROLF H. SABERSKY 


California Institute of Technology 
Pasadena, Calif. 


After receiving his Ph.D. degree in 1949, the author joined the staff of the California Institute of 
Technology, where he is now an Associate Professor of Mechanical Engineering. He has been active 
in teaching and research in the field of heat transfer and thermodynamics, and has written a text- 
book on this latter subject. Rocket propulsion has been one of the author’s special interests. 
He has been associated with the Aerojet-General Corp. since 1943 and is now a consultant for that 
company. Dr. Sabersky is a member of the American Rocket Society. 


N WRITING about recent developments in convective heat 

transfer, it must first be admitted that the word “recent’’ 
lias not been interpreted very literally. The developments 
discussed have been selected largely because they are of cur- 
rent interest, and the exact time at which they were achieved 
lias been considered of secondary importance. An attempt 
has also been made to include enough of the background of 
cach topic so that the account may be of use to the reader 
whose main interest lies outside the field of heat transfer. 
Furthermore, there is no pretense that this short paper should 
give a comprehensive survey of the extensive field of convec- 
tive heat transfer. On the contrary, only a few specialized 
subjects have been selected, and in making this selection the 
writer’s own interests and prejudices could not, of course, be 
fully suppressed. Heat transfer involving high speed flow has 
been left out completely, as it has been taken up in other 
articles of this series. The topics which are discussed in this 
review are transpiration and film cooling, heat transfer in 
rough tubes, nucleate boiling, heat transfer to fluids near the 
critical point, and fiow with internal heat generation. Ref- 
erence has been made to some of the principal sources of in- 
formation. Although the given bibliography is relatively ex- 
tensive, it cannot be considered complete, and the interested 
reader may wish to use the cited references as a further guide 
to the literature. 


Film and Transpiration Cooling 


Let us consider a surface being heated on one side by a 
stream of hot gases. The principle of transpiration cooling 
consists of cooling such a surface by injecting the coolant 
directly through the surface and into the hot stream. The 
coolant is to form a protective layer on the hot side. The in- 
jected fluid acts not only by absorbing heat but also by alter- 
ing the flow pattern so as to reduce the heat transfer rate. 
In order to be able to cool by transpiration, the wall has to be 
porous. Attempts were made originally to produce porous 
walls by sintering powdered metals (1).1. This method 
proved to be feasible, and even such complicated shapes as 
rocket nozzles could be made. The porosity was controlled 
by mixing the metal powder with a filler of a definite grain 
size. The filler would then evaporate during sintering, leav- 
ing the metal with the desired porosity. One of the major dif- 
ficulties was probably that related to the clogging of the 
pores. More recently, the manufacture of porous surfaces has 
been advanced considerably by the invention of a new tech- 


Received Feb. 1958. 
' Numbers in parentheses indicate References at end of paper. 
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nique (2) in which a metal wire is taken as the basic compo- 
nent, and the wire is wound to yield a sheet of the desired 
porosity. The sheet is then sintered. The final product is 
generally stronger than the structure obtainable from pow- 
dered metals. It is relatively easy to produce shapes of vari- 
ous kinds, and by using different winding patterns variations 
in porosity and other special properties can be obtained. 
Early experiments on transpiration cooling have demonstrated 
the feasibility of the method (1). These experiments showed 
that relatively small rates of transpiration cooling would pro- 
tect a surface requiring high heat transfer rates when cooled 
convectively. 

Film cooling differs from transpiration cooling in that the 
coolant is injected through individual orifices spaced relatively 
far apart. The coolant is injected in such a way that a pro- 
tective film is formed adjacent to the wall. Film cooling has 
the advantage of not requiring any special porous materials. 
On the other hand, the distribution of the protective film may 
be less effective. One should actually divide film cooling into 
two separate types, depending on whether the coolant is an 
evaporating liquid or a gas. If the coolant is a liquid, the 
evaporating vapor leaves essentially in a direction normal to 
the surface, and the mass addition is therefore quite similar 
to that in transpiration cooling. If a gaseous film is injected 
tangentially to the wall, the influence on the boundary 
layer formation is different. Film cooling with a liquid cool- 
ant is very advantageous because of the benefit which is 
derived from the latent heat. When reference is made to film 
cooling, it is often this type of cooling which is being con- 
sidered. 

One of the best known applications of film cooling has been 
in the thrust chamber of the V-2 missile (3). Several sets of 
film cooling holes were provided in the thrust chamber of this 
unit for the protection of critical areas, such as the throat. 
In film cooling with a liquid coolant, at least two problems in 
addition to those occurring in transpiration cooling have to be 
considered. The first of these is the problem of injecting the 
coolant in such a way that it attaches to the surface in order 
to form a film. The second is that of determining the area 
which can be covered from a single injection station. 

Extensive experiments have been conducted regarding this 
latter aspect (4). Ina particular set of experiments, the main 
flow consisted of heated air flowing through a duct. The walls 
of the duct were protected by film cooling with water. The 
film cooling rate was varied, and the corresponding change 
in the protected area was observed. The tests were repeated 
for several air flow rates. In each case it was seen that for 
low film cooling rates the protected area was directly propor- 
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tional to the flow rate. Beyond a certain critical value, a 
progressively more rapid increase in flow rate was required in 
order to extend the protected area. Visual observation 
showed that the initial thickness of the film increased as the 
flow rate was increased. At a flow rate close to the critical 
one, the film thickness seemed to reach a value at which the 
surface became unstable. Waves appeared on the surface, and 
fluid particles seemed to be torn from the crest of the waves 
and carried into the main stream. Such particles do not 
contribute to the cooling, and it is because of their breaking 
away that an ever-increasing flow rate per unit of protected 
area is required beyond the critical point. In general, there- 
fore, it will be desirable to cool in such a way that the film re- 
mains stable. This will mean that from each injection point 
a relatively small area only can be covered, and in order to 
protect an extended surface several injection stations may be 
required. 

Attempts have been made to analyze the unstable behavior 
of the liquid film (5). In this connection, it was noted that 
the critical thickness y; for all of the above mentioned experi- 
ments corresponded to a value of 25 for the dimensionless 
thickness y,*, where 


n* = 
v 


and where 


T = wall shear 
p = the density of the film 
v kinematic viscosity of the film 


These data led to the hypothesis that the value y,* = 25 
represented a universal constant for the stability of the film. 
This belief was strengthened by the fact that in the flow of a 
fluid in a pipe, the edge of the turbulent boundary layer also 
occurs for a value of y;* equal to about 25. It appeared, there- 
fore, that the generation of turbulent motion in the flow of a 
single fluid and the onset of instability in a film might both be 
based on the same mechanism. Further experiments (6), 
however, showed a dependence of the dimensionless critical 
thickness on the ratio of the viscosity of the film to that of the 
main stream. It must, therefore, be concluded that the critical 
film thickness cannot be given by a single value of y,*. 
Nevertheless, it may be useful, in the further development 
of the theory of film stability, to consider possible connec- 
tions between the factors limiting the thickness of a stable 
film and those determining the distance to the turbulent core 
in normal pipe flow. 

The problem of attachment of a liquid film has also been 
studied experimentally. In addition to the properties and 
velocities of the fluids, the size as well as the direction 
of the injection orifices must be counted as_ variables. 
For injection normal to the direction of flow of the main 
stream, it was found that the Reynolds number of the gas and 
the ratio of the specific kinetic energy of the main stream to 
that of the injected fluid were important parameters (6, 7). 
In addition, the experiments of one author show that the 
Reynolds number of the liquid as well as a cavitation parame- 
ter (6) are also quantities affecting the attachment of a film. 
The cavitation parameter is defined as 


Pr) 


p, = density of the liquid 

V = average liquid velocity 

p = prevailing pressure 

Pp,» = vapor pressure of the liquid at the given temperature 


The occurrence of such a parameter is certainly plausible and 
is explained by the following reasoning. The cooling fluid in 
order to form a surface film has to turn sharply immediately 
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upon leaving the injection orifice. The pressure at the 
orifice edge must, therefore, be very low. When this pres- 
sure reaches the vapor pressure, detachment of the jet is 
likely. Experiments were also conducted in which the injec- 
tion orifices were inclined at various angles towards the direc- 
tion of the main flow (7). It was shown that the condition for 
attachment became considerably more favorable as_ the 
centerline of the orifice was inclined toward the direction oi 
the main flow. In fact, in the set of experiments under dis- 
cussion, no detachment could be observed when the angle be- 
tween the centerline of the orifice and the surface was les- 
than 45 deg. 

The problem of heat transfer from a stable evaporating filn: 
is essentially the same as that from a fluid layer injected 
through a porous wall. The most important application: 
probably involve turbulent flow of the main stream. Severa! 
theoretical treatments of this problem have been offered (see 
e.g., (8, 9)). In one of the earliest, the velocity and tempera- 
ture profiles were determined for the injection of a gaseous 
coolant into a stream in turbulent motion (8). In this 
analysis the flow was divided into a laminar sublayer and « 
turbulent core. The temperature and velocity profiles in the 
laminar layer were calculated taking into account the rate ot 
injection of fluid through the wall. In the turbulent core the 
Reynolds analogy was assumed to be applicable. Further- 
more, the flow was considered to remain unaffected by the 
mass injection, and the usual relations for the friction co- 
efficient, the velocity profile, and the thickness of the laminar 
layer were assumed to hold. These assumptions were made 
in order to obtain a first approximation, and they will be 
modified as additional information becomes available. 
Nevertheless, the predicted profiles were in relatively good 
agreement with measured values (10). The theory was next 
extended so as to be applicable to an evaporating coolant as 
well as to a coolant which might undergo chemical reaction 
in the boundary layer (11). In a subsequent analysis for 
evaporating coolants (12), the Reynolds analogy was used in 
a somewhat different form in an attempt to take into account 
the effects of mass addition on the flow in the turbulent core. 
In addition, the surface temperature of the liquid was made 
to correspond to the local vapor pressure. The results of this 
analysis have been compared with experimental data and en- 
couraging agreement has been found to exist. 

Further studies of film cooling with both laminar and tur- 
bulent flow have been made (see, e.g., (13, 14, 15 and 16)). 
Many of these also include the effects of compressibility. In 
(15 and 16), which are concerned with turbulent flow, par- 
ticular care has been taken to account for the effect of mass 
injection on the flow in the turbulent core. In the develop- 
ment, assumptions are made regarding the dependence of the 
turbulent exchange coefficients on the velocity profile and the 
distance from the wall. Most of the earlier work was limited 
to the injection of a gas with properties identical to those of 
the main flow. Several sets of experimental data (see, e.g., 
(17, 18, 19 and 20)) have also been obtained for this case, and 
they are in fair agreement with the predicted values. Re- 
cently, emphasis has been given to the problem in which the 
coolant gas differs from that of the main flow (21, 22). The 
computations show the advantages that can be gained from 
the injection of a gas of low molecular weight. 

One other paper should be mentioned here in which an 
analysis is given of the heat transfer in laminar flow to a non- 
porous surface, just downstream of a transpiration-cooled 
section (23). Such a combination of porous and nonporous 
regions may be desirable for structural reasons, and the 
analysis was performed in order to investigate the effective- 
ness of such an arrangement. In the limit, as the porous re- 
gions are reduced, the conditions for gaseous film cooling are 
approached in which the coolant is injected through individual 
slots. A special analysis for this type of film cooling has also 
been performed (24). This analysis is based on an ingenious 
simplification by which the injection slots are regarded as heat 
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sinks. The effect of injection on the hydrodynamic boundary 
layer is neglected. The results of the analysis have been com- 
pared (25) to a set of older experimental data (26), and satis- 
factory agreement is shown to exist. Additional experi- 
mental results have become available in recent years (see, e.g., 
(28, 29)). 

The several analyses mentioned are believed to show fruit- 
ful approaches to the problem of transpiration and film cool- 
ing. They are, however, based on many assumptions which 
have not yet been sufficiently verified in detail. Special ex- 
periments with the purpose of investigating these assumptions 
will probably be required before any significant further ad- 
vances may be made. 

Before concluding this subject, the unique character of film 
and transpiration cooling should be emphasized. It is the 
only method by which a wall can be kept at very low surface 
temperatures, independent of the heating rate or the tempera- 
ture of the hot gas. This is quite different from ordinary con- 
vective cooling, in which the surface temperature is a function 
o: the hot gas temperature, and depends on the relative magni- 
tude of the heat transfer coefficients of the hot gas and the 
coolant, and on the heat transfer resistance of the wall. Film 
or transpiration cooling does, of course, require the expendi- 
tire of mass. From the investigations completed so far, how- 
ever, it can already be shown that this expenditure is feasible 
for many applications (see, e.g., (30)). For large rocket mo- 
tors, for example, the amount of coolant required can be cal- 
culated to be of the order of 1 per cent of the propellant flow, 
assuming that the coolant does not benefit the thrust. Actu- 
ally, however, the coolant does add to the thrust, particularly 
if the coolant is one of the propellants. In that case, film 
cooling with unstable films may also be indicated. The fluid 
particles which are torn off by the main stream can no longer 
be counted as lost, since they will probably react just as the 
remainder of the propellant. Film cooling with a small num- 
ber of injection stations then becomes possible, which is most 
desirable from the point of view of design and control. It may 
well be that film cooling will become an essential factor in the 
design of rocket motors using high energy propellants. Fur- 
ther study of this method of cooling should be of considerable 
interest. 


Heat Transfer in Rough Tubes 


Friction in rough tubes has been considered to be one of the 
fundamental problems in fluid mechanics, and the problem 
has been investigated thoroughly. The most extensive ex- 
perimental results are probably those due to Nikuradse (31). 
One would think that heat transfer in rough tubes should be 
considered an equally basic problem. Nevertheless, it is ap- 
parent that heat transfer in rough tubes has not received the 
same attention as the corresponding friction problem. Some 
studies have, of course, been devoted to this subject in the 
past, and recently it has received somewhat more emphasis. 
The studies in which friction and heat transfer are treated 
jointly are of particular interest, as these two phenomena 
are in general intimately connected. 

In one of the earliest studies (32), experimental data were 
obtained for water flowing in tubes which had been roughened 
by knurling. These data indicated that there was some in- 
crease in the heat transfer coefficient with roughness. The 
friction coefficient, however, showed a much larger increase 
than the heat transfer coefficient. As a consequence, the 
ratio of the heat transfer coefficient to the friction coefficient 
decreased with increased roughness. Although there has been 
some question regarding the accuracy of these data, the gen- 
eral results have been reproduced by subsequent investiga- 
tors. 

In one of the important recent studies (33), a set of experi- 
ments with rough tubes has been performed, using air as the 
working fluid. The roughness in these tubes was produced by 
inserting specially shaped piston rings as roughness elements. 
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These results confirmed that there is a tendency for the fric- 
tion coefficient to increase more rapidly than the heat transfer 
coefficient. In addition, it was shown that a fairly unique 
relationship seemed to exist between the two coefficients, and 
that this relationship was essentially independent of the type 
of roughness element. A few tests were also conducted with 
a naturally rough tube, and the behavior of this tube did not 
seem to differ materially from that of an artificially rough- 
ened tube. In particular, the results of the tests on the 
naturally rough tube were compared to those of an artificially 
roughened tube having the same friction factor. It was shown 
that both tubes also vielded the same heat transfer co- 
efficient. 

There are, however, other experiments which may have an 
important bearing on the above results (34, 35). In these ex- 
periments, flow in annular ducts was studied. The inner 
surfaces of the ducts were provided with transverse fins in 
the form of annular disks oriented normal to the flow direc- 
tion. Various fin sizes and spacings were examined, the fins 
in some cases extending into the duct for about 40 per cent of 
the width of the annulus. Heat was transferred from the 
inner surface to the fluid, and the heat transfer rate and 
friction drop were measured. It was found that the presence 
of the disks would increase the heat transfer coefficient as well 
as the friction coefficient. In (34), it was shown that for most 
spacings the ratio of heat transfer to the friction coefficient 
would decrease below the value corresponding to a smooth 
wall, but that for certain particular spacings, equally high or 
even higher ratios seemed obtainable. The exact numerical 
values in this comparison are somewhat in question, because 
the values for the coefficients for smooth walls were estimated 
rather than measured. There is also some difficulty in making 
analytical comparisons between heat transfer and the friction 
coefficient. The latter has to be derived from pressure drop 
measurements, and it therefore constitutes an average be- 
tween the friction coefficient of the smooth outer wall and 
that of the rough inner surface. Nevertheless, for engineering 
purposes the above result is most significant; it indicates 
that by using a properly roughened surface a given amount 
of heat may be transferred from a smaller length of duct with 
a possible decrease in pumping power. The aim of the ex- 
periments in (34) was actually to produce such an effective 
surface for the heat transfer ducts of the nuclear reactors in 
Calder Hall. Experiments were also performed in order 
to visually examine the flow pattern for the various disk ar- 
rangements (34, 35). For the spacing most favorable to heat 
transfer, a set of stable vortices would develop between adja- 
cent disks. The velocities due to these vortices were ap- 
parently responsible for providing good heat transfer condi- 
tions along all exposed surfaces. For less favorable spacings, 
one might expect relatively quiescent wake regions behind 
each disk. In this case, the heat transfer from the surface 
bounding these. regions might be poor, but the friction effect 
might remain high because of the drag forces caused by the 
disks. 

It must be admitted, of course, that the surface with annu- 
lar disks as described above differs from the usual concept of 
a rough surface. Nevertheless, for the present purpose one 
may well think of it as an artificially roughened surface, and 
one may suspect that ordinary rough surfaces might show a 
similar behavior. It seems, therefore, justified to question if 
the shape and distribution of roughness elements will in- 
fluence the ratio of the heat transfer to the friction coefficient. 
The set of experiments previously cited (33) may possibly not 
have covered a sufficient range of surface arrangements to in- 
dicate such an effect. Further work on this aspect promises 
to be most interesting as well as of practical importance. 

Three additional papers should be mentioned as they also 
contain results of experiments which may be helpful in for- 
malizing a theory for the heat transfer in rough tubes. The 
first contains data for the flow of air through artificially 
roughened tubes (37). The second presents information for 
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the flow of air through a large number of different, naturally 
rough tubes (38). In the third, a set of experiments is re- 
ported concerning the flow of water through a knurled tube 
(39). In these tests the effect of Prandtl number was studied. 
The change in Prandtl number was obtained by varying 
the water temperature. The Prandtl number is one of the 
essential parameters in heat transfer, and theories can often 
be tested by comparing their predictions as to the effect of 
Prandtl number with corresponding experimental results. 

Attempts to analyze the effect of roughness on heat trans- 
fer have also been made. In an appendix to his well- 
known paper on heat transfer and fluid friction, Martinelli 
(40) suggests a simple extension to the relationship which 
applies to smooth tubes. The modified equation fits some of 
the experimental data surprisingly well. In a more recent 
study, an interesting model for analyzing the effect of rough- 
ness was proposed (33). In this proposal, the drag force 
offered by the rough wall is divided into surface drag and form 
drag. A wire grid is then imagined to be suspended in the 
fluid at the edge of the laminar boundary layer. The drag 
of the grid represents the additional shear caused by the 
roughness. The flow is then analyzed on this basis, and the 
results compare favorably with the experimental data ob- 
tained with air (see, e.g., (33, 41)). The model is an interesting 
one and certainly deserves further study. Some of the as- 
sumptions, such as the one that the thickness of the laminar 
boundary layer is unaffected by roughness, will have to be 
checked. Further modification of the model may be required 
if the roughness type should prove to have an effect on the 
relationship between friction and heat transfer. 


Nucleate Boiling’ 


Boiling heat transfer occurs when heat is transferred from a 
surface to an adjacent liquid under such conditions that the 
vapor phase is formed in the neighborhood of the surface. 
Boiling heat transfer occurs in many applications, such as 
the cooling of rocket motors and the transfer of heat from the 
core of a nuclear reactor. The phenomena occurring in 
boiling heat transfer have been subject to intensive study. 
Qualitatively, the process is now rather well understood, and 
the typical experimental results have become quite familiar. 
These results are usually presented in the form of a graph of 
the heat transfer rate vs. the temperature of the heating sur- 
face at constant fluid temperatures. As the surface tem- 
perature is increased above that of the liquid, heat transfer 
will take place by the usual convection process. When the 
temperature reaches a value approximately equal to the 
boiling point of the liquid at the prevailing pressure, bubbles 
appear in the liquid. If the liquid is near its saturation tem- 
perature, the bubbles will detach. If the liquid is at a tem- 
perature well below its boiling point, most bubbles will grow 
to a certain size and then collapse again at the surface. As 
the bubbles appear, a very marked increase in the heat trans- 
fer coefficient is observed, so that a minor increase in tem- 
perature brings about a very large increase in the heat trans- 
fer rate. As the surface temperature is raised further, the 
bubbles become more numerous and the strong effect of 
temperature on heat transfer continues. A maximum heat 
transfer rate is reached just before the bubbles become so num- 
erous that they coalesce periodically, forming an unstable and 
intermittent vapor blanket over the surface. An additional 
increase in surface temperature leads to more vapor formation 
and to a decrease in the heat transfer coefficient. As a 
consequence, the heat transfer rate now decreases with 
increasing temperature. Finally, sufficient vapor is gener- 
ated so that a rather stable vapor film forms over the sur- 
face. Further increase: in surface temperature will now again 


2? For an account of work in boiling heat transfer in Russia, 
the reader is referred to a book by 8S. 8. Kutadeladze which has 
been translated into English (42). 
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lead to an increased heat transfer rate. The region from 
the initiation of the first vapor bubbles to the point of max- 
imum heat transfer rate is known as the nucleate boiling region. 
The region in which the heat transfer rate decreases with 
surface temperature is called the partial film boiling region, 
and the region beyond has been named the stable film boiling 
region. The fact that a maximum occurs at the end of the 
nucleate boiling region is a characteristic of boiling heat 
transfer. In order to exceed this maximum, the surface has 
to pass into film boiling. The corresponding surface tem- 
perature may, however, have to be very high, the exact 
magnitude depending on the properties of the fluid involved 
In many instances this temperature would have to be s 
high that the surface material would melt if attempts wer 
made to exceed the maximum heat transfer rate of the nu- 
cleate boiling region. This maximum is, therefore, often 
called the burnout point. This general description o: 
boiling heat transfer applies regardless of whether the fluic 
is at its saturation temperature or below, and it also applies 
to free convection as well as to forced convection. However 
the existence of a maximum and the transition from nucleate 
boiling to film boiling becomes, of course, less pronounced as 
the pressure approaches the critical pressure of the fluid. 
A large amount of data for various liquids under a variety of 
conditions has been obtained (see, e.g., (43 to 48)). A very 
complete list of information available on rocket propellants 
is contained in (49). 

The nucleate boiling region is of special importance in 
engineering, since it affords a possibility of obtaining very 
high heat transfer rates without necessitating excessive wall 
temperatures. The burnout point takes on a special sig- 
nificance since it presents a limit which in many designs should 
not be exceeded. In analyzing nucleate boiling it is con- 
venient to divide the problem into three parts: The gen- 
eration of the bubbles, the growth and collapse of the bubbles, 
and the interaction between bubble motion and heat trans- 
fer. As to the generation of bubbles, most investigators 
believe at present that bubbles form from pre-existing nuclei 
(see, e.g., (50, 51, 52)). These nuclei probably consist of gas 
or vapor pockets which are stabilized either in cavities on the 
surface or on solid impurities contained in the liquid itself. 
As the temperature of the surrounding fluid increases, vapor 
will evaporate into these nuclei. If the inside pressure be- 
comes sufficiently large to overcome the effects of surface 
tension, the nucleus will grow into a bubble. The next 
phase, i.e., the growth and subsequent detachment or col- 
lapse, is governed largely by the heat transfer between the 
bubble and the surrounding fluid. Finally, the improvement 
in the heat transfer coefficient is believed to be due in part to 
the agitation caused by the bubbles, and in part to the vapor 
transport. In the case where bubbles are actually leaving 
the surface, the effect of vapor transport is a most important 
one. If the bubbles collapse on the surface, most inves- 
tigators believe the effect of agitation to be the essential one, 
although some feel that the rate of evaporation across the 
bubble makes a significant contribution to the total heat 
transfer. 

This general description of nucleate boiling probably still 
applies if the fluid flowing along a pipe is made to spiral. 
Quantitatively, however, the results have been shown to 
differ. Experiments have recently been carried out in which 
water was ducted through a heated tube in such a manner 
that it had a tangential velocity component in addition to the 
usual axial one (53). In this way, a very significant increase 
in the burnout point was achieved. The increase was larger 
than that corresponding to an axial velocity equal to the total 
velocity along the spiral path, and the change could not be 
explained by the increase in surface velocity. The spiral 
motion does, however, give rise to a strong radial pressure 
gradient, and it is believed that the influence of this gradient 
on bubble motion is responsible for the observed increase in the 
burnout heat flux. The fact that spiral motion may lead 
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to an increase in the burnout point may be of great value to 
the designer of heat transfer equipment. 

It is, of course, the purpose of a heat transfer investigation 
to arrive at a sufficient understanding of the phenomenon so 
that the heat transfer rate can be predicted quantitatively as a 
function of the fluid properties and the wall temperature. 
In working towards this goal, several of the individual steps 
in the process of nucleate boiling have been studied in detail. 
The bubble growth has, for example, been analyzed quite 
extensively (see, e.g., (54,55,56)). Most of these analyses have 
been based on a spherical bubble and spherically symmetrical 
conditions. The results are in satisfactory agreement with 
experiments which were devised to fulfill the assumptions of 
the analyses (57). In recent years, an interesting new ap- 
proach has also been pursued which is aimed particularly at 
determining the maximum heat transfer rate which can be 
achieved in nucleate boiling. In this approach, attention is 
focused on the vapor in the region of partial film boiling. 
The stability criteria of this film are studied, and the maxi- 
mum of the nucleate region is regarded as the limit of the 
partial film boiling region (58,59). The concept of film 
stability in connection with the boiling heat transfer may 
prove to be very fruitful. 

A complete solution of a problem as complicated as nu- 
cleate boiling heat transfer will be very difficult to obtain. 
As in many other problems in fluid mechanics or heat transfer, 
nowever, the determination and experimental correlation of 
the significant dimensionless parameters would in itself be of 
great usefulness. Considerable attention has been devoted 
to this approach. The problem is, however, so complicated 
that even this approach proved to be by no means easy. 
The reader who is interested in educational problems might 
use this fact as an example to illustrate that the application 
of dimensionless analysis and the Pi-theorem does not lead 
to the solution of a problem quite as readily as some textbooks 
would have one believe. 

Two principal types of correlations have been developed: 
Correlations for the nucleate boiling region and correlations 
for the point of maximum heat transfer in nucleate boiling. 
For many engineering applications, the latter are of particular 
importance. In selecting the significant parameters, careful 
attention was paid to the theoretical analyses described pre- 
viously. In developing appropriate expressions, many in- 
vestigators followed the concept that the heat transfer is due 
to vapor transport and added agitation. In order to express 
the effect of agitation, appropriate Reynolds numbers were 
selected (see, e.g., (51, 60, 61)). In forming these Reynolds 
numbers, the bubble diameter and bubble growth velocities 
(change of radius with time) have often been used as the sig- 
nificant quantities. A great simplification was introduced 
when it was shown that for a bubble growing in a fluid of 
uniform temperature, the product of growth velocity and 
diameter was a constant (61). Although the temperature 
field is not generally uniform, a correlation for the nucleate 
boiling region based on this simplification (61) has shown 
promising agreement with experiments. It has been shown to 
cover the results of fluids with such wide variations in prop- 
erties as water and liquid metals (62). 

A comprehensive equation for the determination of the 
point of maximum heat transfer has also been published 
recently (60). The equation is intended to be applicable 
when the motion of the liquid is forced as well as when it is 
free, and it is applicable to any fluid at all temperatures in- 
cluding the saturation temperature. Correlations based on 
the concept of film stability have also been developed. It is 
interesting and encouraging to see that one of these is quite 
similar in form to the correlation based on the concept of 
a Reynolds number (59). 

It appears, therefore, that the attempts to develop 
suitable correlations for nucleate boiling are becoming more 
successful, and there is some hope that further experimenta- 
tion will prove some of the proposed relationships to be ade- 
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quate for many engineering applications. 

There is, however, at least one phase which will need further 
intensive study before the present correlations can be used 
with confidence and that is the nucleation process which was 
briefly described previously. It was stated that the initial 
appearance of bubbles was believed to depend on the size and 
number of pre-existing gas or vapor cavities. This number 
will depend on the condition of the surface, the surface ten- 
sion and contact angle between the fluid and the surface, 
and on the amount of dissolved gas and impurities contained 
in the liquid. The fact that these various characteristics 
can have an influence on nucleate boiling has been demon- 
strated experimentally (see, e.g., (51, 52, 63, 64)). In partic- 
ular, these factors may influence the temperature at which 
bubbles will first form, the slope of the curve of heat transfer 
vs. wall temperature, as well as the wall temperature at which 
burnout occurs. 

There is some indication, however, that nucleation is suf- 
ficiently insensitive to the above mentioned surface factors, 
so that variations in these factors can be ignored for “normal” 
conditions and when extreme accuracy is not required. 
There is also some indication, based on the reasoning pre- 
sented in (59) and on the form of the correlations in (59 and 
60), that the burnout heat flux is relatively independent of 
surface conditions. If these indications can be verified in 
general, correlations for nucleate boiling may prove successful. 

In discussing nucleate boiling, it has so far been assumed 
that the main body of the fluid was in liquid form. The case 
in which the main body itself is a mixture of liquid and vapor 
is becoming of increasing importance. Heat transfer to such 
a two-phase mixture is actually not a new problem. On the 
contrary, it occurs in the boiler of every conventional steam 
powerplant. In boilers, however, the heat transfer from the 
combustion gases to the steam is largely controlled by the gas 
film resistance which is usually the highest resistance in this 
process. The exact magnitude of the heat transfer coefficient 
between the pipe wall and the steam is, therefore, often of 
seconds ry importance, and the heat transfer rates are gener- 
ally sufficiently low so that there is no danger of “burnout.” 
In water-boiling nuclear reactors, however, the gas film re- 
sistance is absent, and the transfer process between wall 
and steam becomes of major significance. It is largely be- 
cause of this application that the heat transfer to vapor- 
liquid mixtures is being studied extensively. 

Studies of the flow of two-phase mixtures have been under- 
taken previously (see, e.g., (65, 66)). In some of these it was 
shown that different types of flow are possible, each type 
corresponding to a certain distribution of vapor and liquid. 
In the region of most interest in heat transfer, the distribution 
is probably such that there exists a vapor core surrounded by 
a liquid film. It is quite possible, therefore, that the heat 
transfer between the wall and the film is similar to that pre- 
viously described for nucleate boiling. 

Several correlations have been proposed based on experi- 
mental data. In particular, one of the correlations for the 
burnout point, previously mentioned in connection with 
nucleate boiling in a liquid (60), has also been suitable in 
correlating some of the data for liquid-vapor mixtures. This 
fact, as well as the occurrence of a burnout point in itself, may 
indicate that the process of heat transfer consists of nucleate 
boiling in a liquid film as proposed above. Other correlations, 
limited to water-steam mixtures, have also been suggested. 
In one of these, the thickness of the assumed liquid film ap- 
pears as a factor (67). The thickness is calculated on the 
basis of the quality of the flowing mixture. Another cor- 
relation involves the total enthalpy (68) of the mixture as a 
significant quantity, which again is a function of the quality. 

From the foregoing account on boiling heat transfer, it may 
be concluded that significant progress has been made in this 
field. The phenomenon is, however, a very intricate one, 
and many properties of the fluid as well as of the heating sur- 
face seem to exert an important influence. It therefore 
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must be expected that a considerable amount of time will 
pass before the engineer will have available a reliable general 
correlation suitable for all fluids. Until then, he will have to 
base his design largely on direct experimental data obtained 
for the particular fluid of interest and under conditions simu- 
lating those of the prototype as closely as possible. 


Heat Transfer to Fluids Near the Critical Point 


In several modern applications, heat transfer takes place to 
fluids which are near their critical state. Steam powerplants, 
for example, are now being designed for pressures of 5000 psi 
and higher, and in the process of heating, the steam passes 
through states which are relatively close to the critical point 
(p., = 3206.2 psia, 7... = 705.4 F). Asa second example, the 
cooling of rocket motors may be cited. Fuels based on hydro- 
carbons have critical pressures of the order of 500 to 600 psia, 
and liquid oxygen has a critical pressure of 730 psia. When 
these fluids are used as coolants, they will in many designs 
reach pressures and temperatures close to the critical values. 

Some experimental work on the heat transfer to fluids near 
their critical point is available. The volume of existing data 
is, however, not yet very extensive. Experiments which 
have been performed include work with water (69), oxygen 
(70), nitrogen (71) and carbon dioxide (72). 

From one point of view, heat transfer to fluids near the 
critical point may not require a special classification. The 
principal characteristic of fluids in this region is that their 
properties are strongly affected by pressure and temperature. 
It may, therefore, be that the heat transfer to these fluids 
may be regarded as normal convection to a fluid under con- 
ditions for which the variations of properties through the 
boundary layer are of essential importance. From another 
point of view, however, heat transfer near the critical point 
may be regarded as the limit of boiling heat transfer, and one 
might suggest that some of the characteristics of this type of 
heat transfer may still prevail. 

To study these aspects of the problem, heat transfer co- 
efficients were calculated on the basis of analyses in which 
the change in fluid properties throughout the boundary layer 
was taken into account. The analytical results were com- 
pared to the experimental data in question, and it was seen 
that the trends were in general well predicted by the theory 
(69, 72). However, there was at least one exception, in that 
the heat transfer coefficient for water in the immediate neigh- 
borhood of the critical temperature showed a very marked 
increase not predicted by the theory (69). On the basis of 
this result, it was suggested that the fluid in this region still 
might behave similarly to a fluid in nucleate boiling (see, e.g., 
(69, 73)). Since the pressure is above the critical value, 
one can not, of course, expect the appearance of a second 
phase. If the decrease in density with temperature is very 
pronounced, one may, however, imagine that regions of rela- 
tively low density form near the heating surface. These 
regions might become unstable and grow and collapse or 
detach from the surface. In other words, they may possibly 
behave in a manner very similar to bubbles in nucleate boiling. 
If so, the activity of these low density regions may provide 
the explanation for the improved heat transfer observed for 
water near the critical temperature. Whether or not such 
bubble-like activity will occur may well depend on the varia- 
tion of density with temperature of the particular fluid as 
well as on the temperature gradients existing near the surface. 
It may be quite possible, therefore, that such a phenomenon 
might occur near the critical point in some fluids but not in 
others. Further work to clarify the mechanism of heat 
transfer near the critical point should be most interesting and 
timely. 


Flow With Heat Generation 


One of the special problems which has arisen out of the 
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nuclear reactor field is that involving heat transfer in a fluid 
with internal heat generation. This problem occurs princi- 
pally in reactors in which the fuel is in liquid form and where 
this fluid itself is circulated. This problem has been analyzed 
for turbulent flow in a circular pipe (74) assuming a fully 
developed velocity profile and a constant axial temperature 
gradient across the pipe. In the analysis, the turbulent ex- 
change coefficient for heat transfer was set equal to that for 
friction. The latter was derived from the known velocity 
profile. With these assumptions, the differential equation for 
heat transfer could be solved. The solution was given in 
terms of the temperature distribution which prevails when 
the pipe walls are insulated. It was furthermore shown that 
solutions for other wall conditions could be obtained easily by 
superimposing these temperature profiles on the normal pro- 
files which exist in the absence of internal heat generation. 

In order to check the analyses, an ingenious experiment was 
devised (74). An electrolytic solution was used in the test 
fluid so that the internal heat generation could be obtained 
by conducting electric current directly through the fluid. 
Temperature and heat flow rate measurements were made, 
and the results were in satisfactory agreement with the pre- 
dicted values. The analysis, although involving certain 
simplifying assumptions, is certainly accurate enough to be of 
great value in guiding the designer faced with problems in- 
volving internal heat generation. 


Closing Remarks 


As stated earlier, an attempt has been made in this paper to 
briefly review some of the problems in low velocity convective 
heat transfer. Most of the subjects discussed are directly 
related to design applications of current interest, and, in fact, 
the problems usually have arisen from such applications. 
A full understanding of the various heat transfer mechanisms 
would be most beneficial in further advancing the designs in 
question. 

The subjects reviewed are all quite complicated and most 
difficult to approach analytically. As a consequence, the 
designer as well as the research engineer has to rely heavily 
on experimental information. Most of the problems men- 
tioned are in the stage in which such experimental work is 
being pursued actively and in which initial steps toward an 
analytical solution have been taken. Few of the problems, 
however, can be said to have been fully solved or to be com- 
pletely understood. It should be a genuine challenge to the 
research engineer to pursue these problems to a satisfactory 
conclusion. 
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Hypersonic Shock Tunnel 


A hypersonic shock tunnel has been developed for obtaining fluid mechanic information at the 
high Mach numbers and corresponding stagnation temperatures encountered in flight by long range 
ballistic vehicles and satellites. This report describes the hypersonic shock tunnel and presents 
some of the results obtained in the driven tube and in the nozzle. A stoichiometric mixture of 
oxygen and hydrogen with an excess of helium is ignited in the driver to produce strong shock waves 
in air. A shock velocity in air as high as 55,000 fps with a calculated equilibrium temperature of 
16,000 K has been produced in the driven tube. The effects of high stagnation temperatures upon 
the detached shock wave and the pressure distribution for blunt bodies have been observed in the 
nozzle test section. The detachment distance decreased greatly at high temperatures. The pres- 
sure distribution for the hemisphere was found to be less than that predicted by the modified New- 
tonian theory. Shock wave boundary layer interaction at the leading edge of a flat plate was ob- 
served, and the results agreed with the analytical prediction. A detached shock wave was observed 
for a blunt two-dimensional body at very low densities in the test section with a flow Mach number 


of 19.6. 


IHE PRIMARY purpose of the hypersonic shock tunnel_is 

the investigation of high temperature gas effects on the 
flow around bodies traveling at hypersonic Mach numbers 
encountered by intercontinental ballistic missiles and satellites 
re-entering the Earth’s atmosphere. Representative Mach 
numbers and corresponding equilibrium stagnation tem- 
peratures encountered by such vehicles at an altitude of 
100,000 ft are Mach number 20 at 6500 K and Mach number 
26 at 8000 K, respectively. At such temperatures, the air 
may no longer be considered an ideally behaving mixture of 
diatomic oxygen and nitrogen; the gas molecules will dis- 


and the shock tunnel will be described. The combustion 
technique used in the driver for consistent results will be 
discussed in conjunction with the safety precautions for the 
operation. Extremely strong shock waves produced in the 
driven section of the shock tube will be discussed. Some 
preliminary results obtained in the test section and the flow 
over blunt bodies will be presented. The latter include 
surface pressure measurements as well as Schlieren and 
luminosity photographs of different bodies showing high 
temperature effects. Shock wave boundary layer interaction 
at the leading edge of a flat plate has been observed and is 


sociate, ionize, and chemical reactions will occur. These compared with the analytical prediction. di 
microscopic phenomena are usually referred to as “real gas in 
effects.” is 
At the present time very little information is available Description of the Hypersonic Shock Tunnel st 
concerning the effects of these phenomena upon the inviscid Driver and Diaphragm Section m 
and viscous flow fields of hypersonic bodies. The hypersonic The hypersonic shock tunnel consists of a driver, a * 
diaphragm section, a shock tube of constant area and a nozzle 
simultaneously both the flight Mach numbers and stagnation ay inatalied we dump nak, as shown in Figs. 1 and 2. ar 
Seiapenatuwes squociaied with practical: eypireenie vehicles The driver section is 20-ft long with an ID of 6 in. and a wall ta 
tees (5 to 8) thickness of 4 in. Maximum design driver operating pressure th 
in which only high stagnation temperatures can be achieved. oom of til 
Pir enihaiis)-with ole ub the: working ana; the flow Mach safety. To prevent corrosion, the inside surface of the driver in 
number behind a reasonably strong shock wave is only about mi 
Where Mach is un- driver length was determined 80 that the reflected driver W 
bids tenting, ean expansion would not be the limiting factor in nozzle testing sal 
This paper will present a description of the hypersonic 
the equating flome of the possibilit) of detonation in the driver using a stoichiometric 
instrumentation that has been developed for os hea tid mixture of hydrogen and oxygen with helium as a diluent, vail 
18 modified automotive-type spark plugs are mounted flush lis 
gibi, TAS with the inner surface in a spiral arrangement along the tube. maa 
Ballistic Missile Division, ARDC, U. S. Air Force under contract automotive-ty pe coil to each sp ark plug. All the capacitors cla 
AF 04(645)-24. are discharged by a single solenoid switch to simultaneously 
olytechnic Institute, Lroy, N. Y. : 3esides having ignition take place simultaneously at a 
Ce. & Dept. Philsdeiphie Electric number of locations, the gases must be well-mixed before 
‘ Designer, Gas Dynamics d eat combustion. Therefore, the gases are injected at six loca- de 
5 Numbers in parentheses indicate References at end of paper. tions, each inlet having an orifice of four 0.018-in. diam holes the 
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Fig. 1 Idealized picture of shock tunnel phenomena 


drilled at an angle of 45 deg, to produce spiral motion. The 
driver is evacuated to a low pressure before the oxygen is 
injected.. Then about 50 per cent of the total helium charge 
is injected. Next, the hydrogen is added at high pressure in 
stoichiometric proportions with oxygen. Finally, the re- 
mainder of helium is injected at high pressures to produce a 
swirling motion in the mixture and to clear the charging lines 
of combustible gases. With this loading procedure, the oxygen 
and hydrogen mixture is always diluted by a considerable 
amount of helium during the loading process. This is impor- 
tant, since normal detonation pressures as high as 50 times 
the initial pressure (9) or, in certain cases, as high as 150 
times (10) have been observed. After the gases have been 
injected into the driver there is a waiting period of about 5 
min before the mixture is ignited by the multiple sparks. 
With this method of loading and igniting the gases, extremely 
reliable and reproducible driver combustion results have been 
attained. 

The driver section is connected to the driven tube by a 4-ft 
conical transition piece. At the entrance of the transition 
section, there is a section which holds the diaphragm. Three 
hydraulic pistons operating at 3000 psi pressure are used to 
move the clamping unit of the diaphragm section. Check 
valves and mechanical locks prevent the opening of the 
clamping section during the combustion process. 

Stainless steel, copper, monel and aluminum have been 
used as diaphragm material. For most operations, stainless 
steel has been utilized, since it has reasonably high yield 
strength. The stainless steel diaphragms open uniformly and 
do not shatter as they are slapped against the inside wall of 
the transition piece. To permit the diaphragm to open into 
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petals, two grooves are cut at right angles to each other to a 
depth of one-third to one-half the material thickness. A cutter 
radius of 0.020 in. is maintained in order to keep the stress con- 
centration consistent. The pressure at which the diaphragm 


Fig. 2a General Electric Research Laboratory hypersonic 
shock tunnel 
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bursts is controlled by varying the material thickness and the 
depth of the cuts. Consistent bursting pressures have been 
achieved by this method of fabricating the diaphragms. 

To permit longitudinal recoil of the driver section when 
strong shocks are produced, the driver is anchored to the 
floor, as shown in Fig. 2b, by a large channel section which is 
bolted through rubber inserts to a steel beam embedded 
in the floor. The transition piece at the clamping section is 
designed to permit movement of the driver relative to the 
driven tube. Since the diaphragm is stainless steel with a 
maximum thickness of 3 in., the diaphragm will always be 
much weaker than the end flange of the driver and will 
therefore act as a safety valve. 

The use of hydrogen and the possibility of detonation 
with its extremely high peak pressures rejuire precaution to 
protect the operating personnel and the equipment. During 
the charging of the driver, the area housing the driver is 
ventilated by four exhaust fans on the blowoff roof of the 
building. This eliminates the possibility of hydrogen 
accumulating in the driver room. All lines bringing the 
gases to the driver are designed for an operating pressure of 
30,000 psi. A blast mat was installed between the driver 


Fig. 3 Shock tunnel nozzle and dump tank 
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Fig. 2b Combustion driver and diaphgram section 


and the driven tube, as shown in Fig. 2b, to contain an 

fragments within the driver area. The back wall of the 
driver room was constructed of material that will blow out 
under a small pressure differential. 

The gases are stored in standard high pressure bottles 
which stand along the outside wall of the driver room. All 
lines coming from these bottles to the loading panel are 
designed for an operating pressure of 6000 psi. The loading 
panel is located on the opposite side of the blast mat from 
the driver room. The entire hypersonic shock tunnel 
facility is located in a reinforced concrete building with stecl 
doors. The control room, to which all personnel are 
evacuated for high pressure tests, is adjacent to the shock 
tunnel room and is separated from it by a 12-in. thick rein- 
forced concrete wall with two high strength shatterproof 
observation windows. 


Driven Section of Shock Tube 

The constant area or driven section of the shock tube 
consists of stainless steel tubing with a wall thickness of 
1 in. and an ID of 4 in. In order to minimize. the wall 
viscous effects, the inside surface was bored and honed to 
obtain an extremely smooth finish. This section was de- 
signed for a maximum operating pressure of 10,000 psi. 
Adjacent sections of the tube are joined with a tapered fit 
and sealed by two O-rings. This arrangement eliminates 
both inner surface discontinuities at the joints and misalign- 
ments that might otherwise result from the impact loadings 
experienced when the shock tunnel is fired. A total driven 
section length of 103 ft is attained with ten 10-ft sections 
and one 3-ft section which is connected to the conical nozzle. 
These sections have varying numbers of holes in their walls for 
installation of pressure pickups, heat gages and other instru- 
mentation. 

The entire shock tube is supported on a welded rail and 
I-beam structure which is anchored to the steel beam im- 
bedded in the floor as shown in Fig. 3. At the downstream 
end, the tube is anchored to the floor with large steel channel 
members (Fig. 3). A thick rubber slab is used between the 
channel member and the steel beam imbedded in the floor to 
absorb some of the impact produced when the shock wave 
reflects from the nozzle entrance at the end of the tube. 

The end of the shock tube is fitted with a threaded collar 
for clamping the second diaphragm at the entrance to the 
conical nozzle. O-ring seals permit the tube to move relative 
to the nozzle in order to accommodate thermal expansion as 
well as the recoil caused by the reflected shock wave. 

It is important for shock tube operation to know the 
pressure in the driven tube accurately before firing. To 
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measure the initial pressure in the driven tube, McCloud gages 
are used for the pressure range of 1 » to 2 mm of mercury. 
For intermediate pressure to atmospheric, both Bourdon- 
type and mercury manometers are utilized. For low pressure 
operations, provisions are made for evacuating the driven 
section by means of a large vacuum pump. For operations 
where pressures above atmosphere are required, clean dry 
air is led into the tube from the building air supply. 


Nozzle, Test Section, Dump Tank and Sting Support 


\ conical nozzle of the reflected type with a total expansion 
angle of 30 deg and an exit diameter of 24 in. is connected to 
the end of the shock tube and is partly housed inside a dump 
tank as shown in Fig. 3. The area ratio for the conical 
nozzle is varied by changing a short throat section. The 
throat diameter can be varied from 4 to 0.05 in. Throat 
di:meters of 1, 0.33, 0.19 and 0.10 have been tested to produce 
flow Mach numbers in the test section of approximately 10, 
16. 20 and 26 for nearly perfect gas conditions. 

\t the entrance to the nozzle an aluminum diaphragm, 
cross cut in the same manner as the driver diaphragm, is used 
to permit the shock tube and the dump tank to be evacuated 
to different pressures. Usually the dump tank and the 
nozzle are evacuated to low pressures in order to facilitate 
flow establishment in the test section. For very high flow 
Mach numbers, 20 and 26, the dump tank is evacuated to 
about 7 » of mercury. At lower Mach numbers, it is only 
necessary to evacuate to pressures below about 100 yu of 
mereury. The second diaphragm opens upon the arrival of 
the incident shock wave. The bursting pressure is controlled 
again by using different material thicknesses and cut depths. 

The shock tunnel was designed for three possible nozzle 
arrangements: (a) straight-through configuration with the 
throat diameter equal to the shock tube diameter of 4 in., 
(b) converging-diverging or reflected shock wave configura- 
tions, and (c) a straight-through double nozzle arrangement 
(4) which removes the boundary layer. The reflected nozzle 
arrangement was selected for initial operation because the 
shock tube boundary layer effects are minimized and higher 
nozzle stagnation temperatures may be attained for a given 
driver pressure. Moreover, the large area ratios required for 
high flow Mach numbers can be achieved without a pro- 
hibitively large nozzle exit diameter or the complications of 
the double nozzle arrangement. 

A 200-ft? dump tank houses the test section of the conical 
nozzle as shown in Fig. 3. The large dump tank insures 
moderate steady-state pressure in the facility even for 
maximum driver charges. Observation windows are located 
at the exit of the nozzle to permit optical investigations of the 
flow field around models. The windows are optically selected 
plate glass. The dump tank was designed for operation 
with the long side either horizontal or vertical. A vertical 
arrangement for the dump tank was selected initially to 
permit the shortest length for the model sting. Even in this 
position, the unsupported sting length is 5 ft. 

The sting is a 6-in. OD hollow cylinder which tapers to a 
l-in. OD tube behind the model. The sting support struc- 
ture is separate from the dump tank and the sting rides on 
Teflon bushings so that a considerable degree of vibration 
isolation is achieved. The sting is actuated by a lead screw 
which is driven by a reversible d-c motor. 


Instrumentation 


Pressure Gages 


Piezoelectric pressure gages (11) have been used to measure 
pressures in the shock tube and on the models in the nozzle 
test section. The standard Kistler SLM quartz pressure 
pickups have been successfully used in some applications, but 
for many model investigations their large size and small output 
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make them unsatisfactory. Consequently, small barium- 
titanate pressure pickups have been developed to meet 
these needs. This type of pressure gage, while much smaller 
than the Kistler, nevertheless has considerably greater out- 
put. These pickups are more sensitive to vibration than the 
Kistler and must be shock mounted in rubber. The response 
time of the barium-titanate gages is about 10 microsec as 
compared to 50 microsee for the Kistler gages. Barium- 
titanate crystals are more sensitive to heat than quartz 
crystals; hence, the pickup is designed with a certain amount 
of heat insulation depending on the installation for which the 
gage is intended. 

For pressures above 5 psia, the Kistler gages have very 
linear and consistent calibrations. Below 5 psia, both the 
Kistler and the barium-titanate are sensitive to the instal- 
lation in models. For model pressure measurements, a 
dynamic calibration of the gages in their actual installations 
is made. This is accomplished by mounting the model at 
the end of the constant area driven section in place of the 
nozzle and subjecting it to the pressure behind reflected 
shocks obtained with weak incident shock Mach numbers 
(about 1.5). 

For higher pressures, the Kistler and barium-titanate gages 
are fitted with adapters which decrease the effective area on 
which the pressure acts. With these adapters, it is possible 
to calibrate and use the gages above 10,000 psi. 


Ionization Gages 


A simple air gap with 300 v d-c between the center electrode 
and the outer wall is used as an ionization pickup. The 
ionization of the air due to strong shocks breaks down the 
resistance across the gap, and the resultant voltage change 
is used as a time of arrival signal for shock velocity measure- 
ments. For a gap of 0.03 in., the response time of the gage 
is of the order of 1 microsec. 


Heat Gages 


A thin platinum film, sputtered or painted, on a quartz 
or glass piece acts as a resistance thermometer, and the output 
of such a device, when a current is passed through it, is 
related to the surface heat transfer rate. This type of gage 
has been developed at various laboratories in recent years 
(12 to 15). The resistance of the platinum film is usually 
about 40 ohms. A current of 50 milliamp is used to detect a 
voltage variation produced by the change of film resistance. 
The thin sputtered heat gages have a response of about 1 
microsec, but are not as durable as the painted ones. For 
shock waves too weak to ionize the gas in the shock tube, the 
platinum film heat gages have been used to detect the arrival 
of the shock wave at stations along the shock tube. Similar 
heat gages can be used on models in the test section to measure 
the heat transfer rate at various stagnation temperatures. 
At very high temperatures for which air becomes ionized, 
the heat gage must be coated with a very thin electrically 
insulating coating, such as silicon dioxide as described in 
(15). This coating prevents the ionized gas from short 
circuiting the platinum film at a small loss in response time. 


Determination of Shock Velocity 


The instrumentation used to determine the shock wave 
position as a function of time in the 100-ft long driven tube 
consists of heat or ionization gages which serve as time of 
arrival pickups, a modified Tektronix 535 oscilloscope, a 
crystal controlled raster generator, and a display chassis. 
The raster generator produces up to 500-cm sweep length 
with a nearly vertical sawtooth or raster arrangement in which 
each vertical line represents a given time. The output of 
12 sensing elements in the shock tube wall is fed into the 
display chassis which produces signals that are displayed as 
horizontal markers on the vertical sawtooth pattern. In 
normal operation of the shock tunnel, the shock velocity 
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measurements are supplemented by a Berkeley counter 
which measures the shock wave traversal time over the last 
6} ft of the tube with 1 microsec accuracy. 


Optical Methods 


A single pass Schlieren system with 6-ft focal length para- 
bolic mirrors is used to obtain photographic records of the 
model flow fields. A spark of about 0.2 microsec duration, 
produced by the discharge of a capacitor, is used as the light 
source. For high temperature conditions with a large 
amount of luminosity from the hot gas, an optical shield is 
placed at the knife edge to permit only small amounts of 
luminosity to reach the film. It is thus possible to photo- 
graph the shock wave with the luminous gas region. 

Direct luminosity or “flash’’ photographs have been taken 
by leaving the camera open at the test section window during 
the entire test. In this manner, both color and black and 
white photographs of the luminous region around blunt 
bodies have been obtained. For investigations of the flow 
phenomena as a function of time, direct photographs have 
been taken with a Fastax camera operating at 8000 frames 
per sec. 

For spectroscopic investigations of the heated gas after a 
shock wave, quartz windows may be inserted in the instru- 
mentation holes in the shock tube walls. Some preliminary 
results have been obtained in the shock tube and in the test 
section. 


Preliminary Results From the Shock Tunnel 


Combustion Driver 


Initial combustion investigations (9) were made in a 10-ft 
section of the driver in order to determine combustion per- 
formance using gas mixtures of helium and stoichiometric 
hydrogen-oxygen. The driver “loading’’ pressure was varied 
from 25 to 500 psia, and the mixture ratio from 45 to 80 per 
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cent helium. Tests were performed with one to nine spark 
plugs igniting the combustible mixtures. Two Kistler 
pressure pickups were installed to obtain the pressure time 
history during the constant volume combustion tests. 

The observed peak pressures for smooth burning were 
close to the theoretical values calculated by assuming no 
dissociation (Fig. 4). Richer mixtures and higher loading 
pressures tended to produce detonation in the driver. At 
25 psia loading pressure, detonation occurred in mixtures 
with 45 per cent or less helium, whereas at 300 psia loading, 
detonation occurred with a 65 per cent helium mixture. 
With detonation, peak pressure to loading pressure ratios 1s 
high as 50 were measured at the higher driver loadings; 
ratios as high as 150 have been observed in a long tube (1(). 
It was determined during these studies that a 70 to 75 pr 
cent helium mixture with a balance of stoichiometric hydr- 
gen and oxygen was lean enough to prevent detonation with 'n 
the range of driver loading pressures investigated. 

It was observed that the time to reach peak pressu’e 
increases with the amount of helium in the mixture and wit! 
increased loading pressures as indicated in Fig. 5. With nine 
sparks, the time to peak for 70 per cent helium is 10 to 15 
milliseec. As the detonation mixture ratio at a given loadin: 
pressure is approached, the time to reach peak pressur: 
drops rapidly to about 5 millisec, and the amplitude of th: 
pressure oscillations increases greatly. 

These combustion results indicate that great flexibilit 
can be obtained in the driver characteristics by varying the 
percentage of helium and the number of ignition sparks. 

Following the assembly of the complete 20-ft driver section. 
another series of combustion investigations were conducted 
The results of these tests were quite similar to those in the 
10-ft section except for the longer period of pressure oscilla- 
tions when the mixture approached the detonation limits. 
Thus, with the technique currently in use in the Research 
Laboratory shock tunnel, it has been possible to obtain 
consistent repeatable combustion driver performance. 
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Fig. 5 Time to peak vs. initial pressure 
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Results in the Shock Tube 


At low shock Mach numbers, the agreement between the 
observed shock velocity at the 40-ft station and the predic- 
tions of the idealized shock tube theory (16, 17) was reason- 
able, allowing for shock wave attenuation (Fig. 6). However, 
as the shock Mach number is increased, the trend of the 
experimental results for the pressure ratio across the driver 
diaphragm vs. the shock Mach number varies from that of 
the theory (Fig. 7). For shock Mach numbers above 25, 
the observed velocities were much greater than one would 
predict from the simplified theory. Similar results with a 
cold driver were observed and investigated in more detail 
(18). The results of the investigation indicated that better 
agreement between theory and experiment could be obtained 
by taking into account the finite opening time of the dia- 
phragms and the mixing of the driver and driven gases at 
the contact surface. In addition, with the present case 
where combustion driving is utilized, it is possible to produce 
shocks greater than predicted by theory by letting the dia- 
phragm burst well before the combustion is completed, so 
that an approximately constant pressure combustion situation 
exists during the remainder of the burning of the gases (19). 
This last effect was appreciable in producing the hyper- 
velocity shock wave. A typical variation of the shock Mach 
number with distance is shown in Fig. 6. 

To date, a maximum shock velocity of 55,000 fps has been 
produced with combustion driving and air as the driven 
gas. The calculated gas velocity after the shock wave was 
approximately 52,000 fps, and the equilibrium temperature 
was 16,000 K. Thus, it may be possible to investigate the 
stagnation heat transfer for a blunt body at a flight Mach 
number of about 50. 


Determination of Flow Conditions in the Nozzle 


All tests so far have been conducted with the reflected-type 
nozzle, and thus the conditions behind the reflected shock 
wave were the nozzle stagnation conditions. Both the inci- 
dent shock velocity just before the end of the tube and the 
pressure behind the reflected shock wave are measured for 
each test. The gas behind the reflected shock is assumed to 
be in thermodynamic equilibrium so the temperature can 
readily be computed by numerically solving the normal 
shock equations for both the incident and reflected shocks, 
based upon the thermodynamic properties for air (20 to 22). 

Test section conditions can be estimated from one-di- 
mensional isentropic nozzle calculations and measurements 
of reflected pressure and test section static or impact pressure. 
For low temperatures where air is nearly an ideal gas, the 
standard y = 1.4 results are applicable. For higher tem- 
peratures, a series of calculations based upon equilibrium 
assumptions and carried out for a range of stagnation condi- 
tions are utilized. 


Preliminary Nozzle Flow Results 


The hypersonic shock tunnel may be operated at different 
nozzle flow stagnation temperatures by varying the strength 
of the incident shock at the end of the constant area section. 
The strength of the incident shock is a function of the shock 
tube initial conditions. Initial shock tunnel investigations 
were conducted in the nominal Mach 10 nozzle having a 
throat diameter of 1 in. and an exit diameter of 24 in. The 
reflected stagnation temperature was high enough, about 
1100 K, to avoid condensation in the nozzle, yet low enough 
so that the air may be considered a perfect gas with constant 
ratio of specific heats of 1.4. Thus, the ideal gas calculations 
for air were applicable. These nearly ideal gas results served 
as the reference for the investigations at higher temperatures 
for a given nozzle area ratio. 

Reflected equilibrium temperatures, at the entrance to 
the nozzle, as high as 12,000 K have been produced. For 
stagnation temperatures greater than about 2000 K, the real 
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Fig. 7 Diaphragm pressure ratio vs. shock Mach number at 
40.5-ft station 


gas effects must be considered during the expansion process 
in the nozzle. The flow Mach number in the test section 
decreases as the stagnation temperature is increased, provided 
the flow during expansion in the nozzle is in thermodynamic 
equilibrium. Results of a one-dimensional isentropic nozzle 
calculation for an equilibrium nozzle expansion show a 
decrease in flow Mach number with stagnation temperature. 
Some preliminary investigations into the question of whether 
an “equilibrium” or “frozen” nozzle expansion takes place 
seem to indicate that if the nozzle stagnation pressure is 
high enough the gas will be nearly in thermodynamic equi- 
librium. At low pressures, the indications are that the 
expansion is nearly a frozen one. As used in this paper, a 
completely “frozen” expansion is one in which only the 
translational and rotational energies of the gas change from 
their stagnation value, while all the others remain at their 
stagnation limits as the gas expands in the nozzle. Further 
research is being conducted to obtain more information 
regarding the state of the air during the expansion process. 
Flow establishment times are determined by measurement 
of the impact pressure at the exit of the nozzle as a function 
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Fig. 8 Impact pressure vs. time for Mach 10 nozzle with 7; = 
3000 K 


Fig. 9 Five-in. span 40-deg wedge in Mach 10 nozzle with 
Ts = 1300 K 


Fig. 10 Hemisphere in Mach 10 nozzle with 7; = 1050 K 
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of time. A representative impact pressure time history for 
the nominal Mach 10 nozzle is shown in Fig. 8. For this 
particular test at 3000 K reflected stagnation temperature, 
the starting time of the flow was about 0.7 millisec, and the 
flow duration was about 2.5 millisec. In general, the flow 
starting times were found to vary between 0.6 and 1 millisec. 

To check the flow uniformity in the test section, an impact 
rake with five pressure gages was tested at the usual model 
location in the nominal Mach 10 nozzle. The results ob- 
tained at 1100 K reflected stagnation temperature have 
indicated that, within the accuracy of the measurements, the 
transverse impact pressure distribution is uniform at the 
model location. A 40-deg wedge with a 5-in. span was also 
used to detect any large disturbances in the test section, 
such as might be caused by separation of the flow from the 
nozzle wall. The Schlieren photograph taken at a reflecte: 
stagnation temperature of 1300 K in the Mach 10 nozzle is 
shown in Fig. 9. For the 40-deg wedge, the shock wave 
was approximately straight up to the point where the ex- 
pansion from the shoulder starts interacting with the shock 
wave. 

Flow over a hemisphere was investigated at different 
stagnation temperatures in the nominal Mach 10 nozzle. At 
an equilibrium reflected temperature of 1050 K, the impact 
Mach number was 10.J, and the ratio of the bow shock 
detachment distance to the hemisphere radius 6/R was 0.13, 
as indicated in Fig. 10. At the stagnation temperatures of 
approximately 10,000 K, the detached shock wave moved 
closer to the hemisphere and the 6/R ratio decreased to 
approximately 0.09. The calculated normal shock density 
ratio for this test condition, assuming the flow to be in ther- 
modynamic equilibrium, is approximately 15. These results 
indicate that the detachment distance decreases with the 
increase in the shock density ratio caused by the real gas 
effects at high temperatures, as discussed in (26). 

A luminosity photograph of a 5-in. diam hemisphere at a 
reflected stagnation temperature of 5800 K in the nominal 
Mach 10 nozzle is shown in Fig. 11. At this temperature the 
region of luminosity was very sharply defined at the front, 
indicating that the relaxation time after the shock wave was 
very short. The luminosity disappears as the flow expands 
around the hemisphere and the temperature decreases. 
Similar results for the luminosity for a 50-deg cone-hemisphere 
were obtained for an equilibrium reflected stagnation tem- 
perature of 5800 K as shown in Fig. 12. 


Fig. 11 Luminosity photograph of hemisphere in Mach 10 
nozzle with 7; = 5800 K 
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Fig. 12 Luminosity photograph of 50-deg cone-hemisphere in 
Mach 10 nozzle with 7’; = 5800 K 


Preliminary pressure distributions for a hemisphere and a 
50-deg cone-hemisphere have been obtained in the Mach 10 
nozzle. The pressure distribution for the hemisphere at a 
stagnation temperature of 1300 K is shown in Fig. 13. Over 
most of the hemisphere the pressure was lower than that 
predicted by the modified Newtonian theory. Recent 
analytical results (23 to 25) have indicated lower pressures 
than the modified Newtonian at hypersonic Mach numbers. 
At higher temperatures, the real gas effects cause the 
shock wave to move closer to the body without affecting the 
pressure distribution to any extent. A more detailed discus- 
sion of the high temperature effects upon the pressure distri- 
bution and the shock wave shape is presented in (26). 

Flow over a 5-in. wide flat plate was investigated in the 
Mach 10 nozzle at a stagnation temperature of 1560 K. 
For this condition the flow Mach number was approximately 
10, and the Reynolds number per inch was 60,000. A 
Schlieren picture of the flat plate showing the shock wave 
and the boundary layer is presented in Fig. 14. At the 
leading edge, the shock wave and the boundary layer are 
merged together, and they both vary with distance along the 
plate as X°74, In (27), the analytical prediction was X°-7, 
so that the agreement with the experiment is quite good. 
As the distance from the leading edge increases, the shock 
wave and the viscous boundary layer separate. This has 
been observed by others (28 to 30). 

It has been possible to establish flow in the test section 
for the nominal Mach 16 and 20 nozzles. In the latter case, 
with a throat diameter of 0.19 in., the impact flow Mach 
number in the test section was 19.6 for a stagnation tempera- 
ture of 2140 K. The corresponding calculated static pressure 
in the test section for adiabatic flow conditions was about 
15 of mercury. 

It is seen in Fig. 15 that this relatively rarefied Mach 
19.6 flow produced a sharp detached shock ahead of a two- 


6 Private communications from Dr. Eva Winkler of the Naval 
Ordnance Laboratory at White Oak, Maryland, have indicated 
similar results in the hypersonic wind tunnel for the hemisphere. 
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Fig. 13 Pressure distribution for hemisphere in Mach 10 nozzle 
with 7; = 1300 K 


Fig. 14 Flat plate in Mach 10 nozzle with 7; = 1560 K 
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dimensional blunt body. Rarefied flow phenomena have 
been investigated (31) at lower Mach numbers in ideal gas 
flows where rarefied conditions exist before and after a weak 
shock wave. At higher Mach numbers and stagnation 
temperatures, the change in mean free path across the shock 
should be quite large due to the larger shock density ratio 
caused by the real gas effects. Thus, it seems possible to 
have a rarefied flow ahead of a detached shock and continuum 
flow behind. 


Conclusions 


The combustion driver technique, using a mixture of 
stoichiometric hydrogen and oxygen with an excess of helium, 
is a reliable and useful method for producing strong shock 
waves in air. 

At high pressure ratios across the diaphragm separating 
the driver and the driven tube, it is possible to produce shock 
wave velocities in the shock tube in excess of those predicted 
by the simple shock tube theory. 

A shock velocity of 55,000 fps has been produced in the 
shock tube. The calculated equilibrium temperature be- 
hind the shock wave was 16,000 K with the air after the shock 
moving at a velocity of 52,000 fps and with a Mach number of 
approximately 7. 

Hypersonic flows in the conical nozzle have been established 
with durations of 1 to 4 millisec of fairly uniform flow in the 
test section. A flow Mach number of about 20 has been 
obtained in the test section, and with further investigations 
it will be possible to produce even higher Mach numbers. 

Pressure distribution, shock wave shape, luminosity 
regions and heat transfer rates can be obtained for models 
in the test section at hypersonic flow Mach numbers and the 
temperatures corresponding to those encountered in flight. 

Flows over blunt bodies at high stagnation temperatures 
have been investigated in the shock tunnel. The shock wave 
detachment distance decreases significantly with increasing 
stagnation temperature. 

At very high flow Mach numbers and low densities, rarefied 
flow conditions, it has been possible to observe a detached 
shock wave for a blunt two-dimensional body using a con- 
ventional Schlieren system. 

The hypersonic shock tunnel is a useful research tool 
for investigating the flow phenomena associated with the 
high flow Mach numbers and high stagnation temperatures in 
air under continuum and rarefied flow conditions. 
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Thrust Optimization of a Nuclear 
Rocket of Variable Specific Impulse 
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A simple relation is derived between the specific impulse and the flow rate of a propellant gas 
under nuclear reactor operating conditions which allow molecular dissociation of the gas. Utilizing 
this relationship, the optimum variable thrust program maximizing the kinetic energy of a ver- 
tically ascending vehicle is obtained under the assumptions of constant gravity and zero aerody- 
namic drag. The advantage of this program is demonstrated by comparison with the usual 


constant thrust program. 


» INCE a nuclear rocket is a source of practically unlimited 
J energy, the problem of a nuclear rocket is to maximize 
the energy absorbed by the thrust producing working medium 
during its passage through the reactor. The amount of 
energy which a working medium could absorb depends, of 
course, on the temperature in the reactor core, which in turn 
is a property of the reactor materials used. Once, however, 
the best high temperature materials are developed and in- 
corporated into the reactor, there remains the problem of in- 
suring that the energy communicated to the rocket at burnout 
of the total energy injected into the rocket system be a maxi- 
mum. The present paper is mainly concerned with the latter 
problem. 

At temperatures attainable by some possible reactor core 
materials, certain propellant gases will dissociate to a sig- 
nificant degree. With dissociation, the gas will have dras- 
tically increased specific heats. The enhanced specific heats 
will then yield specific impulse values greater than expected 
by the usual square-root-of-temperature law. Thus, molecular 
dissociation is an important property of a gas as a nuclear 
rocket propellant (1).4 

The amount of dissociation products present depends on the 
temperature and pressure of the gases. Although the tem- 
perature is limited by the materials available, the pressure is 
strictly a design parameter, subject to the designer’s selection. 
However, for a fixed size of engine and a given temperature, 
the operating pressure would determine the thrust level. 
Therefore, the thrust level determines the degree of dissocia- 
tion and hence the value of specific impulse which a particular 
gas is able to produce. In general, since high thrust requires 
high operating pressure which tends to suppress gas dissocia- 
tion, a high thrust operation will result in low specific impulse, 
whereas a low thrust operation will result in high specific im- 
pulse. These facts lead to the problem of determining the 
optimum thrust level for a nuclear rocket. 

It is wel] understood that when the specific impulse is a 
constant, efficient use of a rocket powered vehicle requires full 
thrust operation of the rocket, or, if aerodynamic drag is a 
factor, the use of an optimum thrust program which would 
accelerate the vehicle to the maximum before the drag be- 
comes limiting. This is because low acceleration will waste 
energy by unnecessarily lifting the propellants to high alti- 

Presented at the ARS 13th Annual Meeting, New York, N.Y., 
Nov. 17-21, 1958. 

! Associate Manager, Aerophysics Department. Member ARS. 

2 Member, Technical Staff. 

3 Associate Director, Astrovehicles Laboratory. Member ARS. 

* Numbers in parentheses indicate References at end of paper 
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tudes. Thus, in the ideal case of limitless engine capacity and 
zero drag, the most efficient operation would mean an impul- 
sive thrust. However, for engines of practical size and when 
the specific impulse is not a constant, and specifically if it in- 
creases with the decrease of thrust, there will exist an optimum 
thrust program which represents the most efficient operation 
under the conflicting requirements of high thrust to avoid 
working against gravity and low thrust to realize the ad- 
vantage of increased specific impulse. 

In this paper, a mathematical representation is developed 
for the specific impulse of nuclear rockets in terms of the pro- 
pellant flow rate (and consequently engine thrust) on the 
basis of dissociation considerations. To a first approximation, 
it is found that the specific impulse can be represented by a 
simple power of the propellant flow rate multiplied by a con- 
stant. The exponent and the constant multiplier depend 
upon the particular engine and its operating conditions. Using 
the power law, an analysis is carried out to determine the 
optimum thrust program which will maximize the kinetic 
energy of a vertically ascending nuclear rocket vehicle. To 
simplify the analysis, air drag is neglected, and constant 
gravity is assumed. This approach leads to the very simple 
conclusion that for optimum operation the propellant flow 
rate must be reduced with time in such a way that the ac- 
celeration is maintained constant. The performance of a 
vehicle utilizing the constant acceleration program is then 
compared with that expected from a constant thrust rocket. 
The advantage of the optimum program is clearly demon- 
strated. 

It should be noted that problems of thrust programming 
have been studied by other investigators, see for example (2 
to 4). These investigations, however, are dealing exclusively 
with rockets of constant specific impulse. The problem in this 
paper is new, and the thrust program derived is quite different 
from those derived for constant specific impulse rockets. 


Specific Impulse 


Basically the specific impulse (/,,) is obtained by combining 
the propellant discharge momentum and the force resulting 
from the pressure at the nozzle exit in excess of the ambient 
pressure. It can be shown that the following formula yields 
I,» for given gas conditions at the reactor discharge, an isen- 
tropic expansion ratio e and the ambient pressure p.. 
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The first term represents the momentum relative to the ve- 
hicle per unit mass flow rate. The second term represents the 
force per unit mass flow rate due to the static pressure at the 
exit plane. 

To evaluate Equation [1], it is necessary to determine the 
shifting equilibrium composition of the gas during its ex- 


pansion through the thrust chamber. The computation re- 
quired is, in general, very complex, and a high speed machine 
computing program would be desirable. Such a program has 
been employed at the Space Technology Laboratories, Inc., to 
calculate the J,, for Hz as a nuclear rocket propellant. For 
example, in Fig. 1 are shown the results of shifting equilibrium 
calculations as compared with frozen equilibrium values of 
the specific impulse for a nozzle expansion ratio of 40 and zero 
ambient pressure. 

In order to facilitate analytical investigations, some simpli- 
fications will now be considered in the calculation of J,,. 
Since rockets of high performance operate most of the time at 
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high altitudes, the ambient pressure p, may be taken as zero. 
Furthermore, these rockets usually employ high expansion 
nozzles so that p./p. is very small, and the whole pressure 
term in Equation [1] can be neglected. This leaves the first 
term which will enable one to compute J,, from the enthalpy 
differential (H. — H.), namely 


Ip 2aWH. — He 2] 


We shall now investigate this term for a diatomic gas, for ex- 
ample, He. 

The enthalpy per unit mass of the heated working flui:| 
may be expressed as 


The first term represents the contribution of the atoms to th: 
total enthalpy, and consists of the energy of their translationa 
motion and the energy of molecular dissociation. The secon 
term represents the molecular contribution to the enthalp, 
due to the translational, rotational and vibrational energies o 
the molecules. At the temperature of the reactor discharge 
T = T., the gas molecules may dissociate to a significant de- 
gree, resulting in the first term being the most important. At 
the nozzle exit, 7 = T., the recombination of the atoms is 
usually complete, and the enthalpy expression will consist of 
the second term only. 

Reference to Fig. 2 illustrating the enthalpy value of hydro- 
gen, shows that the monatomic contribution can indeed be 
dominant at a high temperature such as may exist at the reac- 
tor discharge. Assuming this to be the case and neglecting 
the molecular contribution, the change in enthalpy per unit 
mass becomes 


fi 


where h(7.) is a function of T, only. 

The weight fraction f; of the atoms is dependent upon pres- 
sure as well as temperature. By definition, the pressure 
equilibrium constant K, is given by 


where p, is in atm, and K, isa function of T. only. Utilizing 
the relation f; + fe = 1, and assuming that K, < p. (which is 
generally true in the temperature range of interest), one may 
readily derive the following relation for the weight fraction f, 


fi [6] 


Considering Equations [4 and 6] together, we see that the 
enthalpy differential (H. — H.) for a given temperature 7’ is 
inversely proportional to the square root of the reactor dis- 
charge pressure. 

We further note that the weight flow rate G through a unit 
area of the nozzle throat is given by 


G = pg/C* [7] 


Calculations show the characteristic velocity C* to change by 
only a few per cent at temperatures of interest as the pressure 
p- varies from 2 to 20 atm. Thus we may consider G to be 
proportional to p. and, therefore, conclude that for a given 
temperature 7’, the enthalpy differential (H. — H.) is in- 
versely proportional to the square root of G. Since J,» is 
related to (H. — H.) by Equation [2], if we use the total flow 
W instead of G, we obtain 


sp = a( —W) 


RT. = fih(T-) [4] 


a, B, -w>0 [8] 


where a@ is a constant, B = 1/4. 
It should be pointed out that Equation [8] is valid neither 
for an undissociated gas nor for a completely dissociated gas. 
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Thrust Program 


The problem here is to determine an optimum thrust pro- 
gram for a nuclear rocket of variable J,,.. The J,, is assumed 
to vary according to Equation [8]. Although there are many 
trajectories of interest, we shall confine the following discus- 
sion to the trajectory offering simplicity of result and present- 
ing a clear demonstration of the effect of dissociation on nu- 
clear rocket propulsion. Thus the vehicle is assumed to as- 
cend vertically under constant gravity and with zero drag. 
The thrust program is to be such that for a given take-off 
\ eight and a given amount of propellant, the velocity gained 
shall be a maximum. The first step of the analysis will be 
to establish the optimization equation. The second step is 
t» derive the optimum program for the type of reactor under 
consideration by applying the optimization equation. The 
third step is to compare the performance of such a programmed 
rocket with that of the more conventional constant thrust 
} rogram. 


Optimization Equation 


Under the above conditions, the equation of motion of the 
rocket is 


= w ( 
dt at 9] 
since 
dp do dw do 
dt dwdt dw [10] 


quation [9] reduces to 
y= — — } dw 
ca g ( + ;) u {11] 


which by integration gives the velocity 


wb 1 


In order for the velocity to have a stationary value, it is 
necessary that its variation év be zero; thus we obtain the 
optimization equation 


=w {13] 


Optimum Program 


We shall now determine the optimum thrust program for a 
nuclear rocket governed by Equation [8]. To do this, one 
can solve Equation [13 and 8] simultaneously as follows: 

Differentiating Equation [8] yields 


Ol sp 
ap ) [14] 
Substituting Equation [14] into Equation [13] results in the 
following equation 

aB(-w)'-8§ =w 


or 


(= 
a8 (15) 


Integrating Equation [15], the “burnout” time is obtained 
ii-s(1-8 —8/1—8 
t, = (aB)"’ wy — W (16] 


May 1959 


Since the vehicle acceleration is given by 


dt w 


substituting Equations [8 and 15] into [17] and then integrat- 
ing will give the velocity 


i- *) 
ot) =g9(—} 18 
g ( B [18] 
The altitude at time ¢ can be obtained by integrating Equation 
[18], thus 
1 1-8 
= 59 ( [19] 


To obtain the ‘‘burnout”’ velocity and altitude as a function of 
the launch weight and “burnout” weight, one can simply sub- 
stitute Equation [16] into Equations [18 and 19]. 

To determine how the specific impulse varies with time, it is 
most convenient to examine Equation [13] which may be 
written 


wt Ol sp wt dI,, dt 
2 
ow dt dw 
or 
al, (d*w/dt?) 
dt (dw/dt)? 
Upon noting the identity 
d (d*w/dt?) 
dt \w (dw/dt)? 
Equation [21] may be expressed as 
d 
= 23 
dt dt (“) 
Hence 
I.p = t — (w/W) + constant 24) 
Referring to Equation [20] again, [24] becomes 
ly +W + constant [25 
8 consta 2 
P ow 


This is a general result whatever the dependence of specific 
impulse on mass flow rate. From relation [8] 


Ol sp Bi. 
age) 
[26] 
so that 
I.p = t/(1 — B) + constant [27] 


the constant being determined from satisfaction of Equation 
[13] at time of launch. 


Rocket Performance 


Equations [16, 18 and 19] give the performance of the nu- 
clear rocket when operated under the thrust program which 
gives the maximum burnout velocity. Utilizing the following 
parameters, the rocket’s performance was calculated and the 
results are presented in Figs. 3, 4 and 5. 


a = 4500 [28] 
B =+1/4 [29] 
wy = 200,000 Ib [30] 
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Fig. 5 Weight, thrust and acceleration vs. time 


In Fig. 3, the effective launch velocity ver is a measure of the 
total energy per unit mass being defined as 


vert = Vv? + [31) 


For comparison, the performance of two constant thrust ve- 
hicles is also presented. For one, the initial thrust to weight 
ratio and initial specific impulse were assumed to be the samc: 
as for the programmed rocket, while for the other, the sam: 
initial specific impulse was assumed but with a reduced initia! 
thrust to weight ratio of 1.5. It is apparent that the pro- 
grammed thrust rocket greatly outperforms the constan' 
thrust vehicles even at moderate mass ratios. Moreover, a‘ 
the higher velocities the constant thrust vehicles are under- 
going an acceleration an order of magnitude greater than th: 
constant acceleration rocket. 

In Fig. 4 the variation of specific impulse and propellant 
flow rate with time is shown for the programmed vehicle 
while Fig. 5 presents the variation of the thrust and remaining 
weight with time. The last two quantities vary in such « 
way that the vehicle acceleration remains constant. 


Conclusion 


The present analysis has shown that, for a nuclear rocket 
of variable specific impulse, there is definitely an optimum 
way of programming the thrust to fulfill a particular mission. 
This is because the dissociation characteristics of gases are 
such that, for a given temperature, a decrease in the operating 
pressure (and thereby the thrust) is accompanied by an in- 
crease in the specific impulse. In the case considered, the 
optimum thrust program results in a much better performance 
than the usual constant thrust operation. 


Nomenclature 
a = constant V2/g 
A* = throat area, ft? 
C* = characteristic velocity, fps 
C, = specific heat at constant pressure, ft lbf/R lbm-mole 
Dy) = dissociation energy, ft lbf/lbm-mole of diatomic gas 
f = weight fraction 
g = conversion factor, lb/slug; gravitational acceleration, 
fps? 
G = weight flow rate per unit throat area, lbf/ft? sec 
= enthalpy, ft Ibf/lbm 
I;p = specific impulse, lbf/(lbm/sec) 
K, = pressure equilibrium constant 
m = mass of vehicle, slug 
M molecular (atomic) weight 
p = pressure, lbf/ft? 
R= universal gas constant, ft lbf/R lbm-mole 
t = time, sec 
v = velocity 
w = mg, weight of vehicle, lbf 
y = altitude, ft 
a, 8 = constants 
€ = area ratio of nozzle, A;/A* 
Subscripts 
a = ambient condition 
b = burnout 
c = reactor discharge or stagnation condition 
e = exit condition 
0 = take-off 
, 2 = monatomic and diatomic, respectively 
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Similarity Parameters for Surface 


Melting of a Blunt Nosed Body in a 
High Velocity Gas Stream vers: 


Space Technology Laboratories, Inc. 
Los Angeles, Calif. 


By employing a simple model for the two-phase boundary layer produced by surface melting, 
“scaling laws”’ are obtained for the regime in which vaporization does not occur at the gas-liquid 


interface. 


When the liquid is fairly viscous, the gas boundary layer is virtually independent of the 


liquid film. The temperature rise across this film depends not only on the local heat transfer rate, 


but also on the difference in total enthalpy across the gas boundary layer. 


Thus, tests in rocket 


exhaust jets may not simulate vaporization or film instability phenomena that could occur in flight. 
Calculations of the behavior of Corning glass No. 7740 are carried out for some typical flight condi- 


tions. 


One interesting result is the low value of the liquid film Reynolds number, which indicates 


that the liquid film is almost certainly laminar. This liquid-solid system absorbs about 1125 Btu/Ib, 


or about six times the heat capacity of a nonmelting copper heat sink. 


A study of the effect of 


mass addition to the gas phase shows that the amount of heat ‘“‘blocked”’ by this process is sensitive 
to the total enthalpy difference across the gas boundary layer. The possible increase in effective heat 
capacity with vaporization at the gas-liquid interface is quite substantial at hypersonic flight ve- 


locity. 


HE PRESENT paper contains a simplified analysis for 

the two-phase boundary layer formed when a solid surface 
melts and the molten material flows around the body under 
the influence of the aerodynamic forces exerted by the gas. 
When the heat transfer rate to the liquid film is below a cer- 
tain critical value, no appreciable vaporization takes place 
at the gas-liquid interface, and the flow in the gas phase is 
virtually unaffected by the presence of the extremely thin 
liquid film. The present study is limited to this regime, but 
one of the objectives is to bring out the parameters that 
govern the temperature rise across the liquid film, and thereby 
determine the critical conditions for vaporization. The effects 
of vaporization are discussed later.* 

Sutton (1)4 and Roberts (2) have analyzed the differential 
equations governing the motion of the gas boundary layer and 
the molten liquid film at the forward stagnation point of a 
blunt body, in the absence of vaporization. The results of 
their calculations bring out the importance of the liquid vis- 
cosity in determining the interfacial temperature and the 
shielding effect of the liquid film. Roberts (2) in particular 
obtained simple approximate expressions for the melting rate 
and other quantities of interest in terms of a parameter con- 
taining the ratio of the driving enthalpy difference for heat 
transfer to the sum of the heat of fusion and the heat capacity 
of the solid up to the melting point. More recently, Bethe 


Received Dee. 1, 1958. 

1 This Po is a revised version of a report with the same title 
dated July 29, 1957, issued originally by Space Technology Lab- 
oratories for limited distribution. 

? Consultant; also Professor of Aeronautics, Guggenheim 
Aeronautical Laboratory. Member ARS. 

’In certain cases a portion of the heated solid material near 
the surface may be pyrolyzed into a gas. Sutton (1) suggested 
that the aerodynamic effect of this process can be roughly simu- 
lated by assuming that a fixed fraction of the melted material 
enters the gaseous boundary layer as a gas. Then, this case also 
can be included in the present study without much additional 
work. (See section titled Effect of Mass Addition to Gas Bound- 
ary Layer.) 

* Numbers in parentheses indicate References at end of paper. 
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This phenomenon can also be analyzed along the lines of the present paper. 


and Adams (3) gave an approximate treatment of this prob- 
lem for highly viscous, glassy materials in which the viscosity 
is strongly temperature dependent. Their approximation de- 
pends upon the fact that the liquid layer is thin in comparison 
with the “thermal thickness” in the material. Therefore, the 
temperature gradient rapidly approaches the exponential 
variation associated with a material countercurrent of uniform 
velocity toward the interface. By utilizing this exponential 
variation as a first approximation in the continuity and 
momentum equations, Bethe and Adams obtain a relation be- 
tween the melting rate and the interfacial temperature. Sut- 
ton’s (4) comparison of their approximate results for pyrex 
glass with numerical integrations of the same differential 
equations shows that the melting rate is underestimated by 
about 30 per cent for a given heat transfer rate at the inter- 
face. The equilibrium vapor pressure corresponding to the 
interfacial temperature is too low by about the same amount. 
Similar conclusions can be drawn by inserting numerical 
values into the ‘second order” approximations derived in the 
Appendix to reference (3). In other words, the variation in 
countercurrent velocity across the liquid layer plays a sig- 
nificant role. Nevertheless, the Bethe-Adams analysis is very 
useful in estimating the effective heat capacity of the ma- 
terial; it is also worked out for the case where vaporization 
oceurs (3), and it can be extended around the body. 

In an unpublished report Carrier (5) developed an integral 
method for this problem, which also starts with an exponential 
temperature variation, but regards the logarithmic decrement 
as an unknown quantity to be determined by satisfying an in- 
tegral of the heat conduction equation multiplied by a ju- 
diciously chosen weighing function. No numerical results 
utilizing this promising method are available at present. 

The conservation laws should furnish a relation between 
the interfacial temperature and the liquid viscosity that is 
relatively insensitive to the details of the velocity and tem- 
perature distributions across the liquid film. We require a 
simple, approximate relation that is applicable not only at the 
stagnation point, but also around the body, including regions 
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of turbulent flow in the gas phase. The liquid film itself is 
laminar, provided the Reynolds number based on film thick- 
ness and liquid properties evaluated at the interface does not 
exceed a certain critical value of about 150. Above this value 
the film is almost certainly unstable® when the gas phase is 
turbulent. Thus it is essential to be able to estimate this 
Reynolds number in terms of the parameters of the problem 
(see next section). Although detailed calculations are proba- 
bly required for precise numerical values, the approximate re- 
lations for interfacial temperature and melting rate obtained 
here may serve as a useful guide. 


Similarity Parameters for the Liquid Film 


In order to fix our ideas, consider the flow over an unyawed 
blunt nosed body of revolution (or a cylindrical body). In 
this case, the properties of the single-phase boundary layer at 
ambient density levels corresponding to flight altitudes below 
about 200,000 ft are now fairly well established. In a coordi- 
nate system fixed in the gas-liquid interface, the problem is 
quasi-steady (Fig. 1). Suppose we start out with the assump- 
tion that u;/u, is so small that the gaseous phase behaves as if 
u; were zero. (The correctness of this statement will be 
tested later.) Then the shear stress at the gas-liquid interface 
is equal to the value it would have at a solid boundary at the 
enthalpy level h,; = h:(T;) 


(T9)i = (Pette?/2)ey = (TL)i {1] 
where c, is known for u; = 0. Also 
(qo)i = PetteCu( hs, — h,,) [2] 


and in the absence of vaporization or chemical reaction at 
the interface 


(qui = (Goi = Gi [2a] 


For simplicity, take Le = 1 and ca = c,/2(Pr,)—*/* with 
Pr, = 0.71. 

To an observer fixed in the gas-liquid interface the molten 
material appears to be flowing toward him, as well as along the 
body surface (Fig. 1), and this countercurrent of material 
against the direction of heat flux absorbs heat energy at a rate 
given by 


6L 


Based on the simplest assumption that (pv) varies linearly 
with distance across the film, and that (07'/dy) is constant, 
this quantity has the value 0.50mzc,,(T; — Tm). Amore 
careful calculation (Appendix A) shows that the numerical 
factor is actually approximated very closely by the value 0.60 
over the whole body surface. Thus the heat flux at the surface 
gs is given by 


qs = qi 0.60m [3] 


When the heat of fusion is small compared to the heat ca- 
pacity of the solid up to the melting temperature, or Ln < 
z, (T,, — To), the melting rate approaches its steady-state 
value for a semi-infinite slab very soon after melting begins 
(Fig. 2 of (7)). Therefore 

qs 


and by Equation [3] 
mr qi/K [5] 
where # = Ln + ¢, (Tm — To) + 0.60c,, (Ti — Tm). 
In spite of the change in heat flux ¢ across the liquid film, 
the temperature rise across this layer is given by the “thin 


_ § The dynamical stability of the liquid film on a body decelerat- 
ing through the Earth’s atmosphere is analyzed by Feldman (6). 
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GAS-LIQUID 
INTERFACE 


Fig. 1 Schematic diagram of two-phase boundary layer over 
melting blunt nosed body at high velocity 


film’’ relation with a maximum error of about 20 per cent 
(Appendix A), i.e. 


Gi = (ki( Ti — [6 


Since the thermal conductivity of the liquid varies relatively 
slowly with temperature compared to the viscosity, it is suf- 
ficient at this stage to take kx, = (kx);. 

As expected, the inertia forces in the liquid film are negligi- 
bly small (Appendix A), so that 07/Oy = dp,/ds, and the shear 
stress varies linearly with distance across the film when 
dp./ds ~ 0. Although it is certainly possible to account for 
this variation, let us avoid this added complication for the 
present by taking 


= = (mr) (7] 


The final relation connecting the four unknowns, mz, ui, 7; 
and 6, is provided by the mass balance for the liquid film 


6L 
= (d/ds) prurdy) [8] 
Let 


6L 
prurdy = 


where @ = 1/2 for a linear velocity profile (for example). 
Then 

= o(d/ds) 9] 
where & = 0 for two-dimensional flow, and k = 1 for axially 
symmetric flow. 

By utilizing Equation [7] together with Equations [1 and 2], 
the film thickness 6, is expressed in terms of u;, so that 
Equation [9] furnishes a simple relation for the velocity ratio 
u:/ue. Thus 


= (Mi) [10] 
but 
(do) (Pr,) [11] 
(Tos Ue 


from which 


6, = (Pr,) [12] 
Ue qi 


By substituting this last result into Equation [9] and em- 
ploying Equation [5], one obtains 


k 
rokue(hs, any hs;) 
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Now define 
Q = 


where (q;)o is the stagnation point heat transfer rate. Re- 
gardless of conditions over the rest of the body, the gas 
boundary layer at the forward stagnation point is almost 
certainly laminar, so that (8)® 


(gio & (0.50)2"? (Pr,) W(pette)o.0 (he, — W/(due/ds)o 


and? 
ui? (pebte)o.0 hee — he; 
ue (2)? Ho 
s\f/due\ . 
ds 
and 


we 

According to Equation [14], u;/u. depends mainly on the 
driving potential (h,, — h,,,) and on the ratio of the pu product 
in the liquid and the gas. Numerical calculations of u;/u, for 
both laminar and turbulent flow are presented in the next 
section. However, a rough estimate of u;/u, for a highly vis- 
cous glass-like material can be given here. For laminar gas 
flow at the forward stagnation point, VN = 2-*. 

Suppose 


h., = 8000 Btu/lb = 20,000 fps) 

T; = 4500R 

hein >h,, = 1500 Btu/Ib 

Ho = 750 Btu/Ilb 

(pe)o.0 = po = 0.074 lb/ft? (altitude = 60,000 ft) 
PL = 131 lb/ft? 

(Me) 9.0 10-5 lb/ft see 

(Mx); => 1 lb/ft-see 

a = 1/2 


Then Equation [14] yields u;/u, = 5 X 10~4, i.e., the velocity 
ratio is negligibly small even under these rather extreme con- 
ditions, and our original assumption that the flow in the gas 
phase is virtually independent of the presence of the liquid film 
is justified. 

By employing Equations [6 and 12], one obtains 


and this relation is combined with Equation [14] to give the 
following expression for the temperature rise across the liquid 
film 


Ti — Tm V 
Pr, 


hs; ~ = hs, (pe)o.0 
i 16 


According to Equation [16], the temperature rise across the 
film depends not only on the heat transfer rate, which is con- 


® Relatively small corrections to this expression are introduced 
when Le ~ 1 (8, 9). In (9), it is also shown that ¢; depends 
weakly on the value of (pu), at the solid surface, or in this case, 
ON (Potty )i- 

* The factor (hs — hs;,/hs, — hs,) is replaced by unity in this 
approximation. 
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tained in the factor 
WV (pelle) (he, hs;) 


but also of course on the film thickness, which is proportional 
to the factor 


(The importance of the ‘driving potential” is discussed in the 
conclusions.) Once the interfacial temperature is calculated 
from Equation [16], the melting rate is obtained from the ex- 
pression 

mL (0.50)2"? x 


ole 


V (Pelle) — /H] Kdue/d9)9=0°Q 


V Rew 


[17] 


where 
Re = 


When the gaseous portion of the boundary layer is turbu- 
lent, the most important parameter for the stability of the 
liquid film is the liquid flow Reynolds number based on quan- 
tities evaluated at the liquid-gas interface. By utilizing 
Equations [12 and 13], we obtain 


_ 1 LoL (“) 
Re); kd: 18 


or 


[hee — he; 


18 


Rew 
e )g.0 


where 


For some purposes, the liquid film thickness may also be of 
interest; its value is given by the expression 


Ry 0 Ho Reo Q 


For laminar flow at the forward stagnation point, the quantity 
in brackets in Equation [19] becomes 


1 Q—k/2 
WV V Reo 


If the region of laminar gas flow persists for an appreciable 
distance downstream of the forward stagnation point, the 
expressions for u;/u, and (7; — Tm) /T'm given by Equations |14 
and 16] can be considerably simplified. According to (8) 


Qe 1 1 [20] 

Pro mae) *(1/wa)(due/ds)o 
where P = (p/po’)(ue/Uao). Now, by squaring both sides of 
Equation [20] and differentiating with respect to s, we ob- 
tain a first order equation for P, as follows (10) 
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Fig. 2 Velocity ratio and temperature at gas-liquid interface 
with laminar boundary layer at stagnation point; altitude = 
60,000 ft ; Corning Glass No. 7740; melting temperature = 3000 R 


By integration 


(du\ Q ,. 
P=—{—)¢t FO d 22 


Therefore, the quantity 


appearing in Equations [14 and 16] is approximately 2+ x 
(p/po’), for laminar gas flow. In this case we see that u;/u, and 
(T; — Tm)/T » actually decrease with distance around a convex 
blunt nosed body, while the liquid film thickness increases. 


Numerical Calculations of Flow Quantities at the 
Gas- Liquid Interface 


In order to determine the orders of magnitude of u;/u., 
(T; — Tm)/T» and (Rex);, numerical calculations are carried 
out for the following cases: (a) laminar flow at the forward 
stagnation point (from which values in the rest of the laminar 
flow region can be quickly estimated); (b) turbulent gas flow 
at the sonic point on a blunt nosed body. 


Laminar Flow at the Forward Stagnation Point 


In this case the values of u;/u, and (7; — T,,)/T, are given 
by Equations [14 and 16], respectively, with N = 2-*. Of 
course (Rex); = 0 here. The calculation procedure, which is 
utilized also in the turbulent case to be discussed in the next 
section, consists in multiplying Equation [16] by 


(ML): Ln + = To) 


and computing the new left-hand side for various values of 7’; 
inserting a value of h,; corresponding to one of these values 
of T; on the new right-hand side of the modified Equa- 
tion [16], and calculating this right-hand side for a series of 
values of h,,, or flight velocity u. at a given altitude. The 
correct value of uw. is readily obtained by graphical inter- 
polation. Calculated values of u;/u, and (T; — Tn)/T» are 
plotted against flight velocity in Fig. 2. These numerical 
results agree well with those obtained by more elaborate 
methods (1 to 4). The basic data utilized in the calculation 
are as follows: 

1 Flight altitude = 60,000 ft, or pa/po = 9.084 x 107; 
T. = 218 K, and p. = 7.232 X 10~* atm. 

2 Melting temperature of the surface material (7',) is 
3000 R; c,, = 0.33 Btu/lb/F, and L,, + (7 — To) = 750 
Btu/lb. Also, pz = 131 Ib/ft®; (kz); = 1.71 10-3 Btu/ft- 
R-sec, and (uz); = 0.0672 exp [(8720/7';;)!-*!2], where T is in 
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deg R. (These data for Corning No. 7740 glass (1).) 

3 = 0.33 T;, in Btu/Ib. Values of and (T,,),.0 
are obtained from the charts in (11), and values of (u.),,0 are 
computed by the Sutherland formula. (See also Fig. 21 
of (10).) 


Turbulent Gas Boundary Layer at Sonic Point on Body 
Surface 


So far as the temperature rise across the liquid film is con- 
cerned, the most severe conditions are encountered at or near 
the sonic point on the body surface when laminar-turbulent 
transition occurs in the gas boundary layer near the forward 
stagnation point. Because the viscosity of the liquid is 
strongly temperature dependent, the highest film Reynolds 
numbers are also attained near the sonic point, unless the: 
afterbody is unusually long. In order to predict the turbulen' 
heat transfer rate to the liquid film in this region, we must re- 
sort to semi-empirical methods, supported wherever possibk 
by qualitative physical reasoning. If the mechanism of turbu- 
lent shear flow is not much affected by enthalpy level, ther 
the turbulent heat transfer rate is expected to depend mainly 
on the total enthalpy difference across the boundary laye 
and only weakly on the mode of heat transfer (12, 13). In 
fact, the principal effect of the high temperatures appears in 
the physical properties. In the present paper we regard the 
values of these properties at the outer edge of the gas bound- 
ary layer as the most significant for the heat transfer prob- 
lem. 

If this reasoning is correct, supersonic wind tunnel measure- 
ments of turbulent heat transfer rates on blunt nosed bodies 
at low to moderate enthalpy levels can probably be extrapo- 
lated to hypersonic flight conditions with some confidence. 
Such measurements have been made by Libby and Cresci (14) 
at the Polytechnic Institute of Brooklyn Aeronautical Labora- 
tory by means of a new “shroud” technique (15) which pro- 
duces very high pressures and Reynolds numbers around the 
model. Their experiments on a spherically blunted cone show 
that the curve of heat transfer rate vs. distance along the 
body is “flat”? near the forward stagnation point, and then 
rises rapidly to a peak somewhat upstream of the sonic point. 
Of course, the location and height of this peak is strongly in- 
fluenced by the location of laminar-turbulent transition. The 
fully turbulent heat transfer rates correlate best with the local 
low-speed, nearly isothermal, flat-plate empirical relation 
cu = 0.030 Re-®.2 Pr—*/*, over a range of ratios of stagnation 
enthalpy to surface enthalpy of about 1.3 to 2.0. 

Recently Beckwith and Gallagher (16) measured turbulent 
heat transfer rates on a hemispherical nose at Mach numbers 
of 2.0 and 4.15 in the blowdown jet at NASA Langley Aero- 
nautical Laboratory. At a Reynolds number based on nose 
diameter of 3.4 X 10°, the peak turbulent heat rate is about 
2.5 times go, and this peak value is very close to the local flat- 
plate prediction. Again, gis within 20 per cent of go so long as 
0 < 20 deg. 

Shock tube data obtained by Rabinowicz (17) over a range 
3 < (h../h») < 8 and by Rose, Probstein and Adams (13) at 
enthalpy ratios equivalent to satellite velocity, support the 
conclusion that for subsonic and low supersonic Mach num- 
bers, turbulent heat transfer rates are to a remarkable degree 
independent of falling pressure gradients and the ratio of 
stagnation to surface enthalpy. This conclusion suggests the 
following procedure for the present problem of the two-phase 
boundary layer: 

1 Fora spherically blunt nose, let us suppose that g rises 
linearly with distance along the body surface from the value 
qo at s = 0 to g,* at the sonic point, which we take at 0 = 45 
deg, i.e. 

1+ 1) (x/4) 
See Sketch A. 


ARS JouRNAL 


P 
ue \ ds Jo 
| | 
« 
| € 
r 
€ 
: if 


near 
lent 


2 Suppose that the turbulent heat transfer rate at the 
sonic point is given by the usual (local) flat-plate value, i.e. 


qe* peuecu(hs, om hs;) 


where 
0.037 
[(pe*ue*s*) 


Because of the rather conservative first assumption, no correc- 
tion factor for Le ¥ 1 or h,,/h,; >> 1 is applied here. 

3 Take # = Xp in the entire subsonic region up to the 
sonic point. 

With the aid of these three working assumptions, one ob- 
tains the following expressions for the quantities required in 
evaluating u;/u, (Ti — Tn)/Tm, and (Rex); (Eqs. [14, 16 and 
1$a]) 


CH (Pr,)~*/# 


m* = = 0.18 + 0.35Q,* [23] 


s 1/4 1 du, 
Ne =| (— 
[2 06 (a*/ua) (= dé 


(Q,*(0.18 + 0.35Q,*)] [24] 


weiere 


Q* _ 0.091 p*a* 0.8 x 
But p*/pe,,9 = 0.620 + 0.007 over the range 8 S M. S 235, 


(Rew) [25] 


while 
a* 
Uo Peg.o 
and 
Uo db Peg.o 
so that 
N* = 1.11 Q,*(0.18 + 0.35 Q,*) [26] 
and 
0.050/ po 
= (Re )% {27 ] 
2 


In particular, for a nose radius of 2 ft and a flight altitude of 


60,000 ft 
0.15 0.30 
(2) (‘«) |27a] 
Pero) \108 


(Over the range 8 X 10° < uw < 22 X 10%, Q,* increases 
slowly from 3.50 to 4.31.) 

The calculation procedure is very much the same as that 
employed in the discussion of laminar flow at the forward 
stagnation point, except of course that the value of NV on the 
right-hand side of Equation [16], is given by Equations [26 
and 27], or [27a]. Values of u;/u., (Ti — Tm)/T'm, (Rex): and 
surface melting rate at the sonic point are plotted against 
flight velocity in Fig. 3. Evidently the temperature rise 
across the film is much larger than in the laminar ease (Fig. 2), 
but u;/u, is still negligibly small. A somewhat surprising 
result is the relatively low value of (Rez);. These results are 
discussed later. 


* 

~ 


Q: 


Effect of Mass Addition to Gas Boundary Layer— 
Preliminary Study 
So far in this paper, it has been assumed that no transfer of 


mass occurs at the gas-liquid interface. General Electric 
experiments on glass-impregnated materials suggested tc 
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Sketch A Turbulent heat transfer rate distribution over blunt 
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Fig. 3 Liquid film properties at gas-liquid interface and peak 

surface melting rate—turbulent gas boundary layer at sonic 

point; altitude = 60,000 ft; /?) = 2 ft; Corning Glass No. 7740; 
Tm = 3000 R 


Sutton (1) that a certain fraction of the melted substance is 
pyrolyzed into a gas that diffuses across the liquid film and 
enters the gas boundary layer. Presumably no additional 
heat transfer is required. The details of this process are 
probably quite complicated, and no way is known at present 
to determine the fraction f that is gasified, except by experi- 
ment. This phenomenon is quite different physically from 
the process of vaporization at the gas-liquid interface, which is 
uniquely determined by mass diffusion and energy transfer at 
the interface, in conjunction with the curve of equilibrium 
vapor pressure vs. temperature for the liquid. 

Actually this ‘floating parameter’ f is quite useful, because 
it furnishes a simple method of illustrating the effect of mass 
addition to the gas phase. Since u;/u.e < 1, the gaseous por- 
tion of the boundary layer behaves as if the mass addition 
occurred at a permeable solid surface at the enthalpy level 
h.({T);. In Fig. 4 the ratio ez /c, of the Stanton number with 
mass addition to the value with no mass addition is plotted 
against the injection parameter B = m;/p,uecu for two im- 
portant cases: (a) laminar flow at the stagnation point of a 
blunt body, as calculated by Reshotko and Cohen (18)8; (b) 
turbulent flow over a flat plate at M, = 2.7, as obtained ex- 
perimentally by Rubesin et al. (19, 20). (Both cases apply to 
air injection into an air boundary layer; the reduction in cy 
is greater for injected gases of higher specifie heat than air, 
e.g., lighter molecular weight.) Now in the present problem 
the rate of mass addition m; is given by 


mi = = f(qi/Ko) 
8 Reshotko and Cohen (18) utilized the injection parameter 


Vw/Vvwe = B, but B = [(Pr)w/(Nu/V Re)|B, or B = 1.165 B. 
(See Fig. 2 of (18).) 
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according to Equation [5], so that 


PeUeCHy Ho CH, 
or 
cu/cu, = [28a ] 
where 
= 
Ke 


The intersections of the straight lines represented by Equa- 
tion [28a] with the curves of cx/cuo vs. B in Fig. 4 give the 
only possible solutions. 

The reduction in heat transfer rate, or the “heat blockage,” 
can be regarded as an increase in effective heat sink capacity. 
Suppose 


mr = (Gi)o/Kest 
Then the effective heat sink capacity H.+ is defined by the 
relation 


 CHo _ (m1) B=0 


KH 


and Xers/KHy depends only on the parameter 


(he 
Ho 
1.00; T 
Ap 
EXPERIMENTAL ,RUBESIN,*™[m, = 2.7] 


9:6 
7 


NO 


9227 


2.0 3.0 
mj 


Fig. 4 Effect of mass injection on heat transfer rate (air-air) 
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Fig. 5 Effective heat sink capacity with mass addition, as a 
function of driving enthalpy potential and fraction of melted 
material evaporated 
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Once again we see the importance of the ratio of the driving 
enthalpy potential to the original heat capacity (Fig. 5). 

Although it does not seem worthwhile to go into the details 
here, the present scheme is easily modified to take account of 
the fraction of melted material that enters the gas boundary 
layer as a gas. For example, in the relation expressing the 
mass balance for the liquid film (Eq. [9]) the quantity mz is 
replaced by mz (1 —f), and (1 —f) appears on the right-hand 
side of Equation [13] for u,;?/u.?.. What is more important, 
Equation [16] for the temperature rise across the liquid filin 
now contains the mutliplier (cu/Cu)tauV 1 — f on the right- 
hand side. In both Equations [17 and 18a], the simple factor 
(cz/Cu,)Lam appears, so that Kp is replaced by Herr. For a 
given arbitrary value of f, the calculations proceed much zs 
described in the discussion of the numerical calculations, ex- 
cept that for every value of h,, (with a fixed 7';) the value cf 
(Cu/Cu,)Lam is obtained from Fig. 5. When laminar-turbulent 
transition occurs near the forward stagnation point the valu: 
of Q,* (Eq. [25]) is increased by the factor 


/CHy)LAM 


Evidently mass addition by this process can have an impor- 
tant effect in reducing (7; — 7)/Tm and (Rez);. 

When true vaporization takes place at the gas-liquid inter- 
face, the situation is somewhat more complicated in genera! 
(3, 5). However, when 7’; is close to the local boiling tem- 
perature, f > 1, ie., practically all of the molten liquid 
vaporizes. In this limiting case the problem is again con- 
siderably simplified, and the main features can be distin- 
guished without considering all the details. When Lyap > Ho 
a similar simplification occurs even for relatively small 
values of f, as in the case of HO; here Lvap is the heat of 
vaporization. In both cases 9; 2 mi(Lvap + Ho), so that 


tin; — 


= 
(Lvap + Ho) 
also 
CH B 


Now at hypersonic flight velocity, h,, is 3 to 10 times larger 
than h,;, and both (h,, — hs;) and (3o + Lvap) are relatively 
insensitive to the value of the interfacial temperature 7’,, 
which is unknown initially. Therefore B’ can be determined 
rather accurately by inserting a reasonable first estimate for 
T;, and the values of B and m; are then obtained as before 
from the curves of cx#/cx, vs. B in Fig. 4. For Le = 1, the 
analysis of (12) shows that 
(Kei = Brey 

where (Kz); is the element mass fraction of the vaporizing 
material at the gas-liquid interface. In the absence of chemi- 
cal reactions involving the vaporizing material, the species 
mass fraction Kz = Kg, and 


PE 1 Me | 


where pz is the partial pressure of the vapor, and Mg and Mair 
are the molecular weights of vapor and air, respectively. Thus 
a better approximation to the interfacial temperature is de- 
termined by means of these last two relations, and a second 
iteration usually suffices. 

These simple relations furnish a rapid estimate of the flow 
regime for any particular problem; in the “boiling” regime 
they show that 7; “adjusts itself’ to a value determined 
mainly by the energy balance and the “blockage” effect of 
mass addition on heat transfer. Chemical reactions can also 
be included in this scheme in the manner described in (12). 
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Discussion and Conclusions 


1 By employing a simple “model” for the two-phase 
boundary layer produced by surface melting, simple algebraic 
expressions are obtained for the velocity ratio u;/u, at the in- 
terface, the temperature rise across the liquid film, and the film 
thickness and Reynolds number. These relations contain the 
“sealing laws’ for the regime in which no vaporization takes 
place at the gas-liquid interface. 

Even under the most extreme conditions, one finds that 
u./ue < 1 for fairly viscous liquids, so that the original as- 
sumption that the gas boundary layer is virtually independent 
o! the liquid film is justified. 

2 Because the liquid viscosity is strongly temperature de- 
pendent, the temperature rise across the liquid film is an im- 
portant indicator not only for the onset of vaporization, but 
also for the film Reynolds number, or liquid film stability 
when the gas boundary layer is turbulent. This temperature 
rise (7; — Tm)/Tm depends on the product of the local heat 
transfer rate and the factor 


V (Ing — V 


where (h., — h.;) is the driving enthalpy potential, and 5p is 
the heat capacity of the melting substance up to the in- 
terfacial temperature. Thus, it is not sufficient to duplicate 
the local heat transfer rate alone in model experiments. For 
example, in tests made in rocket exhaust jets the ratio 


(hs, h,,)/Ho 


is generally of the order of unity, as compared with a value of 
6 to 8 ata flight velocity 20,000 fps. (7; — Tm)/T'm, in other 
words, might not be high enough in such experiments to 
produce vaporization at the gas-liquid interface, or the liquid 
film Reynolds number might be too low to bring out film in- 
stability that could occur under flight conditions. 

3 In free flight, the temperature rise across the liquid film 
is proportional to Vv Po Uo at the forward stagnation point, 
and to po%® uo** (approximately) at the sonic point on the 
body surface when laminar-turbulent transition occurs near 
the stagnation point. Thus, the critical conditions for the 
liquid film can be determined for any given flight trajectory. 
For ballistic missile re-entry, for example, the well-known 
Allen-Eggers analysis could be applied to this problem. 

In order to illustrate the properties of the liquid film for 
typical flight conditions, calculations were made for a flight 
altitude of 60,000 ft. The material selected is Corning glass 
No. 7740, which melts at 3000 R. As expected, the value of 
(7; — T»)/Tm is much smaller at the forward stagnation point 
with a laminar gas boundary layer than it is at the sonic point 
on the body surface when laminar-turbulent transition occurs 
near the nose (Figs. 2 and 3). For example, the interfacial 
temperature 7’; is 5000 R at the sonic point when the flight 
velocity u,, is 15,350 fps, but 7; does not reach this value 
at the forward stagnation point until u,, = 19,600 fps. To 
put it another way, when u,, = 15,350 fps, the value of 7; at 
the stagnation point is 4100 R, or (7; — T’,,) is about one-half 
the value at the sonic point. Vaporization is unlikely at the 
forward stagnation point, but quite possible at the sonic 
point when transition occurs near the nose. 

4 A somewhat surprising result is the low value of the 
liquid film Reynolds number (Fig. 3). At a flight velocity of 
15,350 fps, (Rex); = 2.5 at the sonic point. In other words, 
the liquid film is quite stable; in fact, vaporization would 
take place before instability could occur, and this mass loss 
should hold the film thickness and Reynolds number to low 
values. 

5 About 750 Btu/Ib is absorbed by this material up to the 
melting temperature. But the countercurrent of molten 
material toward the interface (and against the direction of 
heat flux) absorbs an additional amount of energy equivalent 
to that required to heat about 60 per cent of the melt from 
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the melting temperature up to the interfacial temperature. 
At a flight velocity of 15,350 fps, this energy increment 
amounts to about 375 Btu/Ib, and the total 3¢ = 1125 Btu lb. 
The melting rate at the sonic station is 1.30 Ib/ft?/see. This 
material has about six times the heat capacity of a nonmelting 
copper heat sink. 

6 A study of the effect of mass addition to the gas phase 
shows that the amount of heat “blocked” by this process is 
sensitive to the ratio 


(Reg — he;)/Ho 


to the state of the gas boundary layer, and, of course, to the 
fraction of melted material that is gasified or evaporated. 
For example, in the process considered by Sutton (1), a cer- 
tain fraction (f) of the melted pyrex glass is pyrolyzed into a 
gas and enters the gas boundary layer. With f = 0.08, the 
increase in “effective” heat capacity at a flight Mach number 
of 18 is about 30 per cent at the forward stagnation point, but 
only about 10 per cent for a turbulent gas boundary layer at 
the sonic point. On the other hand, in experiments in a 
typical rocket exhaust the effect at the forward stagnation 
point is reduced to 10 per cent, while at the sonic point it is 
negligible for turbulent flow. 

More spectacular results are expected when true vaporiza- 
tion occurs at the gas-liquid interface. When 7; approaches 
the local “boiling temperature,” f — 1, and the problem is 
considerably simplified. For the same flight conditions con- 
sidered by Sutton (1), the effective heat capacity of the pyrex 
material is about doubled to a value of 2500 Btu/Ib in the tur- 
bulent case, while at the forward stagnation point Ren = 
5000 Btu/Ib (Fig. 5). It would seem desirable to search for (or 
invent) a material that is fairly viscous in the liquid state, 
with a boiling temperature not too far above the melting 
temperature. 


Appendix A 


Estimate of Magnitude of Inertia Forces in Liquid Film 


The change in shear stress across the liquid film required to 
provide the molten material with its momentum parallel to 
the surface is given by the expression 


6L 
Ar = f, prr(dusdy)dy = (piv) u; 


so that 


A us 


Ti Ts dite 


by Equation [11] of the text. 
But (piv) =", = gi/%, so that 


Now u;/u, is about 5 to 7 X 1074 at most, while 


= hy 


Il? 
| 
a 


is 7 to 8 at most; therefore, Ar/7; S 6 X 107%, and the cor- 
rection to 7; from this source is less than 1 per cent. A similar 
calculation shows that the term 6,p1,u(0u/Os) is also of this 
order of magnitude compared with 7; near the forward stag- 
nation point, but is generally much smaller on a conical after- 
body. 


Energy Balance for the Liquid Film 


Because the thermal diffusivity of the liquid is small com- 
pared to its kinematic viscosity, the “thermal inertia” terms 
like pxve,,(O7T/dy) must be taken into account. In fact the 
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contribution of this term to the relative change in heat flux 
Ag/@ across the film is larger than the corresponding contribu- 
tion of pxv(du/dy) to Ar/r by a factor of the order of the 
Prandtl number. At a temperature of 4000 R, Pr, = 445, for 
a typical highly viscous glass; therefore, a more detailed 
calculation is required. On the other hand, the contribution 
of the viscous dissipation to Ag/g is of order u;/u., and the 
contribution of the term u(dp/ds) is also of this order. 

In a region where the gas phase boundary layer is either 
wholly laminar or wholly turbulent, and the pressure is nearly 
constant over the body, the variation in the interfacial tem- 
perature 7’; along the surface is negligibly small. However, 
when the pressure is dropping rapidly, as on a blunt nose, or 
in a region of laminar-turbulent transition in the gas phase, 
the term prucp,(0T/0s) cannot be ignored. Roughly, the 
contribution of this term to Ag/g is proportional to the quan- 
tity 

cp T ip) /Ho 


where 7’,, is the interfacial temperature at the forward stag- 
nation point. This term therefore opposes the main contribu- 
tion from pzvc,(OT/dy) over a blunt nose, except in the 
region of laminar-turbulent transition in the gas phase, where 
T; increases with distance along the surface. It appears that 
the contribution made by this term is about 5 to 20 per cent 
of the main effect, at most. Therefore, it seems justified to 
retain only the principal term in the heat conduction equation 
in order to obtain a useful estimate of the change in heat flux 
across the liquid film. Thus 


oT oT 
\ kt dy? [A-1] 


By integration, one obtains 


log. (4) [A-2] 


and some estimate of the lateral velocity is required. 
According to the equation of continuity 


1 re) 
2) [A-3] 


Suppose that u/u; = 1 — (y/6x1). Then 
(piv) = Ay? — By [A-4] 


( ) ( 
2 To Os 6: 
1 ra) 
B= (‘)(2) (pruiro) 


By substituting this last result into Equation [A-2], we ob- 


0/dy (piv) = — 


where 


and 


tain 
_ (Ay? _ By? 
or 
By* 
qi ky 2 


Cp, By? 


The value of g,/gi is obtained by substituting y = 6, into 
Equation |A-6]. By integrating Equation [A-6] once more, 


we obtain 


The coefficients A and B are readily evaluated in two im- 
portant special cases: (a) near the forward stagnation point; 


(Ti — Tm) = 
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(b) in a region of constant pressure on a conical afterbody. 
In the first case, 6z is constant, and 


by Equation [9] of the text, with a = 1/2. Similarly, B = 


2mz./61. Therefore, 
qs . 2 Cry, 
—~=]— —*§ 5d, A-8 
and 


By Equation [A-9] 


kt 


For our purposes, one cycle of this recursion formula is suf- 
ficient, so that 


qi 


Thus in the most extreme case (Figs. 2 and 3), where 


1 ep (Ti [A-11] 
4 


ep (T: — Tm)/K [08 
the simple “thin-film” relation 
Gi = [ku(Ti — 


is accurate to within 20 per cent; in most cases the error is 


much smaller. 
From Equation [A-8], one finds that 


9 
Gs = 3 [ep + 


18 | [ep,mi(T; — Tm)] [A-12] 


so that the relation 
ds = Gi 0.60cp Ti = Ts) [A-13] 


is sufficiently accurate for the present analysis. In other 
words, the energy absorbed across the film is equivalent to 
that required to heat about 60 per cent of the melted material 
up to the interfacial temperature. 

A similar calculation can be carried out for a region of con- 
stant pressure over a conical afterbody. In this case the con- 
stant 1/4 appearing in Equation [A-11] is replaced by 7/36, 
while the coefficients 2/3 and 1/18 in Equation [A-12] are re- 
placed by 5/9 and 5/108, respectively. In the computations. 
the relation given by Equation [A-13] is used both at the for- 
ward stagnation point and at the sonic point on the body 
surface. 


Ty, m i 


Sketch B Temperature gradient as a function of temperature 
across liquid film 
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Appendix B 


The approximate expressions obtained for the velocity at 
the gas-liquid interface and the temperature rise across the 
liquid film (Eqs. [14 and 16]) can be given a simple interpreta- 
tion by means of a polynomial representation of the tempera- 
ture gradient (—d7T'/dy) = w as a function of temperature 
(see Sketch B). 

u(T) = b(7' — To) 
while for 7, S T S T; 
3 
WT) = >> aT — T,)" 
n=0 


where 


and 
a = Ui 


At the junction point of these two regions, we require that 
dw/dT and d*w/dT? are continuous, or that 7’, v and u are all 
continuous, and that 


Now from Equation |A-1] 


[B-1] 


while the equation of continuity yields the following relation 
ut the forward stagnation point 


i kiw 
dy Je 


Cp, aT? 


Also the relation = y(O0u/dy) yields 


d dw Aor fdue\ Cr, 
dT (« Ue ) ke BL 
We recognize that at T = 7; the value of 7 in Equation [B-3] 


is equal to 
= 
In the absence of vaporization 
(dw/dT) 7, = 0 
and the left-hand side of Equation [B-3] at T = 7; is simply 
w;? (d’w/dT®); 


By evaluating w;, (d?w/dT?)7; and (d’w/dT*)7;, recognizing 
that 6 = —(cp,/kx)pi@, and utilizing the three conditions 
that dw/dT and d*w/dT? are continuous, and that 

Aw = 


at T = T,,, we are able to solve for the three unknowns 7; — 
Tm, Ui/Ue, Pi¥@; the results are identical with the expressions 
given in Equations [15 and 16], with a = 1/2, except that 


Ho = Lan Cy (T'm = + (2/3)ep,( T; 


Nomenclature 

a = local sound speed 

B = ™;/petleCuo), Mass addition parameter 
B’ = Mi/peleCH, Mass addition parameter 
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cs = local skin friction coefficient, (t,)/(*/2peue?) 

CH = local Stanton number, gg/psuchs, — hs;) 

Cy = specific heat at constant pressure 

Ce = average value of specific heat of solid material up to 
melting point 

ni = fraction of melted material gasified or evaporated 

= f: (hse — hs;)/(Ko) 

hs = stagnation enthalpy 

KH = heat capacity of liquid-solid system, 
Lm + —To) + 0.60c, (7: — Tm) 

Hert = effective heat capacity 

k = thermal conductivity, also geometric index 

Ln = heat of fusion 

mi = rate of mass addition to gas phase 

ML = melting rate 

M = Mach number, Vue + v?/a 

mM = average value of normalized melting rate (Eq. [14]) 

N = s/udu,/ds))Qm. (See Eq. [14].) 

p = pressure 

Po" = stagnation pressure behind a normal shock 

Pr = Prandtl number, cpyu/k 

q = local heat transfer rate 

Q = 

Q = average value of Q, (1/roés) £ Q rds 

ro = radius of cross section of body of revolution 

Ro = nose radius 

(Rez); = liquid film Reynolds number, 

Reo = 

Rew = U 

s = distance along gas-liquid interface, measured from for- 
ward stagnation point (Fig. 1) 

T = absolute temperature 

Tm = melting temperature 

To = “internal’’ solid temperature 

u, v= components of velocity parallel and normal to gas- 
liquid interface (Fig. 1) 

y = distance normal to interface (Fig. 1) 

a = velocity profile parameter for liquid film, ” pxuatly = 

6 = boundary layer thickness 

m = viscosity 

v = kinematic viscosity, u/p 

p = density 

0 = angle between flight direction and radius vector from 
center of curvature of nose 

T = shear stress 

Subscripts 

e = flow quantities evaluated at outer edge of gas bound- 
ary layer 

g = gas 

i = gas-liquid interface 

L = liquid 

0 = stagnation point (except 7»), also zero mass addition 
rate 

s = body surface 

t = turbulent flow 

co = free stream ahead of bow shock wave 

Superscript 


(*) denotes physical quantities evaluated at sonic point on 
body surface 
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Some Physical Properties for the System 


Nitrogen Tetroxide-Nitric Oxide 


DONAT B. BRICE? and 
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Mixed oxides of nitrogen may become an important oxidizer with the advent of storable, pre- 
packaged, liquid propellant rocket engines. This paper presents some physical property data for 
the liquid system nitrogen tetroxide-nitric oxide at temperatures between zero and minus 100 F. 
Properties included are vapor pressure, density, viscosity and surface tension. 


APID progress in the development of solid propellants 

has threatened to force into obsolescence liquid fueled 
systems with their notorious complex arrays of hardware and 
their unreliability. However, recent advances in the pre- 
packaging of liquid propellants, with the resulting effect of 
instant readiness previously reserved for solids, have brought 
about a resurgence of interest in liquids, with particular 
emphasis on storable systems. It is thought to be timely, 
therefore, to present physical property data for MON (mixed 
oxides of nitrogen) obtained in 1953 during the course of 
studies concerned with the production of this potential liquid 
oxidant. 

No experimental physical property data for vapor pressure, 
density, viscosity and surface tension were found in the 
literature for mixtures of nitrogen tetroxide-nitric oxide 
(mixed oxides) containing 20 to 30 weight per cent nitric 
oxide at temperatures between 0 and —100F. Such informa- 
tion is desirable in the design of propulsion systems utilizing 
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mixed oxides as the oxidizer. The information is also neces- 
sary in the design of equipment for the manufacture of mixed 
oxides. 

Whittaker et al. (1)4 have determined the vapor pressure of 
mixed oxides for mixtures containing a maximum of 16.85 
weight per cent nitric oxide at temperatures as low as —40 F. 
Sage and co-workers have determined vapor pressure (2) and 
viscosity (3) of mixtures of nitrogen tetroxide and _ nitric 
oxide. Most of these latter data were measured at tempera- 
tures above 40 F and in the range of 5 to 25 weight per cent 
nitric oxide. Results of investigations presented in this paper 
appear to be in substantial agreement with extrapolated 
values of the referenced data. 


Equipment and Results 


Materials 


The mixtures of nitrogen tetroxide-nitric oxide used for 
these determinations were prepared by condensing gaseous 
mixtures of commercial grade nitric oxide and nitrogen tetrox- 
ide, distilled from commercial grade liquid. The gases were 
flowed in the required proportions to obtain the desired com- 
positions of the mixed oxides. The exact percentages of 
nitric oxide in the mixtures were determined by freezing point 
techniques, using the data of Fig. 1. 


* Numbers in parentheses indicate References at end of paper. 
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Fig. 1 Freezing point curve of nitric oxide-nitrogen dioxide 
mixture vs. weight per cent of nitric oxide 


VA 
7) 
| 
MERCURY = 
MANOMETER. 
=a 
e 


STIRRER 


STN STL 
RESERVOIR 
ACE TONE- 


ORY ICE 
BATH 


SWAY STOPCOCK 
TO 
VACUUM PUMP 


Fig. 2 Apparatus for measurement of vapor pressure 
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Fig. 3 Vapor pressure of mixed oxides 
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Vapor Pressure 


For the vapor pressure measurement, the sample was con- 
tained in a stainless steel cylinder. As can be seen from Fig. 2, 
a manometer was connected to the cylinder through a suitable 
valving arrangement. After the lines were evacuated, the 
cylinder was transferred from storage under liquid nitrogen 
to an acetone-dry ice bath and thawed at the lowest tempera- 
ture at which vapor pressure was to be measured. The liquid 
material was allowed to come to thermal equilibrium at each 
successively higher temperature prior to reading the vapor 
pressure indicated by the manometer. The data are pre- 
sented in Fig. 3 as a plot of the log of the vapor pressure vs. 
the reciprocal of the absolute temperature. 
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Fig. 5 Density of mixed oxides 
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Fig. 6 Isotherms of density for mixed oxides 


It can be seen from the shape of the curves of Fig. 3 that 
the log of the vapor pressure deviates from the straight line 
relationship with reciprocal temperature as the liquid tem- 
perature approaches the freezing point. This is explained by 
the possibility that at temperatures near freezing, nitrogen 
trioxide (N,Os) is more stable as a compound, and therefore 
less pressure is exerted. As the temperature is raised, the 
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stability of NO; decreases, and nitric oxide is released to the 
status of a solute and can exert increased pressure. 


Density 


The density of mixed oxides was determined by weighing 
« glass hydrometer float of known volume and mass while im- 
mersed in the sample. As shown in Fig. 4, an acetone-dry ice 
bath was used to obtain the desired temperature. The ap- 
parent weight of the hydrometer was used to calculate the 
density of the mixed oxides. Determinations, which were 
made over a range of temperatures from 0 to —100 F and a 
range of concentrations from 21.4 to 29.6 weight per cent 
nitric oxide, are shown in Fig. 5. Fig. 6, a crossplot of the 
data presented in Fig. 5, shows density isotherms. 


Viscosity 


The viscosities of the mixed oxides were measured with an 
Ostwald viscosimeter, an instrument used to compare the 
iscosities of two liquids. If the coefficient of one is known, 
he coefficient of the other can be determined. This method 
iminates the tedious experimental work involved in the 
lireet measurement of viscosity, which requires the deter- 
nination of the rate of flow of a liquid through a capillary 
‘ube of known dimensions. Comparison was made with di- 
‘thy! ether measured over the same temperature range using 
viscosity data from the International Critical Tables (4). 
Fig. 7 shows the viscosity of mixed oxides in centipoises 
plotted against temperature. It can be seen from these data 
that the viscosities of mixed oxides vary only slightly with 
composition between 20 and 30 weight per cent nitric oxide. 


Surface Tension 


The surface tension determinations were made using a Du 
Noiiy ring tensiometer (Fig. 8) and thermostatting the sample 
for each measurement at each temperature. The ring was 
suspended by a platinum wire through a small opening in the 
cap of the sample flask. Results of the measurements for the 
four compositions of mixed oxides studied are presented in 
Fig. 9. 

Difficulties in handling contributed to the inconsistency of 
the surface tension data. No means were provided for mixing 
of the sample to insure constant surface composition. Also 
despite precautions taken, atmospheric moisture condensed 
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Fig. 9 Surface tension of mixed oxides 


on the inside wall of the sample flask making visibility poor 
and causing further uncertainties in the surface composition. 
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Technical Notes 


Thermostatic Temperature Control of 
Satellites and Space Vehicles 


R. A. HANEL! 


U. S. Army Signal Research and Development Laboratory, 
Fort Monmouth, N. J. 


The temperature of satellites and space vehicles varies 
considerably with orbital conditions. To achieve greater 
reliability and efficiency in long-life instruments, an auto- 
matic temperature control is highly desirable; a radiation 
thermostat performing this function is discussed. The 
suggested method regulates temperature by adjusting the 
effective absorptivity and emissivity. Bimetallic strips 
automatically move a light shield which exposes surfaces 
with high and low absorptivitiy-emissivity ratios. No in- 
ternal power source is needed. Calculations demonstrate 
the effectiveness of the radiation thermostat. The tem- 
perature of the satellite’s instruments can be kept constant 
within a few degrees of the design value, regardless of orbi- 
tal conditions, internally dissipated power and some ero- 
sion of skin coatings. 


HE TEMPERATURE of satellites and space vehicles, or 

of any astronomical body, is determined by the thermal 
radiation properties of the surface and the amount of in- 
ternally dissipated power. Absorptivity and emissivity in 
the visible and infrared region of the spectrum balance the in- 
coming and outgoing radiation; the dissipated power causes 
a temperature rise until a new equilibrium is reached. 

Incoming radiation consists of direct and, when near the 
Earth or another body, reflected solar radiation; thermal 
radiation from the Earth, or the moon in case of a moon 
satellite; and radiation from the stellar background. The 
amount of outgoing radiation is determined by the emissivity 
and temperature of the vehicle’s surface. 

The incoming energy flux is a function of time; it may vary 
considerably with orbital conditions. Thermal emissivity of 
the Earth or moon, albedo and the orbit are known or can be 
predicted with limited accuracy only. These and other un- 
certainties, such as the aspect ratio for nonspherical bodies, 
limit the accuracy of thermal calculations. The temperature 
of an orbiting body can be predicted within a certain range 
only, and all instruments and circuits have to function prop- 
erly in this range. For example, Project Vanguard’s instru- 
mentation must operate between 0 and 60 C (1).? 

Satellites with a longer instrument life will have to contend 
with seasonal changes. An Earth satellite in a near-polar orbit 
might face the sun 60 per cent of each orbit after launching. 
Within 2 or 3 months, it will drift into a 100 per cent sun-il- 
luminated orbit. Since even the temperature limits of 0 and 
60 C create difficulties for batteries and power transistors, 
long-life satellites without a temperature control face difficult 
design problems. Design, construction and test time would 
be reduced considerably if temperature requirements could be 
relaxed. Reliability of long-life instruments would certainly 
improve if the temperature could be adjusted and kept at a 
predetermined level. 


Received July 14, 1958. 
1 Project Scientist. 
? Numbers in parentheses indicate References at end of paper. 


Changing the internally dissipated power is an inefficient 
way of controlling the temperature; a better method is adjust- 
ing absorptivity or emissivity according to the requirements. 
Coatings which change radiation properties within a rather 
narrow temperature range centered around 25 C would be 
ideal, but such coatings are not available. Heat-sensitive 
paints that change color at certain temperatures have been 
known for some time and may serve as a starting point for 
research in this field. 


Principle of Radiation Thermostat 


The solution presented here uses two coatings with opposite 
radiation properties. The first one has a high a/e ratio; the 
Tabor (2) surface with an absorptivity for solar radiation 
of 0.9 and an emissivity of 0.1 at room temperature is suited 
for this purpose. The second coating has an extremely low 
a/e; lead carbonate white paint would be a good choice, since 
a is about 0.05 and e is 0.95. The Tabor surface would give a 
high and the white paint a low equilibrium temperature (3). 
In the recommended configuration, the shell is painted with 
the two different coatings in a zebra-like fashion. 

A perforated hood of highly reflective material can be ad- 
justed so that at times the Tabor surface is completely cov- 
ered and the white paint determines the heat balance, and at 
other times the perforated hood covers the white paint and 
exposes the Tabor surface. The adjustment of the hood is 
accomplished automatically by a bimetallic strip attached to 
the instrumentation can or to the shell. If the satellite enters 
a cooler orbit, the temperature of the can will drop. The bi- 
metallic strips, in good thermal contact with the can, will 
rotate the perforated hood exposing more of the ‘“‘black”’ sur- 
face, which in turn leads to more absorption. If it enters a 
hotter orbit, absorption is reduced. This is an example of 
negative feedback in a thermomechanical system. 


Feasibility of Thermostat 


The feasibility of such a thermostat will be demonstrated 
here by calculating the temperature of an Earth satellite in- 
strumentation similar to the cloud cover experiment of 
USASRDL (4, 5). However, all calculations and conclusions 
are applicable to larger, heavier and more sophisticated in- 
struments as well. 

In these considerations, a spherical or short cylindrical shell 
is assumed to be in good thermal contact with the specially 
coated surfaces, adjustable hood, light baffles for the optic, 
spin stabilizer, solar batteries, antennas and other devices. 
The shell encloses the instrumentation can which contains all 
electronic circuitry, instruments and batteries. In good ther- 
mal contact with the can are bimetallic strips and the optics. 
Shell and can are thermally isolated from each other as much 
as practicably possible, and they can be treated thermally as 
two bodies as shown in the energy flow pattern (Fig. 1). 

Radiation into the optic and low temperature radiation from 
the stellar background are small, and in the following con- 
siderations are neglected. 


Calculation of Energy Flow 


The difference between incoming and outgoing energy flux 
will cause the temperature to rise or fall, depending on which 


3 The cloud cover experiment, scheduled in the IGY Program 
as Group II on the Vanguard vehicle, will detect and track cloud 
formations including tropical storms, such as hurricanes and 
typhoons. 
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Fig. 1 Energy flow pattern 


term preponderates. This is expressed by two equations, one 
for the shell and one for the can 


dT, 
Rsun + Rreti + Rearth + Rive + Reond — Reat = dt 


dT. 


W — R.. — Reoud = Mele = 


If the optic is treated separately, a third equation is required. 
Hach term is analyzed separately 


1 
AOS [as + a2A2(1 k) + + assA,] [2] 


where f,(t) is essentially a step function which indicates the 
amount of time the satellite is in the sun (f; = 1) and in the 
shadow (f; = 0). S is the solar constant: 2cal/min cm? = 
0.14 watt/em*. and are the absorption co- 
efficients for solar radiation for 1) the Tabor surface, 2) the 
white paint, 3) the hood, 4) the rest of the body. A indicates 
the corresponding areas and k the position of the hood in rela- 
tion to the black and white strips; & will be made a function of 
can temperature T. by the bimetallic strips. The factor } 
expresses the ratio between the cross section of a sphere and 
its total area. 


Q 
Rrett fo(t) SAlb[a,,Aik + = k) + + 
[3] 


Q is the solid angle under which the Earth appears. It varies 
with altitude. The function f.(t) can be approximated by a 
halfwave rectified cosine function. Its exact value depends 
on altitude and orbital conditions. Alb is the average albedo 
of our planet. It can vary between 0.2 and 0.5 depending on 
cloud cover and geographical latitude. A good average value 
is 0.35. 


Q 
Rearth = oT Alan Aik + — k) + [4] 


Thermal radiation of the Earth is on the average equivalent 
to the radiation of a black body of 250 to 260 K; a is the 
Stefan-Boltzmann radiation constant: 5.67 X 10~!? watts per 
cm? deg‘. a; indicates the long wave absorptivity and emis- 
sivity values which correspond to a, as explained under 
Equation [2]. 


Ry = ae) [5] 


Radiation between shell and can be determined by the area of 
the can, thermal emissivity and absorptivity of the inside of 
the shell and the outside of the can, and the difference of the 
fourth powers of the corresponding temperatures. 


Reond = p(T. [6] 


Energy flux due to conductivity in construction elements be- 
tween can and shell is proportional to the temperature dif- 
ference. The factor p expresses specific conductivity of sup- 
port materials and geometric factors. In general this term 
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can be kept smaller than R,., the previous term. 
Rat = oT Aik + apA(1 — k) + ayAs + [7] 


Reradiation of the satellite is determined by shell temperature 
T. and long wave emissivity. The terms m,c,(dT,/dt), 
m.c-(dT./dt) represent the thermal energy which is used to 
heat the shell and can. m,, m. represent mass of shell and can, 
respectively; c,, are the specific heats, and (dT7’,/dt), 
(dT./dt) are temperature changes per time unit. W is the in- 
ternally dissipated electrical power. In general, W is also a 
function of time. However, for this calculation a constant 
average value is assumed. Introducing these values, two 
equations are obtained: 
Shell 


1 Q 
( fi + fe S[as:Aik + — k) + as3A3 + + 


Q 
o (2 rs) {an A,k + k) + a23A3 + 


dT: 
= [8] 


Can 


—oa,A(T = p(T. Ts) + W= MeCe 
This set of nonlinear nonhomogeneous simultaneous differen- 
tial equations for 7’, and 7. is not solvable in closed form. 

However it can be expected that the can temperature will 
fluctuate around an equilibrium value. Exponentially de- 
creasing terms will satisfy initial conditions. If we assume 
that the initial temperature at launching 7» is the desired 
equilibrium temperature, the exponential terms will become 
zero, and the can temperature is completely described by 
equilibrium temperature 7) and a periodic time-dependent 
function 7; (t) 


T. = + [9] 
Similarly the shell temperature is 
T, = T. — AT + Tit) 
where AT’ represents the average temperature difference be- 
tween can and shell. The fourth-power terms of 7. and T, 
will be expanded in a Taylor series and only the linear term 
taken into account. This is justified, since AT and the 
fluctuating terms 7 and 7» are expected to be small with re- 
spect to the absolute temperature level 75 
To > | AT |; [71]; | 72! 
T(t) A Tot + 
T.At) To! +47,3[—AT + T,(t)] 


Introducing these values, the equations become linear in 
the time-dependent terms, and separation of fluctuating and 
nonfluctuating parts is possible. For the can 


—I[p + AT + W=0 


dT 
—[p + 3] [Ti(t) — Tx(t)] = mC. 


The time-independent equation gives the average temperature 
difference 


W 
[12] 
p + a,.A.o4T,3 


The AT varies between a few degrees and 10 to 15 C, depend- 
ing on the dissipated power and isolation achievable. A nor- 
mal value of 5 C may be assumed. 
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Temperature Fluctuations of Can 


Temperature fluctuations of the can 7; may be calculated 
as a function of fluctuations in the shell T, by expanding each 
term in a Fourier series. Thus 


= > ane T, = [13] 
n n 
Introducing this in Equation [11] yields 


1 
Ba = On 1 — + p) 


which can be approximated by 
+ p 


Bn = 
inWM Ce 
Therefore 
T, = On einut (14] 


In practice, temperature fluctuations of the can could be kept 
below 1 per cent of the shell’s fluctuations. The can’s 
fluctuations lag approximately 90 deg behind the shell’s. 


Calculation of Average Temperature 


Following the same procedure, the time-independent equa- 
tion for the shell is 


G fit fe ab) SlanAik + — k) + + + 


o (2 TA —Tott+ {15] 
{anAik + — k) + + ands] + W = 0 


The bar on f; and f, indicates average values. 

In order to demonstrate the stabilizing action of the system 
numerically, let us assume the numerical values: A = total 
surface area of satellite 


ay = 0.9 ay = 0.1 A, = 0.25A Tabor surface 

a@2= 0.05 apn = 0.95 A, = 0.25A_ white paint 

ay, = 0.12 A; = 0.25A polished metal hood 

au = 0.418 A, = 0.5A © solar batteries and 
rest of housing 


8 

2 & 

o 


Absorptivity and emissivity for the rest of the housing (a,, and 
a;,) are chosen so that, for an average orbit, k becomes ap- 
proximately 0.5, which is the center of the regulation range. 
This condition can be achieved easily by proper coating of 
areas between solar cells. The term 


1 Q 
i 2 am) s=2 


which is mainly responsible for the incoming radiation is the 
most uncertain one. The Z may vary from 0.021 for a “‘cold’’ 
orbit (60 per cent sun, low albedo), to 0.04 for a “‘hot’’ orbit 
(100 per cent sun, high albedo). 

Solving for the average can temperature JT 


Q W 
Z 


In Fig. 2, 7» is plotted as a function of Z, showing k as a 
parameter. In the numerical calculation T, was assumed 
to be constant and equal to 255 K. AT has been taken equal 
to 5 C, W = 1 watt, and A = 8000 cm?, equivalent to the 
surface of a 20-in. sphere. Coldest and hottest orbits are in- 
dicated in the figure. 

Design limits, however, have to reach farther in order to 
cover uncertainties in numerical assumptions, tolerances in 
the effect of surface coatings, aging, and other possible changes 
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Fig. 2 Temperature as a function of incoming radiation for 
several values of k. Dashed line shows regulation of thermostat 
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Fig. 3 Energy flow pattern with fast acting feedback loop from 
shell 


of emissivity and absorptivity that may occur due to erosion 
in orbit. Design limits equivalent to Z = 0.015 and 0.05 have 
been chosen arbitrarily. 


Thermal Feedback 


In a satellite without a temperature control where k is 
constant, it can be seen from Fig. 2 that the can temperature 
can vary 60 C. By making k a function of temperature, this 
range is reduced. Let k become 


k = 0.5 + 0.1(300 — 7») [17] 


a condition easily achievable by selecting proper bimetallic 
strips. The adjustable hood is rotated from one extreme posi- 
tion (k = 0) to the other one (k = 1) within a temperature 
range of 1OC. This is shown in Fig. 2 by a dashed line which 
also indicates the can temperature for all orbital conditions. 

Regardless of how “cold” or “hot’’ an orbit may turn out 
to be, all sensitive instruments will stabilize within +5 C of 
the design temperature. If necessary, this range could be 
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further reduced by increasing the sensitivity of the bimetallic 
strips. 


Temperature Fluctuation of Shell 


Temperature fluctuations of the shell can be determined by 
solving the time-dependent part of Equation [8] 


(| Sit fe SlaiAifi(k)] — 407 — 


dT, 
+ p)Tx(t) = mcs {18] 
or 
— + = me. — 


In f; and fo, only the fluctuating components will be taken. 

The temperature fluctuations of the can 7; are neglected, 
since | T7;|<|7.|. Expanding Z and 7; in a Fourier series 
again 


Z= 


n 


= > a,einut 
n 
[asiA ifi(k)] 
W 


\ssuming reasonable numerical values, calculations show tem- 
perature fluctuations of the shell on the order of 20 to 50 C. 

These fluctuations are not influenced appreciably by the 
feedback loop, since fluctuations of the can are only a fraction 
of the fluctuation of the shell, and therefore k will be prac- 
tically constant in the time interval of one orbit. 

However, shell fluctuations could be reduced if necessary by 
a fast acting feedback path from the shell to the incoming 
radiation, as shown in Fig. 3. 


Conclusion 


The actual construction and geometry of a temperature 
regulated satellite depends, of course, on the particular pur- 
pose and mission of the instruments incorporated in the 
satellite, but the control of its temperature, as outlined above, 
should not be too difficult a problem with any satellite shape. 
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Use of a Catalytic Probe for Detection 
of Dissociated Nonequilibrium States 
in Hypersonic Flow 


STUART A. HOENIG! 


University of California, Berkeley, Calif. 


It is shown that a wheatstone bridge with one catalytic 
and one noncatalytic arm can be used as a sensitive detec- 
tor for dissociation in hypersonic flow. 


NUMBER of authors have discussed the design and 

operation of hypersonic wind tunnels (1 to 3).2 These 
devices may be driven by a shock tube, an electric arc or a 
plasma jet, but they have one common feature. The working 
fluid, generally air, is heated to almost 104 K and allowed to 
expand through a diverging nozzle to produce a high 
Mach number flow. 

The question of attainment of thermochemical equilibrium 
in the downstream test section region of such tunnels has 
been seriously questioned by Bray (4). Experiments by Weg- 
ner and his group (5) have indicated that for a special gas 
mixture (N. + N,.O,), the flow in a supersonic nozzle may 
have an appreciable frozen component. 

It is the purpose of this paper to suggest a new technique 
for detection of frozen dissociated flow in hypersonic wind 
tunnels. 

The instrument consists of a wheatstone bridge, two oppo- 
site arms of which are platinum wires 0.00254 cm in diameter. 
One wire is rendered almost 100 per cent noncatalytic by an 
application of silicon monoxide (6), while the other wire re- 
mains uncovered. 

The use of platinum for wires is convenient for several 
reasons. The wire is available in many diameters down to 
0.00254 cm, from the Sigmund Cohn Corp. of Mount Vernon, 
N. Y. The temperature-resistance coefficient is effectively 
constant up to 1000 K; the wire is easily fabricated by solder- 
ing, and platinum is an excellent catalyst for most reactions. 

When wire of the smallest diameter is used, it is generally 
possible to operate the detector in the free molecule flow 
regime. In this case, no shock wave is formed on the probe, 
and the mass flux may be simply computed in terms of the 
density, flow velocity and projected area of the probe. 

When the two platinum wires are suddenly exposed to a flow 
containing a component, such as atomic oxygen which re- 
combines on the uncoated wire surface, the resistances of the 
two arms will change to a different degree due to the dif- 
ference in the heat inputs. 

We can gain some idea of the sensitivity of the instrument 
by considering the minimum detectable unbalance to be 50 X 
10-* v. Now we draw the usual bridge circuit as shown in 
Fig. 1, where A and B are the platinum wires, and C repre- 
sents a fixed resistor. 

The voltage drop in arms A and B are 


A B 
Their difference is the voltage detected AV. Hence 
A B 
Al E ct [1] 


To obtain a first order result, take C = A and B = A + 6, 
then Equation [1] reduces to 


Received Dec. 22, 1958. 
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Fig. 1 Wheatstone bridge for catalytic detection probe 


An unbalance of AV = 50 X 10-* v for a value of A = 10 
ohms, V; = 10 v gives a difference in resistance between A 
and B of AR ~ 2 X 10~‘ ohms. 

The resistance-temperature coefficient of platinum wire is 
0.003 ohms per deg C, which indicates that a temperature dif- 
ference of AT = 6.66 X 10-? C is the minimum detectable 
value. We can now calculate how much differential heat in- 
put this represents. 

We shall write the first law of thermodynamics in the form 


Q = C,mAT [2] 


where C, and m refer to the platinum wire, nominally chosen 
as l-em long. Using Equation [2], we see that this tempera- 
ture difference represents a differential heat input of 2.5 X< 
10~7cal. The heat of formation of 1 mole of oxygen is 1.18 X 
105 cal, hence, the above heat input represents the formation of 
1.27 X 10!* oxygen molecules. 

All measurements in hypersonic wind tunnels are of a 
transient nature, and we can consider an appropriate measure- 
ment time to be, say, 100 usec. This is about 1 per cent of the 
usual running time, and is a sufficiently short interval so that 
conduction to the wire supports may be ignored and radial 
temperature gradients in the wires may be neglected (7). 

The detectable number of oxygen molecule formations is 
therefore 1.27 X 10'* oxygen molecules/sec, which represents 
a heat input of 2.44 X 10-% cal/sec. 

We shall assume that every oxygen atom that strikes the 
bare platinum wire recombines, while no recombination 
occurs on the coated wire. Present experiments indicate that 
this is a good approximation to the actual situation, and it 
allows a simple estimate of instrument sensitivity. 

Let us consider the use of a device of this type in a more or 
less typical hypersonic wind tunnel, i.e., that operated by 
Perry at Tullahoma, Tenn. (8). 

The operating fluid is air, which is heated by an are dis- 
charge and then allowed to flow through a diverging nozzle. 
Typical stagnation conditions are: T) = 6000 K, Py = 1155 
atm. Under these conditions, the fraction of oxygen disso- 
ciated may be obtained from (9), and is found to be ap- 
proximately 51.4 per cent even after the effect of pressure has 
been accounted for. 

Conditions in the test section at M = 16 are: V = 4.27 x 
10° cm/sec, T = 180 K, P = 1.8 X 10-4 atm, as taken from 
(8). Under these conditions, the mean free path will be ap- 
proximately 2.75 X 10-*cm. Hence, a wire 2.54 K 107% em 
in diameter will satisfy the conditions for free molecule flow, 
(A/a) > 10. In free molecule flow, no shock wave will ap- 
pear on the wire, and we can calculate the number flux hitting 
the wire from N = nV A, where n is the number of particles 
per cm’, V is the flow velocity, and A the projected wire area. 
The value of n ism = 8.2 X 10% molecules/cm*, and therefore 
N is 140 X 10" molecules/sec. Of this flux, 80 per cent may 
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be considered inert, leaving a flux of oxygen molecules 2.80 X 
10° sec-!. If the initial dissociation fraction of ~50 per 
cent is unchanged, as might be expected from the work of 
Bray (4), we would have roughly 1.4 X 108 molecules formed 
per sec on the catalytic wire, or about 100 times the minimum 
detectable fraction, calculated previously. This indicates that 
even a dissociated fraction as low as 1 per cent in the free 
stream might be detected, hence the enormous sensitivity of 
this simple device. 

We have assumed that the heat input from recombination 
can be separated from that due to convection. This involves 
essentially the assumption that the heat inputs are at least of 
similar orders of magnitude. 

The data on heat transfer that Perry presented was for mod- 
els so large that continuum conditions prevailed. In order to 
estimate free molecule heat transfer at M = 16, we shall ex- 
trapolate Perry’s data (8) down to free molecule flow condi- 
tions, and extrapolate the data of Eberly and Wong (10, 11) on 
free molecular flow upward to M = 16 for comparison. 

Perry’s data yields (Nu/V/ Re) = 0.30, which at a Reynolds 
number of 3.33 gives Nu = 0.547. This represents heat trans- 
fer at the stagnation point, and for an average Nusselt number 
we should divide by about 2, giving Nu = 0.273. Eberly 
presents data for free molecular heat transfer to spheres 
which, when extrapolated to M = 16, Re = 3.33, give an 
average Nusselt number of 0.300. The data of Wong for 
cylinders at Mach numbers from 3 to 5 yield an average value 
of Nu = 0.300. Under the circumstances, we are justified in 
taking an average value of Nu = 0.290 for the free molecule 
flow in the hypersonic tunnel. 

The heat transfer coefficient is therefore 6.1 X 107% x 
cal/em* see C, and the heat input from convection may be 
taken from gq = hAAT, where AT is based on the initial dif- 
ference between stagnation and wire temperatures. (The in- 
crease in equilibrium wire temperature over stagnation tem- 
perature can be ignored in this instance since the Knudsen 
number is only about 7.3, and we seek only order-of-magnitude 
results.) 

Under these assumptions, the convective heat input is q = 
0.296 cal/sec; now previously we derived the minimum de- 
tectable heat input to the catalytic wire from recombination 
to be 2.44 X 107% cal/see. Hence, under the worst condi- 
tions the heat input from recombination is 1 per cent of the 
convective heat input, and under the most favorable condi- 
tions it is of the same order of magnitude. 

It is important to note that the use of a wheatstone bridge 
in this application removes almost all the effect of convective 
heating. Under these circumstances, even an atomic oxygen 
content of 1 per cent should be easily detectable. 

We have considered only the oxygen atom content of the 
flow. If the working fluid is air, the products of the discharge 
will include many more unstable species. Since there is very 
little information about the recombination of such species on 
catalytic surfaces, they have been ignored in this analysis. 
No doubt their presence will result in a somewhat increased 
sensitivity for the detector discussed above. 


Conclusion 


A device has been described which can be used for detection 
of nonequilibrium dissociated flow in hypersonic wind 
tunnels. 
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Velocity of Tank Fragments 


SAUNDERS B. KRAMER! 
Lockheed Aircraft Corp., Sunnyvale, Calif. 


An investigation was carried out to determine the posi- 
tions as functions of time of the sections of an exploded 
tank in the upper atmosphere idealized as a 6-ft diameter 
sphere, since the geometry involved afforded a desirable 
distribution of the fragments and permitted immediate 
derivation of descriptive equations. The sphere was pres- 
surized using helium, and was considered to have prima- 
cord over its surface so that upon detonation a spherical 
scattering of fragments would occur. To simplify calcula- 
tions it was supposed that 10 such fragments of equal mass 
and size were generated in a given detonation. 


Analysis for Spherical Case 


XAMINATION of the thermodynamics of this problem 
leads to an equation which relates the rate of change of 
internal energy in the compartment to the rate of energy flow 
through the gap formed by the separating bodies. The ideal 
gas laws are assumed to apply here; thus we have 


a (U) = —wh [1] 


where 


U = total internal energy at any instant 


w = mass flow rate 

h = enthalpy (heat content) per unit mass of the gas in the 
compartment 

Now 

U = mu (2] 
where 

u = the internal energy per unit mass of the gas (helium) 
contained in the compartment 

m = mass of gas in the compartment 


If the gas flow is assumed to be choked (i.e., a critical pressure 
ratio is reached or exceeded—and this is true since the pres- 
sure outside the compartment is practically zero), we have 


PA 


w = 0.2096 —= [3] 
VT 
where 
A = area through which gas flows 
P = value of the total pressure at any time in the compart- 


ment 
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T = value of the temperature in the compartment (in deg R) 
at any time 


Substituting Equations [2 and 3] in [1], we have 
d Ces 0.2096PA h 
— (muy) = 
am 
Since we have assumed the gas to be ideal 
u=C,T 
h=C,T 
where C,, C, are the specific heats (and are constant). 
The gap area through which the originally enclosed gas 
escapes is formed as the spherical surface separates along the 
primacord boundaries. If the original surface is 44R,? and 
some increase in radius takes place at an incremental time 


after detonation, then the area through which the enclosed 
helium escapes is 


R? R,?) = A (escape) [5] 
This may be written as 
4n(R RoR + Ro) =A (escape) [5a] 


If we call the increase in radius X, the space coordinate, 
Equation [5a], may be described by 


)(X A (escape) [5b] 


Equation [4] then becomes 
(mT) = —(0.2096) WT + 2Ry)P [6] 


where 


C 5 
y= for helium 


Equation [6] relates the rate of change of internal energy 
within the sphere to the rate of energy flow through the gap 
formed by the separating bodies. 

The equation relating the acceleration of the fragments 
follows: Suppose we have 10 fragments of equal mass and 
equal area, then each fragment receives 7/5 of the pressure- 
derived force acting. The equation for the acceleration of 


each fragment is 
4rR,? 
p(—)-2 [7] 
10 m, 


where m, is the mass of one fragment. 
From the equation of state for an ideal gas 


PV = mRT (at any time) [8] 


so that, since R is constant, and assuming the volume con- 
tained in the spherical boundary is approximately constant 
during the time the pressure goes to, say, 10 per cent of its 
initial value, we have 


(9) 


where P, T and m are pressure, temperature and mass of the 
gas at any time, and Pp, JT) and mp are their values at ¢ = 0. 
Substituting in Equation [6] we have 

4nX(X + 2Ro) P 


10 
mT 


d 
= —7(0.2096) VT 


If we use Equation [9] and the relationship 


Mo Po 
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Fig. 2 Velocity attained vs. time for fragments of an initially 


pressurized spherical tank 


the result 


is obtained. 
Let 
mo? 1 [13] 
0 
so that 


Substituting Equation [13] into [10], and [9 and 13] into 
[7], the results are, respectively 


d 
dt (m7) = (X + [14] 
and 
X = m? [15] 
3 


If we now solve Equation [15] for m7, substitute the result into 
Equation [14] and differentiate as indicated, we have 


/2 


K (X) = —X(X + 
1 


[16] 
We shall call the lead constant of Equation [16] Co and recall 
that for helium y = 1.66 ~ 5/3, so that we now have 


CoX = —X(X + [17] 


The initial conditions are 


X=X=0 at 


This problem is solved for a chosen value of m, the fragment 
mass, and for various values of Po, the initial pressure in the 
sphere. Changes in Pp give rise to variations in mo, the initial 
weight of gas contained in the sphere. 

To obtain Fig. 1, the sphere has been chosen as 0.02-in. 
thick (the radius is then 3.00167 ft) and of stainless steel; 
hence the value ascertained for m, is 9.225 lb mass. 

Equation [17], because of its very nonlinear nature, could 
not be handled by any standard analytical techniques (other 
than numerical analysis), and was therefore solved using the 
1103A computer. 


Discussion and Conclusions 


Figs. 1 and 2 describe the results obtained. Here, the posi- 
tions attained by the fragments and the velocity of the frag- 
ments as a function of time are shown for the exploded 
spherical tank. These results should be examined with at 
least the following qualifications: 

1 Movement and position have been considered for a field- 
free space, i.e., with no forces acting after final velocity is at- 
tained. 

2 Velocity values noted are higher than would be the 
actual case since boundary separation by primacord would 
leave a gap through which a considerable portion of gas would 
escape, thus lowering the final velocity of the fragments. 

3 Numerical results show that the volume enclosing the 
gas is not constant during the time the pressure goes to 10 
per cent of its initial value; at present, the error thereby in- 
troduced is not known. 

4 Fragment scattering will not be entirely symmetrical, 
since a certain amount of tumbling will occur, yielding then 
a somewhat indeterminate distribution. 
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Burning Rate of an H-Atom Propellant’ 
HOWARD B. PALMER? 


College of Mineral Industries, The Pennsylvania State 
University, University Park, Pa. 


HE RECENT success of Windsor (1)* in trapping pure H 
atoms as a solid prompts added interest in the possibility 
of using H-atom recombination for propulsion. For this 
reason, an attempt to calculate a priori the one-dimensional 
burning rate of an H-atom propellant seems worthwhile. 
The kinetics are about as simple as one could hope to find for 
« propellant, being those for atom recombination, and the 
physical properties of the system can be calculated rather 
well. In treating the flame problem, one has the advantages 
that there is no need for concern regarding an ignition tem- 
»erature—it will be the initial temperature of the gas—and 
that the chemical reaction has zero activation energy. Our 
physical picture is that the propellant consists of H atoms in 
» matrix, which we take to be solid Hy. The upper limit of 
i{-atom concentration in the matrix is taken to be about 10 
per cent, following the theoretical work of Golden (2). 
The propellant is envisioned as being held at 4 K, with 
vaporization at its boiling point at a rate determined by the 
heat flux from the reaction zone. It is assumed that no com- 
ination of H to form H, takes place prior to vaporization, so 
that we need deal only with a gaseous reaction. Although this 
assumption cannot be exact, it is thought not to introduce 
serious error, since for recombination in the liquid phase it is 
necessary for H atoms to diffuse over a distance of a few 
molecular diameters, and the time required for this is long 
(as will be seen) compared to the time necessary for heating to 
the vaporization temperature. A second assumption that 
appears justifiable ex post facto is that the extraction of heat 
from the flame by the vaporizing liquid has no effect upon the 
temperature gradient in the gas adjacent to the surface. If, 
in addition, we assume constant pressure reaction conditions 
and no diffusion, then the gas reaction can be handled as 
though it were a reaction in a static system at constant pres- 
sure, and the temperature can be calculated as a function of 
time. The temperature gradient at the surface will correspond 
to the initial rate of change of gas temperature. 


k 
For the reaction, 2 H + M set H, + M, the rate constant 
k, is about 10'® ec?/mole* sec, almost independent of tem- 
perature (3), so that 


—dCy/dt = [1] 


where the assumption is made that any M (H or H:) has the 
same efficiency in promoting H-atom combination. 
The rate of change of temperature is 


dT /dt = —(dCy/dt) + Co"Cu) 


where AH, is the heat of reaction, and (, is molar heat capac- 
ity at constant pressure. The initial value of d7’/dt, converted 
to a spatial gradient by dividing by the initial gas flow velocity 
So away from the propellant surface, determines the flux of 
heat conducted to the surface for vaporization. Thus the 
heat balance is 


Ao(dT /dt)o/So = SopoAH vap [3] 


where the subscript zero refers to the surface, and AH yap in- 
cludes the heats of fusion and of heating to the boiling point 
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as well as the heat of vaporization of the liquid. For these 
quantities the approximation is made that the propellant has 
the physical properties of pure molecular hydrogen. 

Although in a practical system the pressure might be ex- 
pected to be several atm, an initial calculation at 1 atm seems 
appropriate. For a propellant having one H atom per 10 H: 
molecules 


To = 20.7 K 

Pc = 1.12 X g/cc 

AHvap = 129.3 cal/g 

Ao = 3.8 X 10-5 cal/em deg sec 
Ne = 5.4 X 107 mole/ce 

CoH? = 5.35 X 10-4 mole/cc 

AH, = —103.3 kcal/mole He 

= 5.9 X 10-4 mole/cc 


)o = 4.96 cal/mole deg 

(C,"2)) = 6.62 cal/mole deg 
From Equation [1], (—dCu/dt)) = 1.72 X 104 mole/ce sec. 
Equation [2] then yields (d7'/dt)) = 2.33 X 10"! deg/sec, and 
solving [3] for So, we get 


S = 7840 cm/sec 


ll 


The rate of change of temperature seems enormous, but we 
note the temperature gradient at the surface is (dT/dr)) = 
2.98 < 107 deg/em. At 20 K, the mean free path in a gas 
is of the order of 10~* cm, so that the temperature gradient is 
roughly 30 deg per mean free path. This may imply that a 
normal thermal conductivity coefficient should not be used; 
however, it may be pointed out that the experimental tem- 
perature profiles in typical gaseous flames (4) show gradients 
of the order of degrees per mean free path, and it is generally 
assumed that normal thermal conductivity coefficients can be 
used in theoretical treatments of such flame systems, aside 
from the anomalies brought about by entropy changes (5). 
A further point is that the temperature gradient is of the order 
of that found in weak shock waves, to which theories of shock 
structure using normal transport coefficients have been applied 
with fair success (6). 

The linear burning rate of the solid will be 


S, = Sopo/p, = 125 cm/sec 


where the solid has been assumed to have the density of liquid 
H.. The mass consumption rate per unit area of the solid is 


m/A = Sp. = 8.76 g/cm? sec 


The effect of pressure upon the burning rate can be esti- 
mated if the effect of pressure upon the AHvap term can be 
approximated. Because the variation of the enthalpy of liquid 
H, with pressure is extremely difficult to estimate, we may as 
a first approximation take AHyap to be a constant. Raising 
the pressure to 2 atm raises the vaporization temperature by 
2.3 K, to 23.0 K. The new Sp is 1.04 & 104 cm/sec, and the 
new S, = 3.02 X 10? cm/sec. For an empirical fit of the type 
S, = ap”, this gives n = 1.28. 

Some of the initial assumptions can be tested on the basis of 
the results. Diffusion coefficients in the liquid phase have 
values around 10-5 cm?/see at room temperature; if we sup- 
pose the activation energy for diffusion to be 5 keal, D for 
liquid Hy, at 20.7 K is of the order of 10~*5 em?/sec. Although 
this figure may be unrealistically small, it would nevertheless 
appear to guarantee that the effect of diffusion in the liquid 
will be nil. In the particular propellant chosen for discussion 
here, there is a mechanism other than diffusion available for 
atom recombination in the liquid phase. If the matrix had 
been Ne, for example, diffusion would have offered the only 
mechanism; however, with a mixture of H and He, recombina- 
tion of H could conceivably occur in the liquid through a sort 
of chain transfer involving formation of H; and its redissocia- 
tion to H, and H. This is just the ortho-para hydrogen con- 
version reaction, which has an activation energy around 7.5 
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keal (7) and which therefore will be negligibly slow at 20 K. 

The rate of heat extraction from the flame by the vaporiza- 
tion of the liquid may be compared to the rate of heat re- 
lease. Per square cm of surface, the vaporization rate re- 
quires about 10° cal/sec. The material vaporized from this 
area has an initial heat release rate of about 10*cal/sec. Thus 
the gaseous reaction is virtually unaffected by the vaporiz- 
ation requirement. 

In conclusion, it is noted that the expression for the linear 
burning rate 


S;? = Ac AH, /[2 ps*AH vap(Cp)o*”] 


where M,*’ is the mean molecular weight (initially), is similar 
to simple gaseous flame theories such as that of Damkéhler 
(8) except that temperatures do not appear in it. Since the 
ignition temperature was specified as the initial gas tem- 
perature in this case, the term expressing the enthalpy change 
for heating of reactants to the ignition temperature has be- 
come just the AHysp. The estimate of the temperature 
gradient at the ignition point in the Damkohler treatment, 
which involves the difference between initial and final tem- 
perature, a correction factor and an average reaction rate, 
has been replaced by a more direct calculation using known 
kinetics, concentrations and thermal quantities.‘ 
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* The adiabatic flame temperature for the 1:10 H-H: propellant 
at 1 atm is only 720 K, almost independent of pressure for pres- 
sures of latmor more. The ideal /,, for this system, operating in 
a vacuum (operation in the atmosphere would, of course, be ex- 
ceedingly inefficient for pp = 1 atm) is 472 sec. For a 1:5 
H-H;, propellant, if it could be made, the 1 atm adiabatic flame 
temperature is 1330 K, and the corresponding /.p is 660 sec. 


Generalized Solution for the Problem of 
Pipe Flow With Friction and Heat 
Addition 


R. W. McJONES! and M. L. CAPLAN? 
Marquardt Aircraft Co., Van Nuys, Calif. 


A generalized chart method is presented for the rapid 
solution of pipe flow problems involving compressible 
fluids, including the combined effects of friction and heat 
addition. A typical chart is shown for the case of subsonic 
air flow through a duct of unspecified but constant wall 
temperature. Similar charts can be constructed for other 
problems. 
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HE CONSTRUCTION of a generalized chart occurred 

to the authors as an extension of Shapiro’s graphical 
presentation of particular calculations on page 246 of (1).* 
It appeared that if enough particular cases were calculated 
from low Mach number to choking and if the results were 
properly plotted, the resulting chart would cover all possible 
intermediate problems. Just as the calculations were com- 
plete and various plotting parameters were being considered, 
previous unpublished work by Frank Rom of the NACA 
Lewis Laboratory was discovered. He had done substantially 
identical work in 1953 and had prepared a completely work- 
able chart. The present calculations are in excellent agree- 
ment with Rom’s, and the present chart is little more than « 
refined version of his. 


Development of Generalized Solution 


Compressible fluid flow with friction and heat transfer in- 
volves interrelated changes in Mach number, total pressur¢ 
and temperature. These changes are clearly described by th« 
influence coefficients of (1) which lead to the following equa- 
tions for flow in constant area ducts with constant mass flow, 
and no work done by the fluid 


AM? _ [1 + + — 


(1 — M?) To 
+ [(y — 1)/21M%} AL 
(1 — M?) af D 
AP, — yM? AT, AL 


Any generalized solutions are clearly dependent upon the 
existence of a fixed relationship between wall friction and gas 
temperature change. The possibility that a suitable relation- 
ship exists is suggested by Reynold’s analogy, and (2) de- 
velops the required equations for both the case of constant 
wall temperature (considered here) and the case of uniform 
heat flux. With some notation changes, Equation [22] of (2) 
becomes 


AL 
D 0.0092 (kw T/T )?(1 — T/T») 


[3] 


Notation here is consistent with reference (2), except that 
the present R, is defined as pV D/u., where p and V are based 
on average values from continuity, but yu. is gas viscosity at 
the wall temperature. The various test reports mentioned 
in (2) have since been consolidated and published as (8). 

Process lines for the present chart were calculated on a digital 
computer using an available program based on Equation [1]. 
A wall temperature of 2000 R was used for the bulk of the 
calculations; however, check calculations using 7,, = 900 R 
and 3460 R verified the validity of the generalized chart. 


Selection of Chart Format 

The various quantities which may be of interest in flow 
problems include total and static pressure, Mach number, 
velocity, total and static temperature, density, etc., and are too 
numerous to include on a single chart. When the channel 
area, wall temperature and mass flow are constant, the con- 
tinuity equation may be reduced to the following form which 
suggests the format finally selected 


VTi/T 


= constant [4] 
P,A*/A 


The gas-to-wall temperature is placed along the abscissa 
together with the corrected tube length parameter; this scale 
is made linear in tube length to give a better physical repre- 
sentation and to facilitate solution of problems where tube 
length is given and temperature ratios must be found by trial 
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Fig. 1 
104 (using viscosity at 7'..). 


Generalized chart for air flow in a pipe with constant wall temperature. 
At inlet, 7’,; = 1200 R, M; = 0.34, P:; = 200 psia. 


Example: Geometric L/D = 63, 7. = 2000 R, R. = 
Calculate L/D X 10/R2-? = 63 X 10/6.3 = 100. Locate 


inlet point at (7,/7'~)i = 0.6 and M; = 0.34; read inlet Pr; = 2.65, and inlet corrected tube length parameter = 110. Follow process 
line to exit corrected tube length = 110 + 100 = 210. Read at exit, M. = 0.61,(7,/Tw). = 0.85. Pr, = 2.03. Calculate exit gas 
temperature = 0.85 < 2000 = 1700 R, exit total pressure 200(2.03/2.65) = 153 psia. 


and error. Relative total pressure is placed along the ordi- 
nate, which is made logarithmic for convenience. Lines of 
constant Mach number are then constructed from Equation 
[4]. Although y = 1.4 is implied in the A*/A vs. Mach num- 
ber relationship, check calculations show negligible overall 
chart error for air temperature approaching 7, = 3460 R. 

The chart thus defined, Fig. 1, is a plane of constant mass 
flow, so the typical process lines computed as described above 
(as well as any other lines which might be drawn) auto- 
matically satisfy continuity. 

The Fanno line plotted along the right-hand edge of Fig. 1 
relates total pressure, Mach number and corrected tube length 
for the case of zero heat transfer. It is also useful in problems 
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where the final gas-to-wall temperature exceeds the chart 
cutoff at T,/7.. = 0.99. Asa matter of interest, at T,/T.. = 
0.5, the temperature term in Equation [2] is only 3 per cent of 
the friction term, and at 7,/T, = 0.99, it is down to 1 per 
cent. 


Chart Construction for Other Cases 


Construction of a chart for the case of heat transfer from the 
gas to an isothermal wall follows exactly the procedure out- 
lined herein, except that the total pressure change due to heat 
transfer will be of opposite sign. The case of uniform heat 
flux is of frequent interest and can be handled by this method 
if a generalized relationship between gas temperature and tube 
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length can be obtained. The results given in Fig. 6 of (2) sug- 
gest that a single curve of relative gas temperature vs. cor- 
rected tube length can be constructed such that it will ap- 
proximate flow with constant heat flux between any two arbi- 
trarily selected end points. Once this curve is obtained, the 
generalized chart can be readily constructed. A problem in- 
volving any arbitrary wall temperature variation can, of 
course, be handled with the present chart by assuming an 
equivalent tube having step changes in wall temperature to 
approximate the given temperature profile. The overall 
solution is then obtained by combining the separate solutions 
for each constant temperature tube increment. 
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On Satellite Librations 


G. M. SCHINDLER! 
University of California, Santa Barbara, Calif. 


When a satellite is not perfectly round but oblong, it will 
oscillate about an attitude of equilibrium somewhat like a 
pendulum. These oscillations are called ‘‘librations.”’ 
The present paper deals with the librations of a satellite 
built in the shape of a dumbbell, which consists of two 
spheres M, and M; of equal mass m separated by a rigid 
but weightless rod of length 2e. The center of gravity of 
the satellite is supposed to move exactly in a circular orbit. 
It is shown that the frequency of such oscillations depends 
very strongly on the angle of excursion and that it ranges 
from a finite positive value to infinity when the excursion 
varies between +77/2 and zero. 


HE RESULTS to be found in literature on the problem 

of determining the librations of a satellite built in the 
shape of a dumbbell differ from each other as far as already 
rather strong simplifications on the basic equations have been 
made. In connection with such simplifications, it is not sur- 
prising that a frequency of the dumbbell was found to be of 
the magnitude of the satellite frequency around the Earth. 
This is true only for a proper choice of the excursion angle, 
and it shall be shown here that by a suitable choice of the 
excursion, it is always possible to perform a libration which 
period is exactly the same as the period of the satellite in the 
circular orbit. 


Discussion 


Let p be the radius between the center of gravity of the 
dumbbell system and the center of the Earth (Fig. 1). At 
time ¢t the rod connecting the two spheres M, and M, of equal 
mass m builds the angle ¢ with the vertical. Now, due to 
their different distances from the geocenter, the two spheres 
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0 
Fig. 1. Geometry of dumbell-shaped satellite around the Earth 


suffer different attractions, viz. 


_ 
(p? + e? — 2 pe cos ¢) 
mgk? mgR? 


ri? (p? + e? + 2 pe cos ¢) 


where 


R = radius of the Earth surface 
g = gravitational acceleration at Earth surface 


Similarly, because they are forced to participate in the same 
orbiting revolution at the common angular velocity Q around 
0, the two masses experience different centrifugal forces 


= mr,Q? = mQ*%Xp? + e? — 2 pe cos ¢)'? 
Ce = = + e? + 2 pe cos 


If the angle ¢ is taken to be positive in the counterclockwise 
direction, then the restoring force is acting on the first sphere 


RR? 
mee = —(G;’ — = sin (= [3] 
1 


while the second sphere suffers the restoring force 
R2 
meg = G,’ — = me sin (= rs) [4] 


The total restoring moment is then given by the sum of both 
forces in Equations [8 and 4] 


1 1 
2meg = mgpR? sin (2 *) [5a] 


or 


l 1 
(p? +e 5b 
+e? + 2 pe cos (p? + e? — 2 pe cos 


From Equation [5a] it follows that the total restoring force 
2me¢ is pointing toward the positive direction of ¢ if (34/2) < 
g < 2m and (2/2) < ¢ < =, and that it points toward the 
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negative direction of if 0 < < (4/2) and < ¢ < (3 77/2). 
This causes an oscillation around the vertical axis ¢ = 0. which 
is the attitude of stable equilibrium. The position g = 7/2, 
however, can only be considered to be the attitude of non- 
stable equilibrium, since the dumbbell always tends to the 
stable position g = 0, if it is deflected from the position g = 

Since in the general equation of motion [5b] the angular 
velocity Q does not occur, there is no relation between the 
periods of the dumbbell and the satellite around the Earth. 

Equation [5b] is valid for arbitrary orbit. We now assume 
the satellite to be in circular orbit, that is, p(t) = constant. 
‘Then we can integrate Equation [5b] with respect to the time 
t. and we obtain 


h [ 1 1 
0(“) L(p? + e? — 2 pecos ¢)'/? (p? + e? + 2 pecosg)'/? 


where K is an integration constant, which is determined by 
the initial condition of t = to 
= Go Alto) = Go [7] 


Kk then turns out to be 


It is hopeless to try to solve this hyperelliptical integral. For 
actual purposes, instead of solving Equation [13] one can 
solve the differential equation for t(¢) 


dt 4e 
re 
with the initial condition 
Ugo) = [14b] 


which may be done by computer application. The time T is 
then given by 


[6] 


T = (0) — t(go) = (4 — k) [14e] 


However, the integral in Equation [13] can be investigated for 


Yo — (2/2) and for g — 0. 


For this purpose it is convenient to use the substitution 


—A cos2X [15] 


cos g = 


[8] 


k= 1 1 
R} 9  (p? + e? — 2 pe cos (p? + e? + 2 pe cos 


Now, the time for one complete cycle swing of the dumbbell 
is given by 


[9] 
where 
1 4 1 
Vp? + e? — 2 pecosy + + 2 pecosy 


= 


where A is given by 


(e+ [16] 
ep 


= - 
re 


[10] 


We then have 


Z(¢) is limited between the boundaries 


2 2 


8eA e 


V p? + e? —2pecosy =n 4/2(p? + e?) cos X 
V p? + + 2 pe == +/2(p? + sin 


T = 


This implies that for sufficient high initial velocity go, the 
constant K can be arbitrarily large, so that the time 7 in 


and 
sin2AdX 
[18] 
Vi — A? cos? 2 [(1/cos A) + (1/sin A) — (1/eos Ap) — (1/sin 
where 
vos Go = —A cos2 Ao 
1 = —Acos2A, 


Equation [9] can be arbitrarily small. However, we are not 
very much interested in this obvious result, since the most 
interesting case is that where we release the dumbbell from 


1 — A? cos?2X 


the position ¢ = ¢ for free oscillation. This would be the 
case where g = 0 and 0 < gv < (2/2). We then have 


[19] 


We further use the expansions 


= (1 — A cos2A)(1 + A cos2X) = 


(Ar + P(A = r)] 


1 1 
+ e? — 2 pe cos (p? + e? + 2 pe cos 


For reasons of symmetry we obtain 


dy 


go 
p? +e— 2 pe cos ¢) + (1/V p? +e? + 2 pecosy) + K 
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ak 
[20] 
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[17] 
] 
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Where P(A; — A) and Q(A — Apo) are regular power series in \ 
for Ay < A < A; which start at least with the first power of 
(A; — A) and (A — Ao), respectively, and 

2A sin 2 

sin Ay cos Ao 

cos? Ay sin? Ao 


We then have to solve the integral 
sin2AdX 
— + PO — — Ao) fdr + QA — Ao)] 


are the same. 
In order to prove this theorem, we first have 
[28] 
Pi 
As can be seen in Equation [16], the constant A does not 
change its value. Since the initial conditions are supposed to 
be the same, the excursion ¢ is the same. This implies that 
the boundaries » and \, are the same in both cases, as can be 
recognized in Equation [19]. Then the integral in [18] docs 
not change at all. If we denote by 7; and 7, the period time 
in both cases, respectively, we obtain from Equation [18] 


V eid: Vim — AMA — Ao) [ (| 


Now, the integral 


WA 


exists, where B(p, qg) is the Eulerian Beta-function. Since P 
and Q are regular power series in the entire integration in- 
terval, the integral J in Equation [22] exists for all values of 
Xo and ), where both c; and d, are unequal to zero. As can 
be recognized from Equation [19], the range of Xp and A, is 
limited by 


Tv us 
24 
p< [24] 


if g varies from zero to (17/2). This implies that the mag- 
nitudes d, and c; are unequal to zero, except in the case when 
Xo tends to 7/4, that is, for g — (7/2). In this limit case, 
d, vanishes, and, therefore, the integral in Equation [22] 
is infinite. This proves the fact that the time 7 tends 
to infinity if the excursion ¢ of the dumbbell tends to 
a/2. Now, A; is greater than 7/4. If Ao tends to di, the 
integral J in Equation [22] tends to (1/Wcid,)e sin 2d. 
The limit period 7p is then obtained from [18 and 19] 

lim T= 1 [25] 
Using here the expansion 


p—e pte p 


the limit period 7» is found in a first approximation (all higher 
terms of the ratio e/p are neglected) to be 


where 7 is the period of the satellite in the circular orbit. 
This result has already been found by Klemperer and Baker 
for small excursions. Since the series in Equation [26] 
consists of positive terms only, the exact limit period 7» is 
even smaller than that given by the approximation in Equa- 
tion [26a] 


1 
T) < <06T 7 
0 V3 [27] 


From the fact that the dumbbell period ranges from a finite 
value, which is even less than 0.6 7, to infinity, it follows that 
it is always possible by a proper choice of the excursion to 
perform a period which equals exactly the satellite period. 

In conclusion, we prove the following Theorem: If for two 
different dumbbell systems py, and ps, €2, where p; and are 
the radii of the circular orbits, respectively, and e; and é: are 
the radii of the dumbbells, respectively, the ratios e,/p, and 
€2/p2 are the same, then the square of the ratio between both 
dumbbell periods is equal to the cube of the ratio between the 
radii of the dumbbells, supposing that the initial conditions 
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Using here Equation [28] we have 


which proves the theorem. This theorem is an analogue tv 


K pler ’s third law 
T2 P2 


where 7; and 72 are the times of one period of the satellites. 
respectively. From Equations [30 and 31] one can conclud 
that there exists a linear relation between 7 and r 

T=C-r [32) 


However, the proportional factor C depends essentially on the 
ratio e/p and on the initial condition go, as shown previously. 
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Time for Escape Using a Succession of 
Small Impulses 


ROBERT E. ROBERSON! 


Autonetics Div. of North American Aviation, Inc., Downey, 
Calif. 


J B. PETERSON? has pointed out that an efficient 
© escape from an orbit can be made with low thrust if the 
thrust is applied in a series of short impulses at perigee. The 
total impulse required is the same as if it were delivered in one 
big burst. He observes that the penalty one pays is in time. 

It is of some interest to see just how bad this time penalty 
is. Consider an initial circular orbit of radius R about a body 
whose mass times the gravitational constant equals K. The 
velocity increment required for escape is AV = V2K/R — 
VK/R. Suppose, however, that thrust limitations prevent 
adding more than Av = AV/N in any (essentially impulsive) 
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ESCAPE TIME/INITIAL ORBITAL PERIOD 
AS A FUNCTION OF THE NUMBER OF THRUST IMPULSES 


7 ~~ AN UPPER BOUND ON ASYMPTOTIC BEHAVIOR 


100 1000 
To 
Fig. 1 Time for escape in units of initial orbital period 
(abscissa)as a function of number of impulses required to pro- 
duce escape velocity (ordinate) 


burst. After each revolution, at perigee, Av is added. The 
problem is to find the time required to escape by this process. 
Let 
T) = period of circular orbit 
T; = period of the elliptical orbit during the ith revolution 
(that is, after application of the 7th small impulse) 
the perigee velocity at the beginning of the 7th revolution 
after the impulse 


It is easy to show that 7,;/T, = (2 — Rv;?/K)~*2.. Now 
», = VK/R + ido (i = 1,2,...,N). Hence 


a 72 |-3/2 
— 2-l)=- — 1)? — 
T, | (Vv ) (v ) 


It should be apparent that escape will occur after V impulses, 
whence the time to escape is 


(The velocity after N impulses is vy = VK/R + Ndv = 
VK/N + AV = V2K/N, precisely escape velocity.) 

The quantity 7'y/7> is shown as a function of N in Fig. 1 for 
NV =1,2,..., 100. For larger an approximate asymptotic 
result is useful. The value of the sum for 7'y/T> is approxi- 
mated by the integral 


1-1/2N 
— 2/2 Dax — (V2 — 


although this integral will give a value a little too high be- 
cause of the rapidly increasing slope of the integrand as x 
nears unity. Evaluating the integral and keeping only the 
most rapidly growing term in N, one obtains 


4. ,% 1 
log ~ 5 log N + 5 (1+ 
This relation is shown as the straight line of the figure, from 
which it appears that the approximation indeed gives a value 
of T'y/To which is slightly too high. 

This result is instructive. It strongly bears out one’s intui- 
tive feeling that the time delay for an escape of this kind is 
prohibitive. As an example, consider the case in which a 
satellite is in a 800 km circular orbit with velocity 7452 m per 
sec, escape velocity being 10,539 m per sec. Suppose that 
only 1 m per sec can be added in any one revolution, whence 
N = 3,087 and Ty/T) ~ 378,500. Since 7) = 6052 sec, 
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Ty =~ 26,520 days—hardly an acceptable break-away time. 
Even if the velocity increment is given in 30 bursts of 102.9 
m per sec each, the total break-away time is about 21 days, 
The latter is an artificial case of course. The only reason for 
going to this kind of escape maneuver in the first place is in 
connection with the use of very low thrust systems in which the 
velocity increment in any one burst is very much smaller than 
the required 3087 m per sec. 


Attitude Stabilization by Means of 
Solar Radiation Pressure 


ROBERT L. SOHN! 


Space Technology Laboratories, Inc., Los Angeles, Calif. 


This paper demonstrates the feasibility of stabilizing the 
attitude of a space vehicle by means of solar radiation pres- 
sure. If a suitable weathervane-type tail surface is at- 
tached to the vehicle, solar radiation pressure, acting on 
the weathervane, will stabilize the vehicle against dis- 
turbing torques which arise from gravitational eccentrici- 
ties (dumbbell effect), meteorite impacts and solar radia- 
tion pressures. 


NE VERY difficult problem confronting the designer 
of long duration space flight vehicles is that of con- 
trolling and stabilizing the orientation of the vehicle for long 
periods of time. For example, it might be desirable to orient 
a communications satellite in a given direction, and maintain 
that direction to close tolerances. Or during lunar or inter- 
planetary flights, a constant attitude might be required for 
various reasons. Now, if disturbing torques, such as solar 
radiation pressures or gravitational torques cause the vehicle 
to rotate away from the desired orientation, a suitable corree- 
tive torque must be applied. At present, a fluid jet of some 
type is usually assumed to supply the restoring moment.? 
Even though the disturbing torques are small, jets cannot be 
used over long periods of time without consuming large quan- 
tities of propellant. 
Solar Radiation Attitude Stabilization System 

Several people have proposed the utilization of solar radia- 
tion pressure as a thrusting device for space vehicles, ‘‘solar 
sailing”’ as it has been called (1). They have demonstrated 
that although solar radiation pressure is low, it nonetheless can 
be significant for thrusting purposes. 

It seems entirely reasonable that solar radiation pressure 
could also be used to good advantage in an attitude stabiliza- 
tion system. The disturbing torques (solar radiation pres- 
sure, gravitational torques and meteorite impact) are very 
small, and, as will be shown, are of the same order of magni- 
tude as the solar radiation pressure restoring torque. These 
effects will be examined in turn to compare the disturbing 
torques with the solar weathervane restoring torques. 


Solar Radiation Pressure Disturbing Torque 


Consider a vehicle moving broadside to the sun in the ab- 
sence of any gravitational field. A disturbing torque will exist 
if the center of pressure is not coincident with the center of 
mass. It is possible to stabilize the vehicle against this torque 
by attaching a weathervane tail surface to the vehicle in such 
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Earth 


Assume 


Mass = 31 slugs 
Weight (on Earth) = 1000 lb 
x 


= 1 ft 
b = 5ft 
S, = 50 sq ft 
Ss = 50 sq ft 
(7 x 7 ft) j 
length = 15 ft 


Fig. 1 Use of a solar stabilizer sail during a lunar flight 


a way that solar radiation strikes the weathervane and rights 
the vehicle (see Fig. 1). 
The disturbing torque of the vehicle 7, is 


T, = 
where 
P = solar radiation pressure, lb/ft? 
S, = exposed surface area of vehicle, ft? 
x = distance from center of mass to center of pressure, ft 


The righting torque of the sail 7, is 


T, = PS,by (for small W) 
where 
S, = surface area of stabilizer, ft? 
b = stabilizer moment arm, ft 
yw = yaw angle, radian 


Equating the torques givesthe equilibrium yaw angle y as 

bS, 

11.5 deg 


where the numerical value is obtained by inserting the design 
numbers of Fig. 1 into the above equation. If the surface is 
movable, the yaw angle could be held to zero deg by deflecting 
the surface through an angle 6 


11.5 deg 


Gravitational Disturbing Torques 


The upsetting gravitational torque is 


T, = Iw? (1 +2 £) sin B cos B (general) 
ra? 
= 3/w? sin B cos B (satellite) 
where 
w = rate of angular movement about attracting mass, radian 
per sec 
r = radius out from attracting mass, ft 
I = moment of inertia, slug-ft? 
8 = flight path angle measured from local gravitational vec- 


tor, radian 


In order to compare the solar radiation pressure stabilizer 
capacity with the gravitational torque disturbance, two cases 
were studied: Vehicle in a 300-nautical mile altitude satellite 
orbit, and the vehicle in an escape shot to the moon. 

In the satellite case, it was found that the vehicle design 
shown previously (with a stabilizer area of 50 ft?) was capable 
of trimming out the gravitational torque under the worst 
conditions, namely when 6 = 45 deg. 

In the lunar ease, a slightly different vehicle was assumed. 
The stabilizer design was held fixed, but the moment of 
inertia was increased by a factor of 10, to J = 300 slug-ft?. 
The sail area required to hold the vehicle in equilibrium along 
the flight path as a function of time after launch is shown in 
Fig. 2. Since the flight to the moon takes approximately 2.5 
days (3600 min), the graph represents a small portion of the 
total flight. After 27 min of flight, the 50-ft? stabilizer can 
trim out all gravitational eccentricities, even with the large 
moment of inertia. 
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Fig. 2 Comparison of required and available stabilizer capacity 
during lunar flight 


By a slight rotation of the vehicle, so as to align its center- 
line with the local vertical, all of the gravitational torque 
could be eliminated. Since there is no compelling reason to 
orient the centerline along the line of flight, at least for the 
first few hours of flight, there is no reason why the solar 
stabilizer could not be used during escape shots. 

If it were mandatory to hold some arbitrary attitude during 
escape, jet boosters could be used for the first 30 min, until 
the solar stabilizer became adequate. These jets would be 
tiny in size and could be attached to the sail framework. In 
this way, the basic sail manipulators would also orient the 
jets. Further, the same servomechanisms used to control the 
sail could simultaneously control the jets. 


Meteorite Disturbances 


Finally, it is possible for meteorites to disturb vehicle 
orientation. These disturbances are extremely small. Fig. 3 
has been prepared to show their effect. The distribution fre- 
quency of the meteorites is based on data from Carter (2). 
Fortunately, the large meteorites are encountered so seldom 
that the vehicle has no difficulty in righting itself before the 
next meteorite strikes it. The smaller meteorites have so 
little impulse that, although they are encountered more fre- 
quently, they cause little disturbance. The meteorites were 
assumed to strike 10 ft from the missile’s center of mass. 
Meteorite disturbances should pose no problem. 


Dynamic Stability 

With the weathervane fixed at some equilibrium angle, the 
vehicle is statically stable, but dynamically neutrally stable. 
If disturbed, the vehicle will oscillate indefinitely with an 
amplitude equal to twice the disturbance angle. This can 
be demonstrated as follows. The excess torque is 


=7,-—17, 


The solution for y can be written as follows, using the subscript 
0 to indicate initial conditions. Time ¢ is in see 


y= tin VAtt cos VA t+ 
PS,b 
A= 7 (square of the natural frequency) 


A rate sensor can be placed in the control system to deflect 
the surface in proportion to the yaw rate (6 = ky, k being the 
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Length of Time Before Encounter (~hrs) 
(99% Certainty) 


Fig. 3 Time to restore vehicle after meteorite disturbance 


Restore Time (~ sec) 


fe) 
5 
= re) 
ov 5 10 15 20 
g Time (minutes) 
-5 


Fig. 4 Time history of yaw following a 5 deg disturbance (with 
rate damper) 


proportionality constant). The vehicle is then damped and 
will return to equilibruim. The equation of motion becomes 
(ignoring the mass of the stabilizer) 


= + 6) — PS.x 


The solution for the above equation is 


(r:/rn) — 1 (r2/n) — 
(Yo/r) en 
(r2/1)) — 1 
Ak | 
2 2 Ak ; 
| 
A typical time history of the motion is shown in Fig. 4 for 
the vehicle described in Fig. 1. It can be seen that the 


vehicle corrects a 5 deg disturbance to within 1 deg in about 
12 min. 
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Conclusions 


The brief calculations presented above serve to indicate 
that a solar weathervane stabilizer is capable of orienting a 
vehicle traveling in solar space. Although the corrective 
torques produced by the stabilizer are small, they are none- 
theless compatible with the external disturbing torques. 
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Technical Comments 


Further Comments on the Powered 
Flight Trajectory of a Satellite 


THEODORE N. EDELBAUM! 
United Aircraft Corp., East Hartford, Conn. 


NUMBER of authors, e.g., (1, 2),? have treated satellite 
launching mechanics by assuming that if 

(a) The distance covered during powered flight is small 
compared to the Earth’s radius, then 

(b) The gravitational field may be taken as uniform during 
powered flight, or, equivalently, the centrifugal acceleration 
due to motion around the Earth may be neglected. 

The assumption is usually made that (b) follows from (a). 
That this is not necessarily true for speeds on the order of 
satellite velocity has been demonstrated analytically in (8). 
In a recent note (4) the author of reference (1) has stated that 
reference (3) is in error and that assumption (b) does follow 
from (a). The author uses two arguments to support this 
thesis. The first argument is essentially a restatement, in 
vector notation, that it is ‘‘clear”’ that (b) follows from (a). 
The second argument uses a higher approximation that allows 
for first order variations in the gravitational field. The result 
of this approximation reduces to the result of approximation 
(b) when, and only when 

(c) The powered flight time is small compared to V R,/2g 
= 95min. (R, = radius of the Earth.) 

If (b) does really follow from (a), then the result of the 
higher approximation should reduce to the first approximation 
when (a) is satisfied and not when some new requirement (c) is 
introduced. 

In fact, (c) is very restrictive. This may be shown by using 
some results for the optimum duration of powered flight. In 
(2) it is shown that if (b) is satisfied, the optimum duration of 
powered flight is 


where 


z: duration of powered flight 
y = altitude 


For this time to satisfy (c), the square root of the distance 
covered during powered flight must be small compared to the 
square root of the Earth’s radius, a much more restrictive as- 
sumption than (a). 

A more important result is the optimum duration of powered 
flight for a realistic inverse square, central gravitational field 
(5). For this case the duration of powered flight is equal to 
approximately half of a satellite’s period, the 1/2 x times the 
time given in (c)! 

In closing, it may be said that (b) is true if both (a) and (ce) 
are satisfied, but that this combination is a very restrictive 
condition which is usually not even approximately satisfied by 
real satellites. 
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Comment on “Frictional Electricity in 
Missile Systems” 


WILLIAM J. GUMAN! 


Rensselaer Polytechnic Institute, Troy, N. Y. 


N READING Mr. Molmud’s article® “Frictional Elec- 
tricity in Missile Systems,”’ I recalled an unusual occur- 
rence of electrification. While performing a test at a Macl: 
number of 2.5 in a continuous flow supersonic wind tunnel in 
which the air was highly contaminated with oil mist, it was 
noticed that a buildup of an electrical charge on the sting 
support discharged to the sting arm, burning through one o! 
the two plastic pressure leads imbedded in the arm. The 
discharge occurred at a sharp corner of the sting support and 
left a pattern of concentric alternately light and dark circles 
on the arm (apparently the residue of burned plastic). Figs. 
1 and 2 show the side view and top view of this pattern (see 
arrows), as well as the burned pressure lead. The discharge 
was subsequently eliminated by rounding off the sharp 
corner. 


FLOW 


DIRECTION 


Fig. 1 Side view 


Fig. 2 Top view 
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Comments on “Photochemistry of the 
Upper Atmosphere as a Source of 
Propulsive Power” 


STERGE T. DEMETRIADES'! 
Aerojet-General Corp., Azusa, Calif. 


fPXHE AUTHOR notes with interest the excellent review 

paper by A. F. Charwat, ‘Photochemistry of the Upper 
Atmosphere as a Source of Propulsive Power’ which appeared 
in ARS Journat (1).? For the record, it is appropriate to 
draw attention to some previous work of this author and his 
co-workers (2, 3, 4) carried out independently at the Astro- 
nautics Laboratory of the Aerojet-General Corp. under 
sponsorship of the Office of Advanced Studies, Air Force 
Office of Scientific Research, contract no. 49(638)-111, and 
summarized in a recent note (5). Charwat’s independent 
conclusions are in agreement with the previous conclusions of 
this author. 

It is unfortunate that the article (6) which described the 
so-called “First Mono-Atomic Oxygen Ramjet” did not 
clearly differentiate between the work done at Aerojet and 
the work done elsewhere. However, this author’s conclusions 
were amply documented in reports available to “qualified 
requestors” (3,4). Therefore the publication of his negative 
results in the open literature, at the expense of more fruitful 
projects, was not considered justified. The work of this 
author indicates that the defects of the fuel-less atomic-oxygen 
powerplant can be summed up by the following statement: 
At orbital speeds, its drag is too high; at suborbital speeds, the 
lift is too low. Although the thrust of a ‘fuelled’ atomic- 
oxygen orbital powerplant (i.e., a powerplant which uses the 
energy of recombination of oxygen atoms to heat and expand 
an expellant or propulsive fluid) can be made greater than its 
drag, there is no real advantage to be gained by this scheme 
for extending the lifetime of low altitude satellites. With the 
same fuel weight one could achieve better results with an or- 
dinary microrocket (firing in pulses if necessary). 

The usefulness of the atomic oxygen in the upper atmos- 
phere for auxiliary power generation is doubtful, since at or- 
bital speeds the stagnation energy of 1 cc of air is at least 10 
times greater than the recombination energy available in the 
same volume. At the same time, at altitudes above 110 km 
solar energy can provide, ideally, more power per unit area 
of receiver than atomic oxygen. Therefore, it appears that 
for practical power-generation purposes, the atomic oxygen of 
the Earth’s atmosphere can be safely neglected. However, 
the measurement of accommodation coefficients, collision 
cross-sections, etc., of oxygen and other atoms make the study 
of the chemosphere and of artificial flows containing large 
fractions of oxygen and nitrogen atoms, a problem of con- 
siderable scientific interest. 

This author feels that it is now time for the study of powered 
flight in the upper atmosphere to focus on an entirely different 
objective. Although it has been recognized for some time that 
great economies can be effected in space travel by splitting 
the problem into two distinct phases, the booster phase (con- 
cerned mainly with escape from the Earth’s gravitational 
field) and the sustainer phase (concerned with providing low 
thrust at very high specific impulse for long periods of time), 
the devices proposed for solving the second phase of the 
problem (ion rockets, colloid rockets, plasmajets, ete., deriv- 
ing energy from a nuclear reactor) still suffer from the dis- 
advantage that the greater part of the weight of the space 
vehicle must be made up of propulsive fluid.’ Since with 
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present techniques it takes scores of pounds of propellant 
and powerplant to put 1 lb into orbit, even greater economies 
can result if the weight of propulsive fluid required for space 
flight can be eliminated from the total take-off mass. 

The essential feature of this scheme is to lift only the energy 
source into orbit at approximately 100 km, and at that point 
to collect the propulsive fluid (air) for continuing the journey 
into space. This can be accomplished by a Propulsive Fluid 
Accumulator, or PROFAC, and it is a rather simple matter 
to show the startling reduction in take-off mass made possible 
by this approach. The PROFAC can maintain a low altitude 
(100 to 120 km) orbit by making use of a magnetogasdynamic 
(MGD) ramjet to overcome drag. The energy requirements 
of the MGD ramjet and PROFAC equipment can be met by 
relatively low power nuclear generators (1 to 10 mw). 
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Comments on “A Ballistic Bomb Method 
for Determining the Experimental 
Performance of Rocket Propellants” 


R. H. WOODWARD WAESCHE! 


Rohm & Haas Co., Redstone Arsenal Research Division, 
Huntsville, Ala. 


N THEIR article (1)? Messrs. Griffin, Turner and Angeloff 
mention that they adapted the gun propellant concept of 
“force” to give experimental values of specific impulse from 
small amounts of novel materials. The method yields impulse 
without need of thrust data or a known configuration. In 
addition, the authors suggest that tests be conducted with 
this method to determine its usage in the solid propellant field. 
In view of this suggestion, it may interest the authors to 
know that a closed bomb method of determining impulse has 
been in use at the Rohm & Haas Co., Redstone Arsenal Re- 
search Division, since March 1955, and that the Naval Ord- 
nance Test Station and Aerojet-General Corp. have been 
using the experimental method since 1956. 

However, this liquid gun method does present a slightly 
different approach from the one used at this division, aside 
from the fact that it is used here for solid propellants. It has 
been found more convenient and accurate at Rohm & Haas 
to use this test as a relative one, utilizing a “control” propel- 
lant, rather than the direct one used at OMCC. The rela- 
tive method utilizes the fact that, neglecting changes in y (the 
ratio of specific heats), the specific impulse of an unknown 
propellant can be obtained from the relation 

(P nknown (1 sp)unknown 


= [i] 


r 
(PV /m) “control (Iep)control 


Received March 4, 1959. 
! Physicist. 
2 Numbers in parentheses indicate References at end of paper. 
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The pressure term is the maximum pressure attained in the 
closed bomb, corrected for heat loss and igniter pressure. 
The volume term is corrected only for co-volume of the pro- 
pellant gases. Correction for y can be made in the term 


] 
appearing (2) in the equation for 


I 
f(y, P) [2] 


In every case, optimum expansion is assumed, generally to 
14.7 psia from 1000 psig. It may be shown by thermodynamic 
considerations that the closed bomb relative force method 
yields the maximum impulse obtainable from a given formula- 
tion. 

The closed bomb relative force method has been used at 
Rohm & Haas to determine impulses for a wide range of solid 


propellants, including samples from NOTS, Thiokol Chemical 
Corp. (both Redstone and Elkton Divisions), Standard Oil 
Co. (Whiting Lab), Grand Central Rocket Co., Allegany 
Ballistics Laboratory, Jet Propulsion Laboratory, Army 
Rocket and Guided Missile Agency, Aerojet-General Corp., 
and a wide variety of Rohm & Haas propellants. Nearly 
every existing type of solid propellant was tested in this sur- 
vey, made in 1956 through the generous cooperation of the 
participating agencies. Good agreement was found with 
motor firing data. 

Drawings of the bomb and some experimental experience 
may be obtained from the author, including precision and 
methods of correcting for heat loss, which can affect results 
significantly in solid propellants. 
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Book Reviews 


Elastic Waves in Layered Media, by W. 
Maurice Ewing, Wenceslas S. Jardetzky 
and Frank Press, McGraw-Hill Book 
Co., New York, 1957, 380 + xiipp. $10. 


Reviewed by Y1-Yuan Yu 
Polytechnic Institute of Brooklyn 


This monograph written by three 
geophysicists is the outgrowth of a plan 
to make a uniform presentation of the 
investigations on earthquake seismology, 
underwater sound and model seismology 
carried on by the group connected with 
Lamont Geological Observatory of Co- 
lumbia University. Since the book pre- 
sents a systematic treatment of elastic- 
wave problems quite general in nature in 
addition to the specific geophysical 
applications, it should be found useful 
also in other fields as the authors hoped. 
In particular, the verification by experi- 
mental data of the theoretical results 
given here and there throughout the 
book should be of interest to scientists 
and engineers whose work is related to 
elastic-wave propagation and vibration 
problems. 

The treatment is almost exclusively 
based upon the Theory of Linear Elas- 
ticity. Thus, the fundamental equations 
of small motion in an elastic body (reduc- 
ible to the special case of an ideal fluid) 
together with some general solutions are 
briefly presented in the first chapter 
(23 pp.). Chapter 2 (50 pp.) is concerned 
with some classical topics on the homo- 
geneous and isotropic half space, such as 
the reflection of harmonic plane waves 
at the free surface, Rayleigh’s theory of 
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surface waves and Lamb’s integral solu- 
tions for an impulsive force applied along 
a line or at a point. Methods of evalua- 
tion of the integrals in Lamb’s solutions 
are described in detail, and experimental 
verification of these solutions through the 
use of models is finally mentioned in 
this chapter. The propagation of plane 
waves and waves from an impulsive 
point source are again discussed in 
chapter 3 (50 pp.) for two semi-infinite 
media separated by a plane interface. 
Here the theory and experiment related 
to the problem of “refraction arrival’ 
in seismology are discussed. A general 
n-layered half space is then dealt with 
in chapter 4 (131 pp.), in which the 
wave propagation due to an impulsive 
point source is given for a two- or three- 
layered liquid half space and for a liquid 
layer on a solid bottom. In the same 
chapter, the propagation of harmonic 
plane waves (Rayleigh and Love waves) 
is investigated for a two-layered solid 
half space and for the same when covered 
by an additional liquid layer, and many 
important geophysical problems are dis- 
cussed. Chapter 5 (26 pp.) considers the 
effects of gravity, curvature and viscosity 
(through the use of the Voigt and Max- 
well models for a_ viscoelastic solid). 
Chapter 6 (47 pp.) is on plates and 
cylinders, in which the well-known prob- 
lem of Rayleigh and Lamb on the 
vibration of a plate and that of Poch- 
hammer on the vibration of a cylindrical 
rod are discussed together with other 
problems. The last chapter (37 pp.) is 
concerned with wave propagation in 


heterogeneous isotropic media with vari- 
able velocity. Two appendixes (14 pp.) 
deal with the method of steepest descent 
and Rayleigh’s principle. At the end 
of each chapter is a very comprehensive 
list of references on the topics discussed, 
which was compiled from the world 
literature, including the Russian, with a 
total of more than 600 entries. 

It is only natural that the selection 
of topics reflects strongly the interest of 
the geophysicist. Seven-eighths of the 
book is devoted to semi-infinite bodies. 
Although the remaining one-eighth of the 
book (chapter 6) is on plates and cyl- 
inders, these also invariably extend to 
infinity. The authors did mention finite- 
size plates and rods, such as on page 281, 
“The plate can have either finite or 
infinite dimensions, and we now restrict 
ourselves to the latter case,’ and on 
page 312, “To extend the results men- 
tioned here to vibrations of a finite rod, 
additional boundary conditions at the 
ends must be taken into account.’’ These 
however do not reveal what would happen 
if such finite-body problems were to 
be solved on the basis of the exact Theory 
of Elasticity; to satisfy all the boundary 
conditions of a problem would be found 
extremely difficult. Such difficulty to a 
very large extent may be avoided if the 
problems are solved by using approximate 
theories, which are available for plates, 
rods and cylindrical tubes (such as 
those given in references 53, 54, 84 and 
87 quoted in chapter 6 of the book) 
and are accurate enough for some practical 
purposes. In spite of its limited use in the 
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vibration studies of finite bodies, the 
exact Theory of Elasticity is nevertheless 
still important in that the accuracy of an 
approximate theory may be checked by 
comparing the results obtained from both 
theories for certain problems of bodies 
extending to infinity. 

Another comment may be made that 
plates and cylinders composed of more 
than one solid layer (such as those used 
in sandwich construction) are not con- 
sidered in the book, although the wave 
propagation and vibration problems of 
such bodies may readily be formulated 
by the use of the exact Theory of Elas- 
ticity. When such layered plates and 
cylinders are of finite dimensions, it is 
again desirable to use approximate the- 
ories. Such theories for layered plates, 
ior instance, may be developed by em- 
ploying the method given in reference 53 
of chapter 6. 

Needless to say, the comments in the 
last two paragraphs are from a viewpoint 
other than that of a geophysicist and are 
meant to supplement rather than to 
criticize the book. In fact, the reviewer 
feels that the book is so admirably 
written that it is a definite contribution 
to the literature on elastic-wave propaga- 
tion. 


Fluid Dynamic Drag, by Dr.-Ing. S. F. 
Hoerner, published by the author, 
Midland Park, N. J., 1958, xiv + 402 
pp. $14.50. 

Reviewed by Ricwarp F. 
Northwestern University 


This book is a considerably revised and 
enlarged successor to ‘Aerodynamic 
Drag,’ published by the author seven 
years ago. The change in title is oc- 
casioned by the addition of two chapters 
concerning hydrodynamic resistance and 
the characteristics of waterborne craft. 
Another new chapter deals with wind 
forces on structures. The remainder of 
the book, dealing with aerodynamic drag 
of aircraft, missiles, landborne vehicles 
and other objects, is largely unchanged 
in organization but considerably en- 
larged by the incorporation of recent 
data, particularly in the transonic and 
supersonic speed ranges. The new book 
is considerably easier to read than the 
earlier edition, being improved in quality 
of printing and clarity of presentation, 
particularly in the figures. 

The product of a dedicated effort in 
collecting, correlating and evaluating 
data, the book compactly presents a 
wealth of information on fluid dynamic 
resistance. The coverage is_ broad, 
including treatments of such diverse 
subjects as drag of aircraft, automobiles, 
railroad trains, the human body, para- 
chutes, hydrofoils, bluff bodies in hyper- 
sonic flow and drag in rarefied atmospheres. 
Primary reliance is made on experimental 
data with theory used for extrapolation 
and correlation. In many instances, 
the methods of correlating and presenting 
the data are original. The sources of the 
author’s information are carefully ref- 
erenced and include virtually all appli- 
cable American and German publications 
through 1957. 

In assessing the book, it must be 
considered that it is the only one of its 
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kind available. For a person unfamiliar 
with the volume, it is not especially 
easy to use as a handbook. However, 
according to the author, it is not intended 
to be. Instead, it is designed to provide 
sufficient information for guiding the 
user to a solution of his problem. Perusing 
the volume, one feels that there is an 
overemphasis on the aerodynamic charac- 
teristics of obsolete configurations, suc! 
as the airplanes of the 1930’s and earl, 
1940’s. In spite of this, the book repre 
sents a collection of information valuabk 
to many engineers. For the practicing 
aerodynamicist, the data compilation and 
evaluation of this book can save con- 
siderable time and effort in making per- 
formance or loading computations. Yet 
it is written simply enough that any 
engineer concerned with objects moving 
through fluids can make use of the in- 
formation. Thus, although not considered 
suitable as a textbook, ‘‘Fluid Dynamic 
Drag” is an important reference work 
recommended to the many _ engineers 
concerned with drag forces. 


Advanced Dynamics, by John E. Younger, 
The Ronald Press| Company, New 
York, 1958. $8.50. 

Reviewed by Burton 
University of Missouri 


Persons who show a real facility with 
the principles first covered in introductory 
engineering mechanics courses are rare. 
Their skill is seldom, if ever, acquired 
without further study, since it calls for 
an overall perspective of mechanics 
along with a proficiency in the several 
parts of the subject. Younger has written 
a textbook that will be welcomed by 
teachers attempting to provide such a 
course of study. 

In rather broad chapters, he deals with 
kinematical analysis, dynamic specifica- 
tion of a rigid body, dynamics of angular 
motion, dynamics of plane motion, 
dynamics of motion in three dimensions, 
dynamics of beams. He includes discus- 
sion of elastic phenomena, Castigliano’s 
theorem, and a pretty thorough intro- 
duction to vibration. Problems appear 
throughout the text. Younger’s method 
of presentation is to break the subject 
down into parts and deal with these in 
some detail. In some ways this presenta- 
tion leaves me unsatisfied. Often he 
seems to make several cases out of what 
might more logically be thought of as 
different aspects of a single case. On 
the other hand, this leaves little to be 
supplied by the student and may be the 
best way to prepare him to use his me- 
chanics. Another objection is concerned 
with the author’s reluctance to include 
discussion of the meaning of what is done. 
For example, there is very little discussion 
of Newton’s laws of motion, and I do 
not believe it is ever made clear just 
why laws 1 and 3 are needed. I feel 
that additional discussion is needed also 
to clarify the presentation of angular 
momentum and the variable mass prob- 
lem. 

These objections should not be allowed 
to detract from the choice of material, 
the well worked out derivations and the 
suitability of the book for a one-semester 
course for engineers. 
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New Patents 


Production of ammonium pentaborate and 
boric oxide (2,867,502). H. Stange and 
S. L. Clark, Tonawanda, N. Y., assignor 
to Olin Mathieson Chemical Corp. 

Method consists of thermally decompos- 
ing ammonium pentaborate from borax. 
Guided missile design (184,208) and Ram- 
jet missile design (184,209). A. J. Bell, 
C. W. Besserer and E. A. Bonney, W. 
Hyattsville, Md., assignors to the U. 8. 
Navy. 


Developed at the Applied Physics 
Laboratory, Johns Hopkins University, 
Silver Spring, Md. in 1949. The guided 
missile, now in service, is the Terrier; the 
ramjet missile, about to go into service, is 
the Talos. Both missiles have undergone 
modifications since the original designs 
shown in the accompanying illustrations. 


Epitor’s Nore: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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Navigation device (2,867,393). R. M. 
Burley, Syosset, N. Y., assignor to Bendix 
Aviation Corp. 

A member periodically interrupting 
rays from a celestial body, and cooperating 
with a sensitive means to produce signals 
when the energy received varies. Revers- 
ible motors operate azimuth and elevation 
controls in response to the signals and 
reference voltages to maintain the craft on 
a fixed course. 

Method for the separation of a 
(2,867,498). J. Roscoe and G. 
Schaeffer, St. Louis, Mo., assignors Re 
Olin Mathieson Chemical Corp. 

Process consists of contacting a mixture 
of lithium aluminum hydride and alumi- 
num hydride with a molar amount, based 
on the moles of the hydrides present, of a 
lower tertiary alkyl amine while the 
reactants are in admixture with a lower 
dialkyl. The precipitate is separated from 
the reaction mixture. 


Conversion of sodium tetrachloroalumi- 
nate to sodium aluminum hydride (2,867,- 
499). C.D. Good and H. B. Batha, Grand 
Island, N. Y., assignors to Olin Mathieson 
Chemical Corp. 

Process consists of reacting the two 
products with a lower dialky] ether main- 
tained in liquid phase, to form an addition 
product of aluminum chloride, and react- 
ing the addition product with sodium hy- 
dride. 

Self-stabilized burner (2,867,977). R. 5S. 
Buck, Glastonbury, Conn., assignor to 
United Aircraft Corp. 

Ramjet power slant having means for 
automatically adjusting a blocking mem- 
ber in any one of various positions. A 
mechanism senses pressure of gases 
approaching the combustion chamber, 
proportionately blocking their passage. 


Auxiliary rotor airfoil (2,868,302). A. C. 
Peterson, Minneapolis, Minn. 

Jet units, at a substantial distance from 
the axis of a mounting pylon, supplied by 
fluid under pressure from a chamber con- 
trolled by a pressure responsive element. 


Apparatus for igniting fuels (2,867,979). 
J. W. Mullen II, Richmond, Va., assignor 
to Experiment, Inc. 

Ramjet with flameholder comprising an 
annular shroud with apertures to the mani- 
fold formed by the hollow of the shroud. 


Aerodynamic wave machine (2,867,981). 
M. Berchtold, Paoli, Pa., assignor to 
I-T-E Circuit Breaker Co 

Machine to produce shaft power for a 
compressor and turbine, operating on a 
forward cycle of instantionary flow phe- 
nomenon. A first stator plate has a hot 
exhaust port and a cold air pick-up port. 
A second stator plate has a cold air intake 
port and a hot intake nozzle. The pick-up 
port has a higher pressure than the nozzle, 
causing a reflected compression wave. 


engines (2,867,982). J. Clarke and 
E. Collinson, W Ontario, 
Canada, assignors to Joseph Lucos (Indus- 
tries) Ltd 

Baffle with openings through which air 
from an annular entrance can flow between 
burner nozzles. Plates mounted on the 
approach side of the baffle set up a swirling 
motion in the air passing through the 
openings. 
Gyroscopic a tus (2,868,021). H. B. 
Sedgfield a he Luttrelle, Holyport, 


George F. McLaughlin, Contributor 


England, assignors to The Sperry Gryo- 
scope Co. Ltd. 

Electromagnets, fixed to a gimbal fraine 
support, having air gaps across which flux 
passes radially in the plane of the fraie 
axes displaced at an angle to each frame 
axis. An energized coil connected to the 
rotor case is bowed outwardly in opposite 
positions to thread the gaps in the electro- 
magnets. 

Device for transforming rotary movemen 
to linear movement (2,868,030). 
Forwald, Ludvika, Sweden, assignor to 
Allmanna Svenska Elektriska Aktiebo!a- 


get. 

Externally screw-threaded sleeve 
mounted to rotate with a driven shaft, and 
having limited axial movement. A nit 
engaging the threads is attached to a 
member extending beyond the housing ‘0 
transmit linear motion from the rotary 
motion of the driving shaft. 

Adjustable mirror support in successive 
oT monochrometer (2,868,063). 

. Weiss, Philadelphia, Pa., assignor to 
Leeds and Northrup Co. 

System including radiant energy dis- 
persing means and a pair of angularly 
related plane mirrors for causing the 
energy to make another pass through the 
dispersing means. Mirrors are supported 
so as to provide a fine focus adjustment of 
the dispersed energy. 

Jet propelled vertical take-off aircraft 
(2,868, 477). H. E. Chaplin, Hayes, 


England, assignor to The Fairey Aviation 
Co. Ltd. 


Jet engine capable of giving a total 
static thrust exceeding the weight of the 
aircraft. The undercarriage is movable to 
a position in which it supports the fuselage 
in a position for vertical takeoff. Power 
driven means, controllable by the pilot, 
moves the cockpit so the pilot remains in a 
horizontal seated position during changes 
of the angularity of the fuselage relative to 
the ground. 
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Now from Du Pont: 


a complete line of photorecording papers and films 


Du Pont now makes photorecording papers and films for 
your various recording needs. Here’s what our complete 
photorecording line includes: 


PHOTORECORDING PAPERS 
Lino-Writ 1. Orthochromatic paper designed for average test 
requirements at lower writing speeds. 
Lino-Writ 2. Moderately high-speed orthochromatic paper 
for use in intermediate speed recording requirements. 
Lino-Writ 3. For high-frequency oscillographic traces with- 
out loss of detail. 
These three papers are available in two types: Type B, stand- 
ard-weight stock when greater opacity is required, and Type 
W, a thin, 100% all-rag stock. 
Lino-Writ 4. It’s the toughest, fastest, whitest paper you can 
use; 20% thinner than Type W, accommodates writing speeds 
as high as 5000 cps at maximum amplitude. Exceptional 
exposing and processing latitude. All-rag stock has unusual 
wet or dry strength. 


*Caonar is Du Pont's trademark for its polyester photographic films 


This advertisement was prepared exclusively by Phototypography 
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PHOTORECORDING FILMS 

CRONAR* Recording Film (High Speed) —relative speed 
(tungsten): 30. 

Lino-Flex 1 (Regular Speed)—relative speed (tungsten): 10. 
Both of these films are on Du Pont’s exclusive CRONAR 
polyester photographic film base, which offers unexcelled 
strength, exceptional dimensional stability and flexibility. 
These fast-drying films have a fine matte surface which 
accepts pencil and ink markings. 


PHOTORECORDING CHEMICALS 


Lino-Writ Rapid Processing Chemicals Kit. Du Pont’s 
high-speed, low-odor chemicals designed for rapid stabiliza- 
tion processing of papers and films. 


For more information on our complete photorecording line, 
write: E. I. du Pont de Nemours & Co. (Inc.), Photo Prod- 
ucts Department, Wilmington 98, Delaware. In Canada: 
Du Pont Company of Canada Limited, Toronto. 
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Senior staff 
openings 
in basic 
research 


THE BOEING SCIENTIFIC RESEARCH 
LABORATORIES are engaged in a 
program of fundamental research 
designed to make contributions to 
the progress of the physical 
sciences. High-level staff positions 
are open now in the fields of 


Gas Dynamics 
Plasma Physics 
Mathematics 

Solid State Physics 
Electronics 
Physical Chemistry 
Geoastrophysics 


Boeing grants scientists the lati- 
tude and independence needed to 
achieve and maintain leadership in 
their special fields. Scientists inter- 
ested in carrying on their work in 
this kind of stimulating research 
environment are invited to com- 
municate with Mr. G. L. Hollings- 
worth, Associate Director, Boeing 
Scientific Research Laboratories. 


P.O. Box 3822 -JTA 
Boeing Airplane Company 
Seattle 24, Washington 
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Technical Literature Digest 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Jet Propulsion Engines 


Rocket Motor Impulse Measurements 
by Means of an Inertia Wheel, by Ben- 
jamin V. Connor, Calif. Inst. Tech., Jet 
Propulsion Lab., Mem. no. 20-167, 
Sept. 1958, 12 pp. 

The Research Combustor of the In- 
stituto di Macchine, Politecnico di Torino, 
by F. Filippi, L’Aerotecnica, vol. 38, 
April 1958, pp. 88-95. (In Italian.) 

On the Influence of the Form of Rocket 
Combustion Chambers, by L. Poggi, 
L’ Aerotecnica, vol. 38, April 1958, pp. 
109-110. (In Italian.) 

Flame Stabilization in a Ramjet Com- 
bustion Chamber, by G. Salvatore, 
L’ Aerotecnica, vol. 38, April 1958, pp. 
111-116. (In Italian.) 

Upper Atmosphere Atomic-oxygen 
Power Plant, by Sterge T. Demetriades, 
os Aero/Space Sci., vol. 25, Oct. 1958, p. 

3. 


Results of Theoretical Rocket Motor 
Performance Calculations for Liquid Pro- 
pellant Systems, by G. W. Elverum Jr., 
J. S. Martinez and R. W. Gillespie, Calif. 
Inst. Tech. Jet Propulsion Lab. Mem. no. 
20-159, Nov. 1957, 34 pp. 

Combustor Performance with Various 
Hydrogen-oxygen Injection Methods in a 
200-Pound-thrust Rocket Engine, by M. 
F. Heidmann and Louis Baker Jr., 
NACA RM no. E58E21, Sept. 1958, 45 


pp- 


Aerodynamics of Jet 
Propelled Missiles 


General Formulation of Powered Flight 
Trajectory Optimization Problems, by 
Burton D. Fried, J. Appl. Phys., vol. 29, 
Aug. 1958, pp. 1203-1209. 

Effect of Favorable Pressure Gradients 
on Transition for Several Bodies of Revo- 
lution at Mach 3.12, by John R. Jack, 
NACA TN 4313, July 1958, 28 pp. diagrs., 
photo. 

Correlation by the Hypersonic Similarity 
Rule of Pressure Distributions and Wave 
Drags for Minimum-drag Nose Shapes at 
Zero Angle of Attack, by Leland H. Jor- 
genson, NACA Res. Mem. A53F12, Aug. 
1953, 23 pp., diagrs. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 128, 7/28/58.) 

Angle of Attack Convergence of a Spin- 
ning Missile Descending through the 
Atmosphere, by Herman L. Leon, J. 
Aero/Space Sct., vol. 25, Aug. 1958, pp. 

480-484. 


The Supersonic Blunt-body Problem— 
Review and Extension, by Milton D. Van- 
Dyke, J. Aero/Space Sci., vol. 25, Aug. 
1958, pp. 485-496. 

The Flow of a Supersonic Jet in a Super- 
sonic Stream at an Angle of Attack, by F. 
Edward Ehlers and Torstein Strand, J. 
Aero/Space Sci., vol. 25, Aug. 1958, pp. 
497-503. 


Eprtor’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged.’ 


Heat Transfer and 
Fluid Flow 


Flow Visualization Techniques Applied 
to Combustion Problems, by E. F. Winter, 
J. Roy. Aeron. Soc., vol. 62, no. 568, April 
1958, pp. 268-276. 

Variable Fluid-pro in 
Free Convection, by Sparrow and 
J. L. Gregg, ASME Ft 57-A—46, Dee 
1-6, 1957, 8 pp. 

The Thermal Conductivity of Nitrogen 
and Argon in the Liquid and Gaseous 
States, by H. Ziebland and J. T. A. Bur- 
ton, Brit. J. Appl. Phys., vol. 9, Fcb. 
1958, pp. 52-59. 

Thermal Resistance of Eddy, by L. (;. 
Clark and W. W. Hagerty, J. AIChE, vol. 
3, no. 4, Dec. 1957, pp. 523-527. 

On the Decay of Homogeneous Tur- 
bulence before the Final Period, by !:. 
Deissler, Phys. Fluids, vol. 1, March- 
April 1958, pp. 111-121. 

Combined Effects of Turbulence and 
Roughness on Transition, by H. Dryden, 
ZAMP, vol. 9B, 1958, pp. 249-258. 


Heat Transfer to Turbulent Boundary 
Layer Downstream of Injection Slot, hy 
R. A. Seban, H. W. Chan and S. Scesi:, 
ASME Paper 57—A-36, Dec. 1-6, 1957, 
11 pp. 

Mass Transfer in Liquid Metal Sys- 
tems, I: Thermal Gradient Transfer, by 
J. W. Taylor, Nuclear Power, vol. 3, Fel). 
1958, pp. 53-57. 

Some Effects of Kinetic Heating on the 
Stiffness of Thin Wings—Analysis of 
Flexural and Torsional Stiffnesses for 
Large Deformation of Thin Solid Wings, 
by 8S. L. Kochanski and J. H. Argyris, 
Aircr. Engng., vol. 30, March 1958, pp. 
82-85. 

Numerical Determination of Transitory 
Temperatures in Wings, by G. Lehr, 
Rev. Document. de la Tech. Aeron. Mon- 
diale, no. 48, Jan. 1958, pp. 3-14. 

The Effects of Kinetic Heating on Air- 
craft Structures, by A. W. Kitchenside, 
J. Roy. Aeron. Soc., vol. 62, Feb. 1958, 
pp. 105-117. 

An Experimental Study of the Effect of 
Mass Injection at the Stagnation Point of 
a Blunt Body, by Howard M. McMahon, 
Calif. Inst. Tech., Guggenheim Aeron. Lab., 
Hypersonic Res. Proj., Mem. 42, May 
1958, 132 pp. 

The Thermodynamics of Steady Con- 
duction-diffusion-reaction States of Multi- 
component Fluid Systems, by W. Byers 
Brown, Trans. Faraday Soc., vol. 54, June 
1958, pp. 772-784. 

Thermal Diffusion in Liquids, by Ste- 
phen Whitaker and R. L. Pigford, Ind. & 
Engng. Chem., vol. 50, July 1958, pp. 
1026-1032. 


Calculation of Reaction Profiles behind 
Steady-state Shock Waves, I: ““_re" 
to Detonation Waves, by Russell E. Duff, 
J. Chem. Phys., vol. 28, June 1958, pp. 
1193-1197. 

High-temperature Gas Viscosities, I: 
Nitrous Oxide and Oxygen, by C. J. G. 
Raw and C. P. Ellis, J. Chem. Phys., vol. 
28, June 1958, pp. 1198-1200. 


“ Ettingshausen Effect and Thermomag- 
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‘MEETS ITS SEVEREST 
SGHALLENGE IN NUCLEAR POWER 
PEANTS FOR FLIGHT 


etaliurgi 


High strength at high temperatures — for long service missions (30 days of continuous 


flight and more are envisaged) is a primary consideration in the selection and develop- 
ment of materials for structures and components of General Electric’s 2nd and 3rd 
generation nuclear propulsion systems for flight application. 


A complex of additional parameters, however, must be taken into consideration by 
G-E engineers and scientists ... radiation effects, nuclear properties for specific appli- 


cation, stringent weight limitations, etc. 


Engineers and Scientists who value the opportunity to do original work with a company 
that fosters free inquiry and initiative, are invited to inquire about positions now open 
in both AppLiED RESEARCH and MATERIALS DEVELOPMENT: 


DIRECT processing and procurement of basic alloys 
for elevated temperature applications. PhD, MS, 3-7 years 
experience. 


CONDUCT liai with technical production and project 
personnel on problems related to quality in fuel element 
production. PhD, MS, 3-7 years experience. 


ADVISE Materials Testing Lab on materials problems 
and techniques. Direct experimental, physical, mechani- 
cal and electrical tests. PhD, MS, 3-7 years experience. 


CARRY ON investigations of reactions in the solid 
state for Materials Testing Laboratories. PhD, Physical 
Chemistry. 


PARTICIPATE in research on advanced moderator ma- 
terials. Investigations in field of ionic, covalent and 
metallic bonding; hydrogen stability in metals at ele- 
vated temperatures; kinetics of advanced moderator 
systems. MS. 


PERFORM laboratory investigations in metallurgical 
melting, casting, fabricating and forming processes. Also 
phase studies; solid solution alloy studies. MS preferred. 


DIRECT studies on materials used in moderator mate- 
rials development. PhD. 


PROGRAM and evaluate specific tests on fuel element 
materials and assemblies. Compile and publish data. 
MS, BS, 5-7 years R&D experience. 


CONDUCT project technical liaison on materials selec- 
tion and processing. Develop design data on materials. 
BS, 5-7 years experience with high temperature alloys. 


CARRY OUT liaison with materials subcontractors. BS 
with 7-9 years experience in development or sales engi- 
neering. 


ASSIST in development of instrumentation for measur- 
ing high temperatures and other instrumentation. MS, 
BS, 1-3 years experience. 


If you qualify, please write in confidence, 
including salary requirements, to: Mr. P. W. Christos, Div. 35-ME. 
AIRCRAFT NUCLEAR PROPULSION DEPARTMENT 


GENERAL ELECTRIC 


P.O. Box 132 « Cincinnati 15, Ohio 


Physical Chemists, 
emical Engineers, PhD, MS, BS 
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THROUGHOUT 
BOLD MINDS have sought to understand 


the forces at work in the universe, and as they devel- 
oped working hypotheses, endeavored to turn all 
knowledge to their own purposes, devising philosophi- 
cal and mechanical systems of their own. 

As old hypotheses become inadequate or unten- 
able, thinking men devise new ones. So the concept 


ENGINEERS — SCIENTISTS 


Is order and organization in 
MACROCOSM and MICROCOSM an 
underlying principle of the universe... 
or does a fundamental state of 
randomness prevail in the cosmos. 

as might be expressed 
by the multinomial 
formula at left? 


of a “flat” world has changed to an oblate orbiting 
spheroid—mere speck in a vast and expanding uni- 
verse; so “empty” formless space is regarded as a 
curved continuum occupied by random knots of tur- 
bulence (creating the new branch of mechanics — 
hydromagnetics). 

Today new horizons of discovery and surmise arise 
before the speculative mind. 


ENGINEERS AND SCIENTISTS AT REPUBLIC FEEL KINSHIP WITH ALL BOLD MINDS OF PAST AND 
PRESENT, AS THEY FACE THE EXHILARATING CHALLENGES OF CREATING VEHICLES TO FLY IN 
ENVIRONMENTS WHERE NEW APPROACHES IN THERMO/AERODYNAMICS MUST BE MADE...AS WELL 
AS APPROPRIATE PROPULSION AND ELECTRONIC SYSTEMS TO POWER AND GUIDE TRANSIT IN SPACE 


Groups of Research, Development and Experimental Engineers and Scientists at Republic Aviation 
are now working on projects over the whole range of aeronautics and astronaulics—from supersonic and 
hypersonic weapons systems, both manned and unmanned, to plasma propulsion and space electronics. 


Imaginative professional men at many levels of experience are invited to inquire 
about opportunities indicated below: 


ELECTRONICS THERMO, PLASMA PROPULSION 

Inertial Guidance & Navigation AERODYNAMICS Plasma Physics 

Digital Computer Development Theoretical Gasdynamics Gaseous Electronics 

Hyper-Velocity Studies and Shock 
rmation ry enomena 

Telemetry-SSB Technique Hydromagnetics 

Doppler Radar Ai Kae fi Physical Chemistry 

Countermeasures irplane/Missile Performance Combustion and Detonation 

Radome & Antenna Design Air Load and Aeroelasticity Instrumentation 

oe Circuitry & Stability and Controls High Power Pulse Electronics 
components Flutter & Vibratio 

Receiver & Transmitter Design Vehicle Dyn pana System NUCLEAR PROPULSION 

Airborne Navigational Systems Designs & RADIATION PHENOMENA 

SES SSE High Altitude Atmosphere Nuclear Weapons Effects 

Miniaturization- Radiation Environment 

Transistorization in Space 

Ranging Systems e-entry Heat Transfer Nuclear Power & Propulsion 

Propagation Studies Hydromagnetics Applications 

Ground Support Equipment Ground Support Equipment Nuclear Radiation Laboratories 


A new $14,000,000 Research Cen- 
ter —to be completed this year — 
is part of Republic's far-ranging 
R&D programs aimed at major 
state-of-the-art breakthroughs in 


Send resume in complete confidence to: 
Mr. George R. Hickman, ENGINEERING EMPLOYMENT MANAGER — DEPT. 10E 


STEP 


every flight regime & environment. FARMINGDALE, LONG ISLAND, NEW YORK 
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netic Cooling, by B. J. O’Brien and C. S. 
Wallace, J. Appl. Phys., vol. 29, July 
1958, pp. 1010-1012. 

Electrically Driven Shock Tube, by D. 
E. Bloxsom Jr., J. Appl. Phys., vol. 29, 
July 1958, pp. 1128-1129. 

Flows in Partly Dissociated Gases, by 
Manfred Heil, J. Aero/Space Sct., vol. 25, 
July 1958, pp. 459-460. 

The Injection of Air into the Dissociated 
Hypersonic Laminar Boundary Layer, by 
Sinclaire M. Scala, J. Aero/Space Sci., 
vol. 25, July 1958, pp. 461-462. 

On the Aerodynamic Heating of Blunt 
Bodies, by E. R. van Driest, ZAMP, 
vol. 9B, March 25, 1958 (Jakob Ackeret 
anniv. vol.), pp. 233-248. 


Combustion, Fuels and 
Propellants 


Calculated Heats of Formation and 
Combustion of Boron Compounds (Boron, 
Hydrogen, Carbon, Silicon), by Aubrey P. 
Altshuller, NACA Res. Mem. E55G26, 
Oct. 1955, 28 pp. 

Aluminum Borohydride as an Ignition 
Source for Turbojet Combustors, by 
David M. Straight, Edward A. Fletcher 
and Hampton H. Foster, NACA Res. 
Mem. E53G15, Sept. 1953, 19 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 129, 7/17/58.) 

The Heat of Combustion of Tetraethyl- 
diborane, by Stanley Tannenbaum, NACA 
Res. Mem. E53E11, Jan. 1957, 5 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 129, 7/17/58.) 

Fundamental Flame Velocity, Heat of 
Combustion, and Spontaneous Ignition 
Temperature of Dimethylaminodiborane, 
by Robert R. Hibbard and Gordon L. 
Dugger, NACA Res. Mem. F52L23, Dec. 
1956, 6 pp. (Declassified from Confiden- 
tial by authority of NACA Res. Abstr. 129, 
7/17/58.) 

Continuous Combustion Problems in 
Aero Engines, by A. Capetti, L’ Aerotec- 
nica, vol. 38, Feb. 1958, pp. 3-7. (In 
Italian. ) 

Flame Stabilization in Air-breathing 
Jet Engines, by F. Filippi, L’ Aerotecnica, 
vol. 38, Feb. 1958, pp. 19-38. (In Ital- 
ian. ) 

Liquid Propellants for Rocket ae 
by E. Macioce, L’ Aerotecnica, vol. 
Feb. 1958, pp. 44-47, (In Italian. ) 

On Turbulence in Flames, by A. G. 
Prudnikov, Akad. Nauk SSSR, Izvestia, 
Otdel. Tekh. Nauk, no. 7, July 1958, pp. 
130-133. 

Preparation and Thermal Stability of 
Lithium Titanium Fluoride, by George J. 
Janz, Max R. Lorenz and Charles T. 
Brown, J. Amer. Chem. Soc., vol. 80, Aug. 
20, 1958, pp. 4126-4128. 

Detonation Velocity of Gaseous Ozone, 
by A. G. Streng, C. 8. Stokes, and L. A 
Streng, J. Chem. Phys., vol. 29, Aug. 1958, 
pp. 458-459. 

The Ignition of Primary Explosives by 
Electric Discharges, by R. M. H. Wyatt, 
P. W. J. Moore, G. K. Adams and J. F. 
Sumner, Proc. Roy. Soc., London, vol. A 
246, July 29, 1958, pp. 189-196. 

The Explosion of Silver Azide in an 
Electric Field, by F. P. Bowden and 
A. C. Melaren, Proc. Roy. Soc., London, 
vol. A 246, July 29, 1958, pp. 197-198. 

Structure and Stability of Inorganic 
Azides, by B. L. Evans, Proc. Roy. Soc., 
London, vol. A 246, July 29, 1958, pp. 
199-202. 

Thermal Decomposition, Inflammation, 
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and Detonation, by W. E. Garner, Proc. 


Roy. Soc., London, vol. A 246, July 29, 
1958, pp. 203-215. 

The Initiation of Explosion by Neutrons, 
A-particles and Fission Products, by F. P. 
Bowden, Proc. Roy. Soc., London, vol. A 
246, July 29, 1958, pp. 216-218. 

The Effect of Nuclear Radiation on Ex- 
plosives, by J. V. R. Kaufman, Proc. Roy. 

oc., London, vol. A 246, July 29, 1958, pp. 
219-232. 

Preparation and Properties of ’Crystals 
of KN; and NaN;, by R. W. Dreyfus and 
P. W. Levy, Proc. Roy. Soc., London, vol. A 
246, July 29, 1958, pp. 233-239. 

The Initiation of Explosion by Light, by 
J. Eggert, Proc. Roy. Soc., London, vol. 
A 246, July 29, 1958, pp. 240-247. 

Ignition of Explosives by Light, by 
J.H.L. McAuslan, Proc. Roy. Soc., London, 
vol. A 246, July 29, 1958, pp. 248-253. 

The Growth of Explosion in Solids, by 
A. D. Yoffe, Proc. Roy. Soc., London, vol. 
A 246, July 29, 1958, pp. 254-256. 

Some Considerations on the Mechanism 
of Initiation of Detonation Explosives, by 
K. K. Andreev, Proc. Roy. Soc., London, 
vol. A 246, July 29, 1958, pp. 257-267. 

The Initiation of Explosives by Shock, 
by G. P. Chachia and E. G. W hitbread, 
Proc. Roy. Sec., London, vol. A 246, July 
29, 1958, pp. 268-273. 

Sympathetic Detonation and Initiation 
by Impact, by R. J. Eichelberger and M. 
Sultanoff, Proc. Roy. Soc., London, vol. A 
246, July 29, 1958, pp. 274-280. 

Deflagration to Detonation Transition in 
Solid and Liquid — by M. A. 
Cook, D. H. Pack and W. A. Gey, Proc. 
Roy. ’Soc., London, vol. M 246, July 29, 
1958, pp. 281-283. 

The Initiation of Condensed Explosives 
by Shock Waves from Metals, by W. R. 
Marlow and I. C. Skidmore, Proc. Roy. 
Soc., London, vol. A 246, July 29, 1958, 
pp. 284-287. 


Materials of Construction 


Fatigue of Structural Materials at 
High Temperatures, by B. J. Lazan, 
NATO, AGARD, Rep. 156, Nov. 1957, 
27 pp. 

Influence of Hot-working Conditions 
on High-temperature Properties of a 
— Alloy, by John F. Ewing 

W. Freeman, NACA Rep. 1341, 
1957, 52 pp. (supersedes TN 3727). 

Limit Design for Economical Missile 
Structures, by L. A. Reidinger, Aviation 
Age, vol. 30, Oct. 1958, pp. 32-39. 

301 Stainless Modified for 800 Deg F 
and Up, by D. B. Roach, A. M. Hall, BR. 
Sergesen and A. R. Stargardter, Aviation 
Age, vol. 30, Oct. 1958, pp. 58-63. 

Friction and Wear with Reactive Gases 
at Temperatures up to 1200° C, by Gordon 
P. Allen, Donald H. Buckley and Robert L. 
Johnson, NACA TN 4316, Sept. 1958, 
26 pp. 

A Phenomenological Relation between 
Stress, Strain Rate, and Temperature for 
Metals at Elevated Temperatures, by 
Elbridge Z. Stowell, NACA Rep. 13438, 
1958, 6 pp. (supersedes TN 4000). 

Structures and Heat at Supersonic and 
Hypersonic Speeds, by . Dorleac, 
NATO, AGARD, Rep. 149, Nov. 1957, 
78 pp. (In French.) 

On Ceramic-metallic Materials, Ce- 
ramic Materials and High Temperature 
Protecting Lines for Jet Engines, by R. M. 


ENGINEERS 
SCIENTISTS 


Here is your opportunity to grow 
with a young, expanding subsidiary 
of the Ford Motor Company. Out- 
standing career opportunities are 
open in Aeronutronic’s new RE- 
SEARCH CENTER, overlooking the 
Pacific at Newport Beach, and the 
facility in Glendale, California. 
You will have all the advantages of 
a stimulating mental environment, 
working with advanced equipment 
in a new facility, located where you 
_ enjoy California living at its 
inest. 


PhD and MS RESEARCH SPECIAL- 
ISTS with 5 to 7 years’ experience in heat 
transfer, fluid mechanics, thermodynamics, 
combustion and chemical kinetics, and 
thermoelasticity. To work on theoretical 
and experimental programs related to re- 
entry technology and advanced rocket pro- 
pulsion. Specific assignments are open 
in re-entry body design, high temperature 
materials studies, bound ary layer heat 
transfer with erm: 
stress analysis, and high temperature 
thermodynamics. 


PROPULSION ENGINEERS with 5 

experience in liquid and solid rocket 

and test. Familiarity with heat st 
bl in desirable. To work 

on program of wide scope in R & D of ad- 

vanced concepts in rocket engine com- 

ponents, and for missile project work. 


APPLIED MATHEMATICIANS, 3-5 

years’ recent experience required and A.B. 

or M. A. degree. Experience in numerical 
and ter work in 

with rockets and rocket propellants. 


STRESS HIGH TEMPERATURE MA- 
TERIALS. Mechanical or Metallurgical 
Engineer, must be familiar with thermal 
stress and shock as well as elasticity and 
plasticity. Application to nose cones and 
thrust chambers. 


FLIGHT TEST & IN: STRUMENTATION 
ENGINEERS with 5 to 10 years’ experi- 
ence in laboratory and flight test instru- 

tec ill d tech- 
niques utilizing advanced instrumentation 
associated with space vehicles. 


THEORETICAL AEROTHERMODY- 
NAMICIST. Advanced degree and at 
least 5 years’ experience in high-speed 
aerodynamics. Knowledge of viscid and 
gas flows required. To work on program 
leading to advanced missile configurations. 
Work involves analysis of the re-entry of 
hypersonic missiles and space craft for 
determining optimum configuration. 


CERAMICIST. M.S. or Ph.D. required 
and 3-5 years’ recent experience with high 
temperature materials, structures and cer- 
mets. 


ENGINEER or PHYSICIST. With ex- 
perience in the use of scientific instruments 
for making physical measurement. Work 
related to flight test and facility instrumenta- 
tion. Advanced degree desired with mini- 
mum of 3 years of related experience. 


Qualified a plicants are invited to send 
ead! to Mr. R. W. Speich, 


Systems, 


AERONUTRONIC 


a subsidiary of Ford Motor Company 


2234 Air Way Bldg. 19, Glendale, Calif. 
CHapman 5-6651 


Newport Beach, Glendale, Santa Ana, 
and Maywood, California 
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Creative people at IBM are hard at work on a problem that has vexed man since the 

dawn of history—communication with his fellow beings. The rapid transmission of precise data 
is vital in this fast-moving age. One group of engineers and scientists at IBM is working 

on multiplexing techniques to channel data. This data has to be made available at maximum 
speed at desired times and places. Another group is developing computer systems that 

will control complete warehousing or industrial operations. For the military, IBM is designing 
a computer that will make available in real time, data concerning troop strength and 
dispositions. This important work needs imaginative people. If you want to be a part of 

this world of discovery and development, there are many opportunities at IBM. 


A NEW WORLD OF OPPORTUNITY. Both technical and administrative engineering careers offer 
parallel advancement opportunities and rewards at IBM. You will enjoy unusual professional 
freedom, comprehensive education programs, the assistance of specialists of diverse disciplines, 
and IBM’s wealth of systems know-how. Working independently or as a member of a small 
team, your individual contributions are quickly recognized and rewarded. This is a unique 
opportunity for a career with a company that has an outstanding growth record. 


For details, write, outlining 
background and interests, to: 
Mr. R. E. Rodgers, Dept. 572E 
IBM Corporation 

590 Madison Avenue 

New York 22, N. Y. 


® 
INTERNATIONAL BUSINESS MACHINES CORPORATION 
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7 OU CAN EXPlOre NEW Areas at umn 


CAREERS AVAILABLE 
IN THESE AREAS... 
CIRCUIT DESIGN & RESEARCH 
COMMUNICATIONS 
COMPONENT ENGINEERING 
COMPUTER ANALYSIS & DESIGN 
CRYOGENICS 
ENGINEERING PHYSICS 
INERTIAL GUIDANCE 
INDUSTRIAL CONTROL 
INFORMATION THEORY 
LOGICAL DESIGN 
MATHEMATICS 
OPTICS 
PROGRAMMING 
RADAR 
REAL TIME STUDIES 
SEMICONDUCTOR DESIGN 
SOLID STATE DEVELOPMENT 
SYSTEMS DEVELOPMENT 


Qualifications: B.S., M.S., or Ph.D. 
in Electrical or Mechanical Engineer- 
ing, Physics, or Mathematics — and 
proven ability to assume a high de- 
gree of technical responsibility in 
your sphere of interest. 


May 1959 


GENERAL DESCRIPTIONS 
OF SOME ASSIGNMENTS: 


SYSTEMS ENGINEERS to analyze and 
design computer systems. Backgrounds 
required include analog to digital con- 
version; analyses and preparation of 
diagnostic programs; development of 
complex devices in servo-mechanisms 
or radar for advanced systems. 


LOGICAL DESIGNERS with two to four 
years’ experience outlining logical 
blocks. Previous computer experience 
desirable. Some experience in transis- 
torized equipment specification helpful. 
Must have strong interest in the theo- 
retical and practical aspects of check- 
ing means and various codes, redun- 
dancy, error detection and correction, 
information flow, and other factors. 


SOLID STATE ENGINEERS AND SCIEN- 
TISTS to do applied research on pre- 
cision linear circuitry employing solid 
state devices using analog to digital 
conversion techniques and sample 
data. Experience in feed-back amplifier 
design desirable. Also opening for 
engineers with experience in precision, 
low-level linear circuits employing 
solid state, to work on analog to digi- 
tal conversion techniques. 


INDUSTRIAL CONTROL ENGINEERS to 
perform precision AC and DC electric 
measurements. Assignments in ampli- 
fier design, relay circuit logic, test 
equipment development, analog to 
digital conversion, and noise reduction. 


CIRCUIT DESIGNERS to design tran- 
sistor amplifiers, relay lines, transistor 
tube conversion circuits. Develop sys- 
tems circuit specifications, perform 
circuit evaluation experiments and reli- 
ability criteria. Other openings in cir- 
cuit design for magnetic devices. 


MATHEMATICIANS to do digital com- 
puter programming, handle analysis- 
of-variance and multiple-regression 
type problems. Design experiments for 
wide variation of engineering applica- 
tions. Knowledge of application of 
probability or game theory desirable. 


Carelli, L’Aerotecnica, vol. 38, no. 2, 
April 1958, pp. 71-87. (In Italian.) 


Instrumentation, Data 
Recording, Telemetering 


Kinematic Drift of Single-Axis Gyro- 
scopes, by R. H. Cannon Jr., J. Appl. 
Mech., vol. 25, no. 3, Sept. 1958, pp. 
357-360. 

Telemetry Standards for Guided 
Missiles, Hlectronics, vol. 31, Oct. 24, 
1958, pp. 96-98. 

Five Hole Spherical Pitot Tube, by P 
C. Pien, David Taylor Model Basin, Rep., 
1229, May 1958, 17 pp. 

High-pressure Gas-to-liquid Trans- 
ducer, by Don O. Coffin, Rev. Sci. Instr., 
vol. 29, Oct. 1958, p. 896. 

Microwave Manometer, by Arnold G. 
Kramer and Philip M. Platzman, Rev. 
Sci. Instr., vol 29, Oct. 1958, p. 897. 


Guidance Components and 
Systems 


How to Beat Instrument Blindness, by 
George A. LeHew and James R. Andres, 
Aviation Age, vol. 30, Oct. 1958, pp. 64-69. 

Microelectronics: Small Circuits Set for 
Big Role in Future Systems, by James 
Holahan, Aviation Age, vol. 30, Oct. 1958, 
pp. 20-21, 166-170. 

Guidance for the Space Age, by C. J. 
Mundo and I. W. ‘tae. Aero/Space 
Engng., vol. 17, no. 9, Sept. 1958, pp. 
30-33. 

Some Effects of Wind and Temperature 
Gradients on the Design of Missile 
Flight Control Systems, by Melvin W. 
Lifson, Aero/Space Engng., vol 17, no. 
9, Sept. 1958, pp. 49-52. 

Effects of Gyrosconic Cross Coupling on 
Missile Automatic Flight Control Systems, 
by J. G. MeCracken, Aero/Space Hngng., 
vol. 17, Nov. 1958, pp. 70-73. 


Flight Vehicle Design and 
Testing 


On the Optimum Trajectory of a 
Rocket, by Robert R. Newton, J. Franklin 
Inst., vol. 266, Sept. 1958, pp. 155-187. 

Topics in Dynamic Programming for 
Rockets, by Angelo Miele and James O. 
Cappellari Jr., Purdue Univ., School of 
Aeron. Engng., Rep. A-58-5 (AFOSR-TN- 
58-685; ASTIA AD 162218), July 1958, 
58 pp., 7 figs. 

Missile Noise as a Factor in Reliability, 
by Paul S. Veneklasen, Aero/Space 
Engng., vol. 17, no. 9, Sept. 1958, pp. 
44-48. 

Problems in the Development of a 
Guided Weapon, by J. Clemow, J. Roy. 
Aeron. Soc., vol. 62, no. 575, Sept., 1958, 
pp. 663-673. 

On Flight Trajectories in the Neighbor- 
hood of a Known Trajectory, by R. M. 
Rosenberg, J. Franklin Inst., vol. 266, no. 
2, Aug. 1958, pp. 109-128. 


Space Flight 


Rocket Propulsion Requirements for 
Satellites and Moon Missions, by John L 
Sloop, Aero/Space Engng., vol. 17, Nov. 
1958, pp. 36-44, 50. 

An Application of Solar Radiation to 
Space Navigation, by S. James Press, 
Aero/Space Engng., vol. 17, Nov. 1958, 
pp. 51-54. 
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TAS. 


His full name is Yuan Chu’en Lee, but to most of us at 
Aerojet, he’s known as “Y.C.” As Director of Corporate 
Research, Y.C. is in charge of our astronautical, other- 
planetary endeavors. 


Y.C. and his associates are grappling with the knottier 
questions of space technology and advanced propulsion. 
They’re studying advanced magnetohydrodynamics, 
extraterrestrial physics, and relativistic problems. They’re 
looking for new fuels based on ions, gaseous free radicals, 
and nuclear plasmas. 


AZUSA AND SACRAMENTO, CALIFORNIA + 


AEROJET-GENERAL CORP.’ 


SUBSIDIARY OF THE GENERAL TIRE & RUBBER COMPANY 


Lee 


Fundamental research has kept Aerojet at the forefront 
of the rocket field. Thanks to people like Y.C. and his 
team, we'll stay there. 


That team, by the way, is hand-picked. It’s one of 
several such groups at Aerojet that offer opportunities of 
unusual scope to the scientific and engineering imagination. 
Let us see a resume if this is your kind of challenge. The 
address: AEROJET-GENERAL CORPORATION, P.O. Box 296, 
Azusa, California or P.O. Box 1947, Sacramento, California. 
Attn: Director of Scientific and Engineering Personnel. 
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Air-breathing Power Plants in the 
Space Era, by Peter G. Kappus, Aero/- 
Space Engng., vol. 17, Nov. 1958, pp. 
62-65, 69. 

On the Mechanisms of Meteoric Abla- 
tion, by Sin-I. Cheng, Princeton Univ., 
Dept. Aeron. Engng., Rep. 427 (AFOSR- 
71'N-58-649 ; ASTIA AD 162181), April 
1958, 21 pp., 8 figs. 

Bibliography of Space Medicine, Public 
Health Service, Publication 617 (Bibliog- 
raphy series 21), 1958, 49 pp. 

Satellite Orbits in Theory and Practice, 
by D. G. King-Hele and R. H. Merson, 
/. Brit. Interplanet. Soc., vol. 16, no. 8, 
July-Aug. 1958, pp. 446-453. 


Russian Technical 


Articles* 
About Noise Stability of Pulse 
Frequency Telemetering, by N. V. Pozin, 


{utomatika 7 Telemekhanika, vol. 19, no. 
10, Oct. 1958, pp. 968-979. (In Russian.) 

Pulse Servosystems Including Two Pulse 
Elements with Unequal Repetition Periods, 
by Fan Chun-wui, Automatika i Tele- 
mekhanika, vol. 19, no. 10, Oct. 1958, pp. 
917-930. (In Russian.) 

Determination of Periodic Behavior of 
Automatic Control Systems Having a Non- 
linear Part with Piece-linear Character- 
istic, by L. A. Gusev, Automatika 7 
Telemekhanika, vol 19, no. 10, Oct. 1958, 
pp. 931-944. (In Russian.) 

About Equivalency of Discontinuous and 
Continuous Control Systems, by V. A. 
Rubtsov, Automatika i Telemekhanika, 
vol. 19, no. 10, Oct. 1958, pp. 945-957. 
(In Russian.) 

On the Theory of Gas Jets, by L. N- 
Stretensky, Doklady Akad. Nauk SSSR; 
vol. 119, April 21, 1958, pp. 1113-1114- 
(In Russian.) 

A Spectrographic Investigation of the 
Condition of a Gas Behind a Shock Wave, 
by N. N. Sobolev and others, /zvestiia 
Akad. Nauk SSSR, Ser. Fiz., vol. 22, June 
1958, pp. 730-736. (In Russian.) 

On Certain Exact Solutions of the Gas 
Dynamics Equations, by S. 8. Grigorian, 
(translated from Doklady Akad. Nauk 
SSSR, vol. 121, Aug. 1, 1958, pp. 606- 
609). Morris D. Friedman, Ine., Need- 
ham Heights 94, Mass., Translation no. 
G-156. 

One-Dimensional Self-Similar Motions 
of a Conducting Gas in a Magnetic Field, 
by V. P. Korobeinikov (translated from 
Doklady Akad. Nauk SSSR, vol. 121, 
Aug. 1, 1958, pp. 613-615). Morris D. 
Friedman, Inc., Needham Heights 92, 
Mass., Translation no. K-187. 

Mass Transfer During Absorption with 
Bubbling, by A. A. Ravdel (translated 
from Zhur. Priklad. Khim., vol. 31, no. 7, 
1958). Associated Tech. Services, P. O. 
Box 271, East Orange, N. J., Translation 
no. RJ-1500. 

The Solution of Gas Bubbles in a Liquid, 
by A. A. Ravdel (translated from Zhur. 
Priklad. Khim., vol. 26, no. 7, 1953, pp. 


* The editors of Technical Literature 
Digest are making a systematic search for 
pertinent Russian articles, both in the orig- 
inal Russian and in translated form. For 
the balance of 1959, these will appear in a 
separate section with this heading, mainly 
to draw the readers’ attention to the new 
listing. After that, they will appear item 
by item under the proper subject headings. 
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752-756). Associated Tech. Services, P. 
O. Box 271, East Orange, N. J., Transla- 
tion no. RJ-1501. 


The Use of Frontal Analysis in Gaseous- 
Liquid Chromatography of Radioactive 
and Non-radioactive Gases, by M. I. 
Ianovsky and G. A. Gaziev, Doklady 
Akad. Nauk SSSR, vol. 120, June 1, 1958, 
pp. 812-814. (In Russian. ) 


Chemical Actions of Atoms Having En- 
ergies Comparable with the Activation 
Energy, by B. N. Provotorov, Doklady 
Akad. Nauk SSSR, vol. 120, June 1, 1958, 
pp. 838-840. (In Russian.) 


Quasi-stationary Thermal Regime of 
Explosive Reaction Processes, by A. G. 
Merzhanov and F. I. Dubovitsky, Dok- 
lady Akad. Nauk SSSR, vol. 120, June 3 4 
1958, pp. 1068-1071. (In Russian.) 


A Boundary Problem in the Thermal 
Explosion Theory, by V. V. Barzykin and 

. G. Merzhanov, Doklady Akad. Nauk 
SSSR, vol. 120, July 21, 1958, pp. 1271- 
1273. (In Russian. ) 


On the Theory of Thermal Chain Flame 


Propagation, by L. A. Lovachev, Doklady 
Akad. Nauk SSSR, vol. 120, July 21 
1958, pp. 1291-1293. (In Russian.) 


Relaxation Reaction Between Shock 


Waves and the Burning Zone, by S. M. 
Kogarko and V. I. Skobelkin, 
Akad. Nauk SSSR, vol. 120, July 21, 
1958, pp. 1280-1283. (In Russian). 


On the Stability of Flame Propagation, 
by G. I. Barenblatt and Ia. B. Zeldovich 
(translated from Priklad. Matemat. 1 
Mekhan., vol. 21, no. 6, 1957, pp. 856-859). 
Morris D. Friedman Inc., Needham 
Heights 92, Mass., Translation no. B-135. 


Exact Solution of the Nonlinear Prob- 
lem of an Explosion in a Gas with Var- 
iable Initial Density, by V. P. Korobeini- 
kov (translated from Doklady Akad. Nauk 
SSSR, vol. 117, no. 6, 1957, pp. 947-948). 
Morris D. Friedman, Ine., 
Heights 92, Mass., Translation no. K-170. 


On Spin Near the Gas Detonation 
Limits, by Ia. K. Troshin and K. I. 
Shchelkin (translated from [zvestiia Akad. 
Nauk SSSR, Otdel. Tekh. Nauk, no. 8, 
1957, pp. 142-143). Morris D. Friedman, 
Inc., Needham Heights 92, Mass. Trans- 
lation no. T-129. 


On Spin Detonation, by B. V. Voitse- 
khovskii (translated from Doklady Akad. 
Nauk SSSR, vol. 114, no. 4, 1957, pp. 
717-720). Morris D. Friedman, Inc., 
Needham Heights 92, Mass., Translation 
no. V-125. 

Errors in the Measurement of Flame 
Temperature in a Stream with Thermo- 
couples, by K. P. Vlasov and N. V. Koku- 
shkin (translated from Jzvestiia Akad. 
Nauk SSSR, Otdel. Tekh. Nauk, no. 8, 
1957, pp. 137-141). Morris D. Friedman, 
Inc., 
lation no. V-128. 


Investigation of Certain Singularities in 


the Combustion of Motor Fuel Sprays, by 
D. I. Raibov and Iu. B. Sviridov, Akad. 
Nauk SSSR, fzvestiia, Otdel. Tekh. Nauk, 
4, April 1958, pp. 133-137. (In Russian. ) 


Some Instances of Motions of a Viscous 
Gas in Narrow Spaces of Variable Thick- 
ness, by Ya. M. Kotlyar, Akademii Nauk, 
SSSR, Izvestita, Otdel. Tekh. Nauk, 5, 
May 1958, pp. 35-39. (In Russian.) 


Effect of a Slightly Blunt Leading Edge 
upon its Streamlined Flow at Hypersonic 
as by G. G. Chernyi, Akademii 

uk SSSR, Izvestiia, Otdel. Tekh. Nauk, 
no. 4, April 1958, p. 54-66. (In Russian.) 


‘Doklady | 


Needham | 


Needham Heights 92, Mass., Trans- | 


Environment 
for Dynamic 


Career Growth 
for 
Engineers and Scientists 


Ability is important in getting ahead in 
engineering and science. But of almost 
equal importance is the environment 
in which you work. A dynamic environ- 
ment that provides challenge and scope, 
unexcelled facilities, and plenty of room 
at the top, can accelerate your advance 
to higher levels of responsibility and 
income. 

This is the kind of environment you'll 
find at Boeing...pictured and described 
in Boeing’s new 24-page booklet, “En- 
vironment for Dynamic Career Growth.” 

The booklet, in addition, reports on 
engineering and science assignments in 
connection with current Boeing proj- 
ects—from advanced military and com- 
mercial jet aircraft to space vehicles, 
guided missiles and intercontinental 
ballistic missiles. It also outlines Boeing 
research and development activities, 
and documents the dynamic Boeing 
environment that fosters career growth 
of engineers and scientists. 

There are openings at Boeing, now, 
for engineers and scientists of all cate- 
gories, all experience and educational 
levels. A Boeing assignment in Research, 
Design, Production or Service could be 


the answer to your future. 


Write today for your free 
copy of the 24-page book 
‘Environment for Dy- 
namic Career Growth.” 
Indicate your degree(s) 
and field of interest. Ad- 
dress: Mr. Stanley M. 
Little, Boeing Airplane 
Co.,P.O. Box 3822 - JPB 


Seattle 24, Washington. 
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CHARACTERISTICS BoEING AIRPLANE Co 
ANALYSIS RECOMMENDED USE Calkins & Holden, Seattle, Wash. 
i For $s operating under high 
Stainless Steel Ball and Race { soln ts (800.1200 degrees F.). 
Chrome Alloy Steel Ball For types operating under high radial CALLERY CHEMICAL Co 
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neers welcome an opportunity of studying individual requirements 
and prescribing a type or types which will serve under your demand- E. I. pv Pont pg Nemours & Co., Inc 
ing conditions. Southwest can design special types to fit individual N.-W. Ayer & Son, Inc., Philadelphia, Pe. 
specifications. As a result of thorough study of different operating 
conditions, various steel alloys have been used to meet specific 
needs. Write for Engineering Manual No. 551.Address Dept. JP-59. GenerAL Exectric Co., ArrcrAFT NUCLEAR 
PROPULSION DEpT 


SOUTHWEST PRODUCTS CO. Deutsch & Shea, Inc., New York, N.Y. 
1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


GENERAL Motors Corp., OLDSMOBILE Div.......... 
D. P. Brother & Co., Detroit, Mich. 


CH EM ISTS GrovE VALVE & REGULATOR Co 3rd cover 


ANALYTICAL PhD or equivalent to perform research on new L. C. Cole & Co., Inc., San Francisco, Calif. 
analytical methods and techniques. 
PhD or equivalent synthetic 

istry. Industrial experience required—prefera in fuels, 
oxidizers and propellant technology. INTERNATIONAL BusINESS MACHINES Corp... ... . .388, 389 
PHYSICAL—PhD or M.S. SPeiene © experience and back- Benton & Bowles, Inc., New York, N.Y. 
ground in exp and 
cybernetics required. Te supervise research programs in 
—- of physical and ballistic properties of liquid and solid 
propellants. 
PHYSICAL—BS. with education or ex LockHEED MissiLEs & SPACE Div 380, 381 
experimental thermodynamics, ~ tome g capacity Hal Stebbins, Inc., Los Angeles, Calif. 
measurements. 


VanS , Md. 
DEVELOPMENT ENGINEERS with experience in design anSant Dugdale & Co., Inc., Baltimore, Md 
analysis, and operation of bench scale units and pilot plants. 
SENIOR PROCESS ENGINEER with broad rather than 


specialized pilot plant and di — er 
and economic studies. Technical supervisory experience Newsrook Macuine Corp 
Melvin F. Hall Agency, Inc., Buffalo, N.Y. 


THER RaMo-Woo.pRipGE, A Div. or THOMPSON 


APPLICATIONS RESEARCH CHEMIST or ENGINEER with AMo-WooL DGE INc 

industrial experience and interest in several technological fields. Ramo-WooLDRIDG . x 

To explore new industrial products or applications of existing The McCarty Co., Los Angeles, Calif. 

products. 

PROPELLANT TESTING SPECIALIST—Chemist, Physicist 

or Engineer, B.S., M.S., or PhD with a minimum of 3 years’ 

pred in small scale testing and evaluation of liquid and RepuBuic AVIATION Corp 

programs for Deutsch & Shea, Inc., New York, N.Y. 

Permanent posit: with prod: of Boron chemicals and high 
energy fuels located 25 miles north of Pittsbargh’s Golden Triangle. 


Gd. 


Please send resume including salary req ts in to: 


PERSONNEL MANAGER 


LLERPY Space TECHNOLOGY LABORATORIES, INC............. 
Callery Gaynor & Ducas, Inc., Beverly Hills, Calif. 
CHEMICAL COMPANY Pennsylvania 


SouTHwEst Propucts Co 
O. K. Fagan Adv. Agency, Los Angeles, Calif. 


ARS JouRNAL 


INT. 
PLAIN TYPES ROD END 
TYPES 
ge 
* | 
€ 
q 888 
s 
ee 
|| 
i} 
e 
. e 
2 
| 
| 
: | fer 
| 
392 


a 


